A SPECT reconstruction method for extending parallel to non-parallel geometries Junhai Wen and Zhengrong Liang
|
|
- Brook Matthews
- 5 years ago
- Views:
Transcription
1 This content has been downloaded from IOPscience. Please scroll down to see the full text. Download details: IP Address: This content was downloaded on 03/04/2019 at 09:09 Please note that terms and conditions apply. You may also be interested in: Measurement and Instrumentation Principles Alan S Morris Design and instrumentation of a magnetic nerve stimulator D McRobbie Simple boxcar integrator covering a large input frequency range P Williams and H W Lam New books on physics and related sciences Elena V Zakharova New books on physics and related sciences Elena V Zakharova A SPECT reconstruction method for extending parallel to non-parallel geometries Junhai Wen and Zhengrong Liang Voigt wave propagation in biaxial composite materials Tom G Mackay and Akhlesh Lakhtakia New books on physics and related sciences Elena V Zakharova
2 IOP Publishing Orthogonal tensors in computational engineering mechanics R M Brannon Chapter 1 Introduction Each practicing engineer must, at some point, characterize reorientations of a body or find vector and tensor components with respect to a rotated coordinate system. Additionally, rotations and reflections can contribute to overall engineering system response in combination with other predominant physical effects. Consequently, many important theorems (such as the polar decomposition and the time derivative of logarithmic strain) involve rotations as a part of a larger problem. A vast literature is available for the most common rotation-related tasks such as simple orthogonal Cartesian coordinate changes. Most reference books, however, tend to cover one or two methods for handling rotations, which makes it difficult for a researcher to translate between different methods used in different books. Resource information for less common rotation-related tasks is even scarcer. For example, mesoscale material simulations (in which crystals are explicitly modeled in a representative volume element) require methods to establish a uniform distribution of initial crystal orientations, which is a non-trivial task that is rarely explained in textbooks. Obscure rotation-related research applications can be found in disparate journal articles, but self-contained comprehensive review articles that cover both elementary and advanced engineering concepts written in a manner comprehensible to engineers are rare. Morawiec [1] offers an outstanding list of references on the topic of rotation and covers many of the topics reviewed here with greater mathematical rigor. Available rotation review articles and textbooks seem to be aimed more at mathematicians than applications specialists. While it is true that this book will cover several mathematically advanced topics, our target audience is practicing engineers, especially those who lack a guru-level background in abstract algebra, topology, differential geometry, hyper-complex numbers, and group theory, that unfortunately seems to be a prerequisite for understanding existing resource publications. Hence, most of our discussions are limited to ordinary three-dimensional (3D) vectors of the type discussed in first calculus courses. Detailed proofs are provided only when they are non-trivial (and not readily available in the literature) doi: / ch1 1-1 ª IOP Publishing Ltd 2018
3 or when the proof itself clarifies a concept. In this book, the assumed level of the reader s mathematical background varies considerably. Some chapters are quite elementary (first-year graduate student level), while other chapters are aimed at a target audience of mathematically adept researchers who are working on the arcane topics of contemporary interest. This book s topics and examples tend to favor materials mechanics issues. Rather than organizing this book using increasing levels of complexity, the topics are grouped according to subject area. Each chapter (and each section within a chapter) opens with an elementary (more intuitive) approach to the subject at hand, and then it progresses into greater levels of sophistication and rigor. Consequently, if you have only limited background in mathematics, you are encouraged to drill into each chapter only to the depth appropriate to your background. Each chapter is (as much as practical) self-contained, so it is not necessary to read preceding chapters to learn about subjects in later chapters. One aim of this book is to catalog the many different ways to define and quantify rotations and reflections. A rotation may be described mathematically by a special kind of orthogonal tensor. By using the special kind proviso, we are implying that rotation tensors are always orthogonal, but orthogonal tensors are not necessarily rotations. A tensor Q is orthogonal if and only if T Q Q = I, where Ĩ is the identity tensor 1. The tensor Q is orthogonal if and only if its component matrix [Q] is orthogonal. A real matrix [Q] is orthogonal if and only if where [I ] is the identity matrix. In indicial form, T [ Q] [ Q] = [ I], (1.1) QkiQ kj = δij (1.2) where δ ij is the Kronecker delta 2 and (using standard tensor notation conventions summarized in chapter 2) repeated indices are understood to be summed from 1 to 3. An elementary identity in tensor analysis states that a tensor is zero iff3 Z u Z v = 0 vectors u and v (1.3) Accordingly, two tensors, Ã and B, are equal iff u A v = u B v vectors ũ and v. A tensor Q is therefore orthogonal iff u Q T Q v = u I v vectors u and v (1.4) 1 This book uses a self-defining notation in which scalars, vectors and second-order tensors are written with zero, one, and two under-tildes, respectively (as in s, ṽ, and ). The composition of two tensors, denoted symbolically by A T B, has ij components given by AikBkj in which the repeated index is implicitly summed. The notation AB : = AB ij ij is the tensor inner product. The notations A u and v B are vectors with ith components of, respectively, A u ik k and v k B ki. These definitions immediately give the often-used identity A u = u A T. 2 I.e., components of the 3 3 identity matrix. Thus δ ij equals 1 if i = j and 0 if i j. 3 I.e., if and only if. 1-2
4 To interpret this assertion, consider the star vectors in figure 1.1, which are defined to be orthogonal transformations of non-starred fiducial vectors: u* = Q u and v* = Q v. (1.5) Noting that T u Q = Q u, equation (1.4) then reduces to u * v*= u v vectors u and v. (1.6) As a special case, taking the two arbitrary vectors to be identical to each other proves that an orthogonal Q must preserve lengths u * = * u and v = v. (1.7) Denoting these magnitudes by u and v, respectively, recall from elementary vector analysis that a vector dot product satisfies u θ v = uv cos, where θ is the angle between the vectors. Consequently, equation (1.6) shows that an orthogonal transformation Q preserves the magnitude but not necessarily the sign of the angle between two vectors: cos θ*= cos θ θ*=± θ (1.8) In two dimensions, the angle θ between ũ and ṽ can be arbitrarily declared to be positive if the vector ũ must be rotated counter-clockwise by θ to become aligned with ṽ. As indicated in figure 1.1, a 2D orthogonal transformation is called proper if it preserves the angle in both magnitude and sign; it is otherwise an improper transformation that will preserve the magnitude of the angle, but not the sign. An improper-orthogonal transformation is seen in the example of figure 1.1 to not only rotate the vectors but also to reflect them into different relative ordering (from clockwise to counter-clockwise or vice versa). These statements may be generalized to higher dimensions based on properties of the [Q] orthogonal transformation matrix itself. Figure 1.1. Effects of an orthogonal transformation. (a) The word fiducial simply means basis of comparison, which refers to the pre-transformation starting vectors (b) shows a change in orientation without a change in length or angle and (c) shows a change in both orientation and numerical sign of the angle. 1-3
5 It is fairly common for writers to refer to the set of all tensors as lin, and the subset of all orthogonal tensors is called orth [2]. Second-order engineering tensors have nine independent components, so they are said to belong to a 9D space. Moreover, because any linear combination of tensors is itself a tensor, lin is a 9D linear manifold and therefore it is possible to set up a Cartesian coordinate system for all members of lin. To rigorously prove that a second-order tensor on 3D space is, in fact, also a 9D vector, one must simply confirm axioms (additive commutativity, associativity, existence of an additive inverse, etc) that must be met for something to be called a vector. Just as the Cartesian representation for a vector x = x1e + e + e 1 x2 x (in which the ẽi are orthonormal basis vectors) may be represented in array form as x x 2 = x1 0 + x2 1 + x3 0, (1.9) x the Cartesian expansion of a tensor, T = T11e1e1+ T12e1e2+ + T33e e 3 3 (in which the ee i j are basis dyads), may be represented in matrix form as T 11 T12 T T 21 T22 T23 = T T T (1.10) T T T This expansion has nine terms involving nine orthonormal basis tensors. The space of ordinary 3D vectors is Euclidian (and thus has a Cartesian basis), but it can contain lower-dimensional spaces, such as the surface of a sphere, that are not Euclidean which means that a linear combination of points in that curvilinear space might not produce a point in the space. A point on the unit sphere, for example, is a vector of unit magnitude. Adding unit vectors does not generally produce a unit vector, so the set of all unit vectors is non-euclidean. This does not mean they are not vectors. They simply belong to a 2D curvilinear space contained within the larger 3D Euclidean space of general vectors. Similarly, we will find that the set of all orthogonal tensors belongs to a 3D curvilinear space (a hyper-sphere) embedded within the larger 9D Euclidean space of general second-order tensors. Unit vectors have only two independent components (given any two components, the third is constrained to produce a unit vector). Even though the components of a unit vector are not independent, it is nevertheless most intuitive to work with them using all three components. Similarly, even though rotations have only three independent components, it is nevertheless most intuitive to work with them in their full 3 3 form. A reduction to fewer components is useful mostly to conserve computer memory, but typically at the penalty of increased CPU cost. The fact that orth is made of tensors having only three independent components follows because columns (or rows) of an orthogonal matrix must form an orthonormal triad. The first column has two independent components (because it is a unit vector), and the second column has only one independent component because it must 1-4
6 be a unit vector and also have a zero inner product with the first column. Once the first and second columns of a rotation matrix are set, the third is the cross product of the first two, giving no additional independent components (other than sign if reflections are allowed). Because a linear combination of unit vectors is not generally a unit vector, it follows that a linear combination of orth tensors is not generally an orth tensor itself. Consequently, orth is a 3D nonlinear (curvilinear) subset of lin. As a nonlinear space, orth is not a flat higher-dimensional plane; instead, it will be seen to belong to a subset of a hyper-sphere in lin. Before doing any high-math to prove this, you should be able to guess the radius of this sphere: just as the magnitude of a vector v is vv k k, the magnitude of a tensor à is AA ij ij. If the tensor happens to be orthogonal Q, its magnitude (and hence the radius of the hyper-sphere) is T QQ = Q Q = δ = 3, (1.11) ij ij ji ij which applies if [Q] isa3 3 orthogonal matrix. In general, an N N orthogonal matrix has magnitude N. An orthogonal matrix [Q] with a positive determinant (equal to +1) is called a rotation or proper orthogonal matrix [3]. The set of all proper-orthogonal tensors is often referred to as orth +. An orthogonal matrix [Q] with a negative determinant (equal to 1) is said to be an improper-orthogonal tensor and is often misleadingly called reflection but it may actually be a reflection in combination with a rotation (see chapter 6). The set of improper-orthogonal tensors is often referred to as orth. A proper-orthogonal tensor is the most useful means of characterizing a rotation even though it is inefficient to store as nine components when only three are independent. Euler angles are also popular, but of limited value. These and many other lesser-known representations will be reviewed. Orthogonal tensors emerge as the superior quantification of a rotation or reflection 4, so each definition of an alternative representation of a rotation will be accompanied by instructions for converting from that representation to an orthogonal tensor and vice versa. Algorithms and source code are provided (and may be downloaded from In addition to cataloging the many ways to quantify rotation, this book covers several rotation-related applications. The topics are arranged as summarized here. For a smaller list that highlights the most significant contributions of this text, see the closing remarks in chapter 28. ij 4 Having nine components, an orthogonal tensor is not superior to other methods if minimizing storage is the metric. Instead, it is the structure and ease of computation that tends to make orthogonal tensors convenient. In the same way, a unit vector ñ is still most conveniently described in terms of three components ( n 1, n 2, n 3 ) satisfying one constraint, n1 2 + n2 2 + n3 2 = 1, even though they may be described in terms of two spherical angles as n1 = cos φ sin ϑ, n2 = sin φ sin ϑ, n3 = cos ϑ. This representation, though efficient in storage, has an ugly imbalance as well as numerical inefficiency in rebuilding the components. Similar assertions apply to rotations. The most compelling advantage of building a full three-component unit vector is that the set of unit vectors forms a 2D curvilinear subset (surface of a sphere) contained within a larger 3D Euclidean space. Euclidean representations are needed to compute Jacobians (needed to quantify cluster density of unit normals on a sphere or, in a higher-dimensional analog, texture of crystallographic orientations in materials). 1-5
7 Chapter 2 outlines notation and key mathematical prerequisites, as well as explaining the linear fractional transformation notation commonly used in robotics/controls. Chapter 3 reviews how orthogonal direction-cosine matrices are used to transform vector and tensor components from one orthonormal basis to another, and how these transformations are used to define engineering vectors or tensors as closure tests in the abstract mathematician s definitions that make no reference to length or direction. The distinction between a position vector s components and coordinates is made clear. Chapter 4 points out the distinction between rotation and coordinate transformation. Analyzing rotation demands awareness of your desired perspective. You can rotate an object, while you stay still, or you can keep the object fixed while you rotate yourself. It is important to be aware of which of these perspectives applies for your problem of interest. The distinction between these fundamentally different transformations goes beyond one being the same as the other with an opposite rotation angle. Chapter 5 describes nonlinear rotations (using torsion and vortex as examples), how to convert an axis and angle of rotation into a rotation tensor and vice versa, properties of skew tensors associated with rotation axes, and spinor vector representations of rotations. Chapter 6 reiterates that a rotation is characterized by a proper-orthogonal tensor i.e., one having a determinant equal to +1. An orthogonal tensor with a determinant equal to 1 should be regarded (in general) as a reflection, possibly in combination with a rotation. Chapter 7 presents the representation of a rotation in terms of a unit quaternion. Quaternions, which (being generalizations of complex numbers) are predecessors to modern vectors, and they are enjoying a resurgent vogue in rotationrelated applications. Even though a rotation tensor can be constructed uniquely from a unit quaternion (having one real component and three hyper-imaginary parts), it is explained that this does not imply that the smallest Euclidean space containing all possible rotation tensors is itself 4D. Chapter 8 describes how any linear vector-to-vector operator (in 3D) can be expressed as a sum of three dyads, and this result is specialized to rotations. This chapter provides yet another way to describe rotations that is based on a colorful nautical analogy between embedded triads in an object and embedded directions in a sailboat (mast, starboard, bow). Any orientation of these embedded directions can be specified through the boat s latitude, longitude and heading. Applications of dyadic representations of tensors include user interfaces for controlling rotation of a 3D object on a computer image, finding basis orientations that most simplify a tensor s component matrix, and constructing a rotation that will transform a specified vector into another specified vector. Chapter 9 shows that any general rotation may be expressed as a (noncommuting) sequence of three rotation operations. When defining a rotation in this way, it is essential to know whether subsequent rotations are defined with respect to fixed axes or follower axes. Rotations applied sequentially 1-6
8 about the fixed laboratory basis are different from rotations applied sequentially about axes of a triad that moves with the body as quantified by Euler angles, but the two descriptions are intimately related. Chapter 10 provides the series expansion representation of a rotation, which turns out to be a tensor analog to the series expansion of an exponential. This chapter also describes the Cayley transformation that constructs a rotation tensor from a skew tensor, as well as the Rayleigh transformation (which is used to apply a superimposed rotation to a tensor of any order). Chapter 11 derives the eigenvalues and eigenvectors of any rotation tensor. Chapter 12 presents the polar-decomposition theorem, which shows how a general deformation can be broken into two distinct (non-commuting) steps, one being a rigid rotation, and the other embodying pure material distortion and/or volume changes with no net material rotation. Chapter 13 describes how the polar decomposition can be used to strip away overall material rotation from material deformation, quantified by various strain measures commonly used in engineering applications. The spatial stretch which generates spatial strain measures is altered in the presence of rotation. The reference stretch which generates reference strain measures is unaffected by rotation. Chapter 14 discusses candidate strategies for advecting rotations in Eulerian physics codes or in Lagrangian codes that remap the mesh. The issues are similar to the challenge of mixing or interpolating unit vectors. Chapter 15 shows how the rate of a rotation tensor is related to the conventional angular rotation vector, vorticity tensor, and vorticity vector. This chapter includes derivatives of a rotation with respect to the rotation angle and with respect to the rotation axis as well as numerical algorithms for accurately integrating a rate of rotation to update material orientation in a physics code. The famous Dienes algorithm for computing polar spin from vorticity is included, along with a debunking of the common (but erroneous) assertion that eigenvectors of the symmetric part of the velocity gradient rotate according to the vorticity tensor defined by the skew part of the velocity gradient. Chapter 16 shows how to compute the time rate (or other derivative) of a principal function of a symmetric tensor, with emphasis on finding the time rate of logarithmic strain. A principal function (such as the Hencky logarithmic strain used in materials modeling) requires casting the argument (a symmetric tensor) in its principal basis to apply the function to the eigenvalues, after which the result is cast back to the physical basis. Finding the rate of such functions is crucial for linearizing numerical solvers, yet is extraordinarily complicated by the fact that principal directions not only evolve but are not even unique when there are repeated eigenvalues. An algorithm for taking rates of principal functions is presented and shown to be valid even when eigenvalues are repeated. The results are furthermore used to express the rate of reference Hencky strain as a linear transform of un-rotated rate of deformation (defined by a Rayleigh un-rotation of the symmetric 1-7
9 part of the velocity gradient and proved in this chapter to not be a pathindependent rate of any measure of deformation). Chapter 17 shows how to generate a uniformly random rotation tensor, which is useful for generating grain orientations for microscale simulations. We also discuss how to find the average of a tensor over all possible (uniform) rotations of that tensor. 5 This chapter includes extensive tutorials on elementary statistical sampling methods, as well as clear explanations about how to construct a non-square Jacobian matrix for transforming random numbers to random realizations and to then use it to prove whether or not an algorithm produces a uniform sampling. A case study, for example, proves that uniform sampling of Euler angles will not produce a uniformly random rotation, as is needed (for example) in mesoscale polycrystalline constitutive modeling. Chapter 18 defines isotropic tensors, which are very special tensors whose components do not change upon an orthogonal change in basis. A distinction is made between proper isotropy and strict isotropy, both of which are useful concepts. The set of all isotropic tensors of a given class 6 is a linear manifold, which means that it must possess a finite number of isotropic tensors (called primitive isotropic tensors) that form a basis for all isotropic tensors of that class. In three dimensions, there is only one primitive isotropic second-order tensor the identity tensor. There are two primitive proper-isotropic tensors in two dimensions. Fourth-order tensors in three dimensions have three primitive isotropic tensors. Chapter 19 provides an elementary introduction to the principle of materialframe indifference (PMFI), which requires that material constitutive models must be invariant under rigid rotation. An important point is that there do exist scalar invariants (which, by definition, will not change under a basis rotation) that will change upon a change in frame (i.e., a superimposed rotation). Likewise, a vector or tensor never changes under a basis rotation (the components change, but the sum of components times basis vectors or basis dyads remains invariant); whether or not a vector or tensor will change under a superimposed rotation depends on its physical definition. The discussion of PMFI includes a review of objective rates of spatial tensors, which are often employed in materials modeling (though we argue against this approach in favor of genuine rates within un-distorted or un-rotated material reference frames). Chapter 20 explains the difference between tensor symmetry and material symmetry, and then proceeds to prove that isotropic tensors are not always limited to mere multiples of the identity. Transverse symmetry (for which a tensor s components are unchanged upon rotation about a fixed axis) is covered in detail, including a rather brute force (hence accessible to a broader audience) explanation of why double-symmetric fourth-order transverse 5 Spoiler alert: it is proved to be the isotropic part of that tensor. 6 Where class refers to a tensor order (also called rank ) and the dimension of the underlying physical space. 1-8
10 stiffness tensors have five independent components, but only four distinct eigenvalues. The corresponding eigenprojectors are derived and a closed-form solution is derived for necessary and sufficient conditions for two such tensors to have the same eigenspaces (useful, for example, to force elastic and plastic stiffness tensors to have the same spectral symmetries). Application to jointed rock is also discussed. Chapter 21 defines scalar invariants of tensors, which are scalar-valued properties of tensors whose values do not change upon a change of basis. Not all real numbers are scalars. Chapter 22 returns to the topic of frame indifference as it applies to incremental constitutive models, first explaining the effect of rotation rate using an easily understood spring model, then moving towards various co-rotational rates in the literature, finishing with the obligatory derivation of oscillatory stress using certain rates (and explaining why it probably is not the rate that is the core issue) along with some guidance about how to run frame indifference testing in computational mechanics. Chapter 23 specializes the general continuum laws of mechanics to rigid-body mechanics. By applying the full power of tensor analysis, this presentation distinguishes itself from the derivations normally found in elementary dynamics textbooks. For example, with direct-notation tensor analysis, Euler s equations and the parallel-axis theorem can be presented in a single line. Chapter 24 returns to the problem of deformable media, but considers the case that bulk rigid motion is the primary motion, with relatively small deformations with respect to an appropriately defined rigidly moving observer frame (e.g., turbine blade, human bodies in a braking automobile, etc). Phantom body forces are derived to allow solving the problem in a non-inertial material frame. Chapter 25 covers some basic equations for 3D computer graphics visualization, in particular providing equations governing mapping from a virtual 3D body and providing recommendations for 3D visualization that would not require a person to wear 3D goggles. Chapter 26 defines conventional Voigt components for arranging the six independent components of symmetric tensors into a 6 1 array and demonstrates that these are in fact components of a 6D vector whose associated basis is not normalized (which is why manipulations using Voigt components always contain unruly factors of 2 that are actually the metrics associated with the non-normalized Voigt basis). Normalization of the Voigt basis produces the Mandel basis with associated Mandel components that are not plagued by factors of 2 and which therefore obey standard linear algebra, permitting more efficient use of computer library functions such as eigenvalue solvers for fourth-order tensors. Using Mandel components, the 6 6 component matrix associated with a tensor rotation operation in physical space is derived. 1-9
11 Chapter 27 completes the main text by providing a brief discussion of rotations in higher dimensions, which is applicable to constitutive material modeling. Ordinary physical tensor rotation operations are shown to correspond to double-plane rotations in stress space, where part of the tensor rotates by an amount equal to the rotation angle, another part rotates by twice this angle (as in Mohr diagrams), and a third part of the tensor does not rotate at all in tensor space. Appendix A provides FORTRAN and Python source listings that perform most of the computations presented in this book. Additional source code in these and other programming languages are available. References [1] Moraweic A 2004 Orientations and Rotations: Computations in Crystallographic Textures (Berlin: Springer) [2] Bigoni D 2012 Nonlinear Solid Mechanics (New York: Cambridge University Press) [3] Gurtin M E 1981 An Introduction to Continuum Mechanics (New York: Academic) 1-10
Here, the parameter is the cylindrical coordinate, (a constant) is the radius, and the constant controls the rate of climb.
Line integrals in 3D A line integral in space is often written in the form f( ) d. C To evaluate line integrals, you must describe the path parametrically in the form ( s). s a path in space Often, for
More informationRotational motion of a rigid body spinning around a rotational axis ˆn;
Physics 106a, Caltech 15 November, 2018 Lecture 14: Rotations The motion of solid bodies So far, we have been studying the motion of point particles, which are essentially just translational. Bodies with
More information. D CR Nomenclature D 1
. D CR Nomenclature D 1 Appendix D: CR NOMENCLATURE D 2 The notation used by different investigators working in CR formulations has not coalesced, since the topic is in flux. This Appendix identifies the
More informationChapter 1. Preliminaries
Chapter 1 Preliminaries 1.1 The Vector Concept Revisited The concept of a vector has been one of the most fruitful ideas in all of mathematics, and it is not surprising that we receive repeated exposure
More informationCE-570 Advanced Structural Mechanics - Arun Prakash
Ch1-Intro Page 1 CE-570 Advanced Structural Mechanics - Arun Prakash The BIG Picture What is Mechanics? Mechanics is study of how things work: how anything works, how the world works! People ask: "Do you
More informationin this web service Cambridge University Press
CONTINUUM MECHANICS This is a modern textbook for courses in continuum mechanics. It provides both the theoretical framework and the numerical methods required to model the behavior of continuous materials.
More informationCovariant Formulation of Electrodynamics
Chapter 7. Covariant Formulation of Electrodynamics Notes: Most of the material presented in this chapter is taken from Jackson, Chap. 11, and Rybicki and Lightman, Chap. 4. Starting with this chapter,
More informationDirectional Field. Xiao-Ming Fu
Directional Field Xiao-Ming Fu Outlines Introduction Discretization Representation Objectives and Constraints Outlines Introduction Discretization Representation Objectives and Constraints Definition Spatially-varying
More informationIntroduction to Vector Spaces
1 CSUC Department of Physics Mechanics: Class Notes Introduction to Vector Spaces I. INTRODUCTION Modern mathematics often constructs logical systems by merely proposing a set of elements that obey a specific
More informationChapter 11. Special Relativity
Chapter 11 Special Relativity Note: Please also consult the fifth) problem list associated with this chapter In this chapter, Latin indices are used for space coordinates only eg, i = 1,2,3, etc), while
More informationClifford Algebras and Spin Groups
Clifford Algebras and Spin Groups Math G4344, Spring 2012 We ll now turn from the general theory to examine a specific class class of groups: the orthogonal groups. Recall that O(n, R) is the group of
More informationPrincipal Component Analysis
Principal Component Analysis Laurenz Wiskott Institute for Theoretical Biology Humboldt-University Berlin Invalidenstraße 43 D-10115 Berlin, Germany 11 March 2004 1 Intuition Problem Statement Experimental
More informationEngineering Sciences 241 Advanced Elasticity, Spring Distributed Thursday 8 February.
Engineering Sciences 241 Advanced Elasticity, Spring 2001 J. R. Rice Homework Problems / Class Notes Mechanics of finite deformation (list of references at end) Distributed Thursday 8 February. Problems
More informationIntroduction to Group Theory
Chapter 10 Introduction to Group Theory Since symmetries described by groups play such an important role in modern physics, we will take a little time to introduce the basic structure (as seen by a physicist)
More informationIntroduction to tensors and dyadics
1 Introduction to tensors and dyadics 1.1 Introduction Tensors play a fundamental role in theoretical physics. The reason for this is that physical laws written in tensor form are independent of the coordinate
More informationCourse No: (1 st version: for graduate students) Course Name: Continuum Mechanics Offered by: Chyanbin Hwu
Course No: (1 st version: for graduate students) Course Name: Continuum Mechanics Offered by: Chyanbin Hwu 2011. 11. 25 Contents: 1. Introduction 1.1 Basic Concepts of Continuum Mechanics 1.2 The Need
More information1 DIFFERENTIATION OF VECTORS
1 DIFFERENTIATION OF VECTORS The discipline of dynamics deals with changes of various kinds, such as changes in the position of a particle in a reference frame and changes in the configuration of a mechanical
More information! EN! EU! NE! EE.! ij! NN! NU! UE! UN! UU
A-1 Appendix A. Equations for Translating Between Stress Matrices, Fault Parameters, and P-T Axes Coordinate Systems and Rotations We use the same right-handed coordinate system as Andy Michael s program,
More informationconventions and notation
Ph95a lecture notes, //0 The Bloch Equations A quick review of spin- conventions and notation The quantum state of a spin- particle is represented by a vector in a two-dimensional complex Hilbert space
More informationA short review of continuum mechanics
A short review of continuum mechanics Professor Anette M. Karlsson, Department of Mechanical ngineering, UD September, 006 This is a short and arbitrary review of continuum mechanics. Most of this material
More informationPART ONE DYNAMICS OF A SINGLE PARTICLE
PART ONE DYNAMICS OF A SINGLE PARTICLE 1 Kinematics of a Particle 1.1 Introduction One of the main goals of this book is to enable the reader to take a physical system, model it by using particles or rigid
More information1. Background. is usually significantly lower than it is in uniaxial tension
NOTES ON QUANTIFYING MODES OF A SECOND- ORDER TENSOR. The mechanical behavior of rocks and rock-like materials (concrete, ceramics, etc.) strongly depends on the loading mode, defined by the values and
More informationFoundations of Matrix Analysis
1 Foundations of Matrix Analysis In this chapter we recall the basic elements of linear algebra which will be employed in the remainder of the text For most of the proofs as well as for the details, the
More informationDiscriminative Direction for Kernel Classifiers
Discriminative Direction for Kernel Classifiers Polina Golland Artificial Intelligence Lab Massachusetts Institute of Technology Cambridge, MA 02139 polina@ai.mit.edu Abstract In many scientific and engineering
More informationINDEX 363. Cartesian coordinates 19,20,42, 67, 83 Cartesian tensors 84, 87, 226
INDEX 363 A Absolute differentiation 120 Absolute scalar field 43 Absolute tensor 45,46,47,48 Acceleration 121, 190, 192 Action integral 198 Addition of systems 6, 51 Addition of tensors 6, 51 Adherence
More informationLinear Algebra (Review) Volker Tresp 2018
Linear Algebra (Review) Volker Tresp 2018 1 Vectors k, M, N are scalars A one-dimensional array c is a column vector. Thus in two dimensions, ( ) c1 c = c 2 c i is the i-th component of c c T = (c 1, c
More informationKinematics. Chapter Multi-Body Systems
Chapter 2 Kinematics This chapter first introduces multi-body systems in conceptual terms. It then describes the concept of a Euclidean frame in the material world, following the concept of a Euclidean
More informationMathematics that Every Physicist should Know: Scalar, Vector, and Tensor Fields in the Space of Real n- Dimensional Independent Variable with Metric
Mathematics that Every Physicist should Know: Scalar, Vector, and Tensor Fields in the Space of Real n- Dimensional Independent Variable with Metric By Y. N. Keilman AltSci@basicisp.net Every physicist
More informationNotes on Mathematics Groups
EPGY Singapore Quantum Mechanics: 2007 Notes on Mathematics Groups A group, G, is defined is a set of elements G and a binary operation on G; one of the elements of G has particularly special properties
More informationAlgebra I Fall 2007
MIT OpenCourseWare http://ocw.mit.edu 18.701 Algebra I Fall 007 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms. 18.701 007 Geometry of the Special Unitary
More informationLecture Notes - Modeling of Mechanical Systems
Thomas Bak Lecture Notes - Modeling of Mechanical Systems February 19, 2002 Aalborg University Department of Control Engineering Fredrik Bajers Vej 7C DK-9220 Aalborg Denmark 2 Table of Contents Table
More informationMathematical Preliminaries
Mathematical Preliminaries Introductory Course on Multiphysics Modelling TOMAZ G. ZIELIŃKI bluebox.ippt.pan.pl/ tzielins/ Table of Contents Vectors, tensors, and index notation. Generalization of the concept
More informationElementary Linear Algebra
Matrices J MUSCAT Elementary Linear Algebra Matrices Definition Dr J Muscat 2002 A matrix is a rectangular array of numbers, arranged in rows and columns a a 2 a 3 a n a 2 a 22 a 23 a 2n A = a m a mn We
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 informationTensor Analysis in Euclidean Space
Tensor Analysis in Euclidean Space James Emery Edited: 8/5/2016 Contents 1 Classical Tensor Notation 2 2 Multilinear Functionals 4 3 Operations With Tensors 5 4 The Directional Derivative 5 5 Curvilinear
More informationPhysics 110. Electricity and Magnetism. Professor Dine. Spring, Handout: Vectors and Tensors: Everything You Need to Know
Physics 110. Electricity and Magnetism. Professor Dine Spring, 2008. Handout: Vectors and Tensors: Everything You Need to Know What makes E&M hard, more than anything else, is the problem that the electric
More informationStatistics 992 Continuous-time Markov Chains Spring 2004
Summary Continuous-time finite-state-space Markov chains are stochastic processes that are widely used to model the process of nucleotide substitution. This chapter aims to present much of the mathematics
More informationHONORS LINEAR ALGEBRA (MATH V 2020) SPRING 2013
HONORS LINEAR ALGEBRA (MATH V 2020) SPRING 2013 PROFESSOR HENRY C. PINKHAM 1. Prerequisites The only prerequisite is Calculus III (Math 1201) or the equivalent: the first semester of multivariable calculus.
More informationWe wish to solve a system of N simultaneous linear algebraic equations for the N unknowns x 1, x 2,...,x N, that are expressed in the general form
Linear algebra This chapter discusses the solution of sets of linear algebraic equations and defines basic vector/matrix operations The focus is upon elimination methods such as Gaussian elimination, and
More informationCourse Syllabus: Continuum Mechanics - ME 212A
Course Syllabus: Continuum Mechanics - ME 212A Division Course Number Course Title Academic Semester Physical Science and Engineering Division ME 212A Continuum Mechanics Fall Academic Year 2017/2018 Semester
More informationTexture, Microstructure & Anisotropy A.D. (Tony) Rollett
1 Carnegie Mellon MRSEC 27-750 Texture, Microstructure & Anisotropy A.D. (Tony) Rollett Last revised: 5 th Sep. 2011 2 Show how to convert from a description of a crystal orientation based on Miller indices
More information02 - introduction to vectors and tensors. me338 - syllabus. introduction tensor calculus. continuum mechanics. introduction. continuum mechanics
02 - introduction to vectors and tensors me338 - syllabus holzapfel nonlinear solid mechanics [2000], chapter 1.1-1.5, pages 1-32 02-1 introduction 2 the potato equations is the branch of mechanics concerned
More informationChapter 8. Rigid transformations
Chapter 8. Rigid transformations We are about to start drawing figures in 3D. There are no built-in routines for this purpose in PostScript, and we shall have to start more or less from scratch in extending
More informationVectors, metric and the connection
Vectors, metric and the connection 1 Contravariant and covariant vectors 1.1 Contravariant vectors Imagine a particle moving along some path in the 2-dimensional flat x y plane. Let its trajectory be given
More informationLecture 10: A (Brief) Introduction to Group Theory (See Chapter 3.13 in Boas, 3rd Edition)
Lecture 0: A (Brief) Introduction to Group heory (See Chapter 3.3 in Boas, 3rd Edition) Having gained some new experience with matrices, which provide us with representations of groups, and because symmetries
More informationConstitutive Equations
Constitutive quations David Roylance Department of Materials Science and ngineering Massachusetts Institute of Technology Cambridge, MA 0239 October 4, 2000 Introduction The modules on kinematics (Module
More informationLecture 2: Linear Algebra Review
EE 227A: Convex Optimization and Applications January 19 Lecture 2: Linear Algebra Review Lecturer: Mert Pilanci Reading assignment: Appendix C of BV. Sections 2-6 of the web textbook 1 2.1 Vectors 2.1.1
More informationIncompatibility Paradoxes
Chapter 22 Incompatibility Paradoxes 22.1 Simultaneous Values There is never any difficulty in supposing that a classical mechanical system possesses, at a particular instant of time, precise values of
More informationIntroduction to Finite Element Method
Introduction to Finite Element Method Dr. Rakesh K Kapania Aerospace and Ocean Engineering Department Virginia Polytechnic Institute and State University, Blacksburg, VA AOE 524, Vehicle Structures Summer,
More information,, rectilinear,, spherical,, cylindrical. (6.1)
Lecture 6 Review of Vectors Physics in more than one dimension (See Chapter 3 in Boas, but we try to take a more general approach and in a slightly different order) Recall that in the previous two lectures
More informationGAME PHYSICS SECOND EDITION. дяййтаййг 1 *
GAME PHYSICS SECOND EDITION DAVID H. EBERLY дяййтаййг 1 * К AMSTERDAM BOSTON HEIDELBERG LONDON NEW YORK OXFORD PARIS SAN DIEGO SAN FRANCISCO SINGAPORE SYDNEY TOKYO MORGAN ELSEVIER Morgan Kaufmann Publishers
More informationCS 143 Linear Algebra Review
CS 143 Linear Algebra Review Stefan Roth September 29, 2003 Introductory Remarks This review does not aim at mathematical rigor very much, but instead at ease of understanding and conciseness. Please see
More informationClassical differential geometry of two-dimensional surfaces
Classical differential geometry of two-dimensional surfaces 1 Basic definitions This section gives an overview of the basic notions of differential geometry for twodimensional surfaces. It follows mainly
More informationMathematical Background
CHAPTER ONE Mathematical Background This book assumes a background in the fundamentals of solid mechanics and the mechanical behavior of materials, including elasticity, plasticity, and friction. A previous
More informationSymmetry and Properties of Crystals (MSE638) Stress and Strain Tensor
Symmetry and Properties of Crystals (MSE638) Stress and Strain Tensor Somnath Bhowmick Materials Science and Engineering, IIT Kanpur April 6, 2018 Tensile test and Hooke s Law Upto certain strain (0.75),
More informationLinear Algebra. Preliminary Lecture Notes
Linear Algebra Preliminary Lecture Notes Adolfo J. Rumbos c Draft date May 9, 29 2 Contents 1 Motivation for the course 5 2 Euclidean n dimensional Space 7 2.1 Definition of n Dimensional Euclidean Space...........
More informationVectors. September 2, 2015
Vectors September 2, 2015 Our basic notion of a vector is as a displacement, directed from one point of Euclidean space to another, and therefore having direction and magnitude. We will write vectors in
More informationQUATERNIONS AND ROTATIONS
QUATERNIONS AND ROTATIONS SVANTE JANSON 1. Introduction The purpose of this note is to show some well-known relations between quaternions and the Lie groups SO(3) and SO(4) (rotations in R 3 and R 4 )
More information1 Euclidean geometry. 1.1 The metric on R n
1 Euclidean geometry This chapter discusses the geometry of n-dimensional Euclidean space E n, together with its distance function. The distance gives rise to other notions such as angles and congruent
More informationElementary maths for GMT
Elementary maths for GMT Linear Algebra Part 1: Vectors, Representations Algebra and Linear Algebra Algebra: numbers and operations on numbers 2 + 3 = 5 3 7 = 21 Linear Algebra: tuples, triples... of numbers
More informationA Tutorial on Data Reduction. Principal Component Analysis Theoretical Discussion. By Shireen Elhabian and Aly Farag
A Tutorial on Data Reduction Principal Component Analysis Theoretical Discussion By Shireen Elhabian and Aly Farag University of Louisville, CVIP Lab November 2008 PCA PCA is A backbone of modern data
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 informationTensor Calculus. arxiv: v1 [math.ho] 14 Oct Taha Sochi. October 17, 2016
Tensor Calculus arxiv:1610.04347v1 [math.ho] 14 Oct 2016 Taha Sochi October 17, 2016 Department of Physics & Astronomy, University College London, Gower Street, London, WC1E 6BT. Email: t.sochi@ucl.ac.uk.
More information14 Singular Value Decomposition
14 Singular Value Decomposition For any high-dimensional data analysis, one s first thought should often be: can I use an SVD? The singular value decomposition is an invaluable analysis tool for dealing
More informationPHYS 705: Classical Mechanics. Rigid Body Motion Introduction + Math Review
1 PHYS 705: Classical Mechanics Rigid Body Motion Introduction + Math Review 2 How to describe a rigid body? Rigid Body - a system of point particles fixed in space i r ij j subject to a holonomic constraint:
More informationContravariant and Covariant as Transforms
Contravariant and Covariant as Transforms There is a lot more behind the concepts of contravariant and covariant tensors (of any rank) than the fact that their basis vectors are mutually orthogonal to
More informationMath 302 Outcome Statements Winter 2013
Math 302 Outcome Statements Winter 2013 1 Rectangular Space Coordinates; Vectors in the Three-Dimensional Space (a) Cartesian coordinates of a point (b) sphere (c) symmetry about a point, a line, and a
More informationELASTOPLASTICITY THEORY by V. A. Lubarda
ELASTOPLASTICITY THEORY by V. A. Lubarda Contents Preface xiii Part 1. ELEMENTS OF CONTINUUM MECHANICS 1 Chapter 1. TENSOR PRELIMINARIES 3 1.1. Vectors 3 1.2. Second-Order Tensors 4 1.3. Eigenvalues and
More informationSome Notes on Linear Algebra
Some Notes on Linear Algebra prepared for a first course in differential equations Thomas L Scofield Department of Mathematics and Statistics Calvin College 1998 1 The purpose of these notes is to present
More informationLecture 1 - Vectors. A Puzzle... Introduction. Vectors: The quest begins! TA s Information. Vectors
Lecture 1 - Vectors Puzzle... The sum starts at 0. Players alternate by choosing an integer from 1 to 10 and then adding it to the sum. The player who gets to 100 wins. You go first. What is the winning
More informationElements of Continuum Elasticity. David M. Parks Mechanics and Materials II February 25, 2004
Elements of Continuum Elasticity David M. Parks Mechanics and Materials II 2.002 February 25, 2004 Solid Mechanics in 3 Dimensions: stress/equilibrium, strain/displacement, and intro to linear elastic
More informationNOTES ON LINEAR ALGEBRA CLASS HANDOUT
NOTES ON LINEAR ALGEBRA CLASS HANDOUT ANTHONY S. MAIDA CONTENTS 1. Introduction 2 2. Basis Vectors 2 3. Linear Transformations 2 3.1. Example: Rotation Transformation 3 4. Matrix Multiplication and Function
More informationinertia of a body, principal axes of inertia, invariants of an inertia tensor, and inertia triangle inequalities are illustrated and discussed.
Preface This book belongs to a series of three books written simultaneously (the remaining two are titled Classical Mechanics: Dynamics and Classical Mechanics: Applied MechanicsandMechatronics). This
More informationLinear Algebra. Preliminary Lecture Notes
Linear Algebra Preliminary Lecture Notes Adolfo J. Rumbos c Draft date April 29, 23 2 Contents Motivation for the course 5 2 Euclidean n dimensional Space 7 2. Definition of n Dimensional Euclidean Space...........
More informationCSE 554 Lecture 7: Alignment
CSE 554 Lecture 7: Alignment Fall 2012 CSE554 Alignment Slide 1 Review Fairing (smoothing) Relocating vertices to achieve a smoother appearance Method: centroid averaging Simplification Reducing vertex
More informationChapter 2. Linear Algebra. rather simple and learning them will eventually allow us to explain the strange results of
Chapter 2 Linear Algebra In this chapter, we study the formal structure that provides the background for quantum mechanics. The basic ideas of the mathematical machinery, linear algebra, are rather simple
More information02 - tensor calculus - tensor algebra tensor calculus
02 - - tensor algebra 02-1 tensor the word tensor was introduced in 1846 by william rowan hamilton. it was used in its current meaning by woldemar voigt in 1899. was developed around 1890 by gregorio ricci-curbastro
More informationLecture 7. Quaternions
Matthew T. Mason Mechanics of Manipulation Spring 2012 Today s outline Motivation Motivation have nice geometrical interpretation. have advantages in representing rotation. are cool. Even if you don t
More informationGeometry for Physicists
Hung Nguyen-Schafer Jan-Philip Schmidt Tensor Analysis and Elementary Differential Geometry for Physicists and Engineers 4 i Springer Contents 1 General Basis and Bra-Ket Notation 1 1.1 Introduction to
More informationMSMS Vectors and Matrices
MSMS Vectors and Matrices Basilio Bona DAUIN Politecnico di Torino Semester 1, 2015-2016 B. Bona (DAUIN) MSMS-Vectors and matrices Semester 1, 2015-2016 1 / 39 Introduction Most of the topics introduced
More informationReadings for Unit I from Ryan (Topics from linear algebra)
Readings for Unit I from Ryan (Topics from linear algebra) I.0 : Background Suggested readings. Ryan : pp. 193 202 General convention. In most cases, the background readings for a section of the course
More informationPreface. Figures Figures appearing in the text were prepared using MATLAB R. For product information, please contact:
Linear algebra forms the basis for much of modern mathematics theoretical, applied, and computational. The purpose of this book is to provide a broad and solid foundation for the study of advanced mathematics.
More informationCOMPOSITE KNOT DETERMINANTS
COMPOSITE KNOT DETERMINANTS SAMANTHA DIXON Abstract. In this paper, we will introduce the basics of knot theory, with special focus on tricolorability, Fox r-colorings of knots, and knot determinants.
More informationReview of Linear Algebra
Review of Linear Algebra Definitions An m n (read "m by n") matrix, is a rectangular array of entries, where m is the number of rows and n the number of columns. 2 Definitions (Con t) A is square if m=
More informationElastic Wave Theory. LeRoy Dorman Room 436 Ritter Hall Tel: Based on notes by C. R. Bentley. Version 1.
Elastic Wave Theory LeRoy Dorman Room 436 Ritter Hall Tel: 4-2406 email: ldorman@ucsd.edu Based on notes by C. R. Bentley. Version 1.1, 20090323 1 Chapter 1 Tensors Elastic wave theory is concerned with
More informationStatistical Geometry Processing Winter Semester 2011/2012
Statistical Geometry Processing Winter Semester 2011/2012 Linear Algebra, Function Spaces & Inverse Problems Vector and Function Spaces 3 Vectors vectors are arrows in space classically: 2 or 3 dim. Euclidian
More informationAPPENDIX A. Background Mathematics. A.1 Linear Algebra. Vector algebra. Let x denote the n-dimensional column vector with components x 1 x 2.
APPENDIX A Background Mathematics A. Linear Algebra A.. Vector algebra Let x denote the n-dimensional column vector with components 0 x x 2 B C @. A x n Definition 6 (scalar product). The scalar product
More informationCaveats concerning conjugate stress and strain measures for frame indifferent anisotropic elasticity
This is a modified version of the paper by the same title that appears in Acta Mech, Vol. 129, No. 1 2, 1998. In this version, the tensor order of a variable equals the number of underlines beneath it.
More informationChapter 3 Stress, Strain, Virtual Power and Conservation Principles
Chapter 3 Stress, Strain, irtual Power and Conservation Principles 1 Introduction Stress and strain are key concepts in the analytical characterization of the mechanical state of a solid body. While stress
More informationLinear Algebra and Robot Modeling
Linear Algebra and Robot Modeling Nathan Ratliff Abstract Linear algebra is fundamental to robot modeling, control, and optimization. This document reviews some of the basic kinematic equations and uses
More informationSTEP Support Programme. STEP 2 Matrices Topic Notes
STEP Support Programme STEP 2 Matrices Topic Notes Definitions............................................. 2 Manipulating Matrices...................................... 3 Transformations.........................................
More information4.2. ORTHOGONALITY 161
4.2. ORTHOGONALITY 161 Definition 4.2.9 An affine space (E, E ) is a Euclidean affine space iff its underlying vector space E is a Euclidean vector space. Given any two points a, b E, we define the distance
More informationCHAPTER 1: Functions
CHAPTER 1: Functions 1.1: Functions 1.2: Graphs of Functions 1.3: Basic Graphs and Symmetry 1.4: Transformations 1.5: Piecewise-Defined Functions; Limits and Continuity in Calculus 1.6: Combining Functions
More informationLinear Algebra and Eigenproblems
Appendix A A Linear Algebra and Eigenproblems A working knowledge of linear algebra is key to understanding many of the issues raised in this work. In particular, many of the discussions of the details
More informationA Primer on Three Vectors
Michael Dine Department of Physics University of California, Santa Cruz September 2010 What makes E&M hard, more than anything else, is the problem that the electric and magnetic fields are vectors, and
More information1 Mathematical preliminaries
1 Mathematical preliminaries The mathematical language of quantum mechanics is that of vector spaces and linear algebra. In this preliminary section, we will collect the various definitions and mathematical
More information1 Differentiable manifolds and smooth maps
1 Differentiable manifolds and smooth maps Last updated: April 14, 2011. 1.1 Examples and definitions Roughly, manifolds are sets where one can introduce coordinates. An n-dimensional manifold is a set
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 information1 Tensors and relativity
Physics 705 1 Tensors and relativity 1.1 History Physical laws should not depend on the reference frame used to describe them. This idea dates back to Galileo, who recognized projectile motion as free
More informationMobile Robotics 1. A Compact Course on Linear Algebra. Giorgio Grisetti
Mobile Robotics 1 A Compact Course on Linear Algebra Giorgio Grisetti SA-1 Vectors Arrays of numbers They represent a point in a n dimensional space 2 Vectors: Scalar Product Scalar-Vector Product Changes
More information