Supplementary Technical Document
|
|
- Clara Dickerson
- 5 years ago
- Views:
Transcription
1 Supplementary Technical Document Qi Guo, Xuchen Han, Chuyuan Fu, Theodore Gast, Rasmus Tamstorf, Joseph Teran 1 FEM Force computation We compute forces on the control points x p by fp KL = ΨS F KL,Etr x KL x KL p = V 0 ψf KL,Etr q x q q x KL q p = q V 0 q ψ F FKL,Etr KL q x q : FKL,Etr q x KL x q, p where x q s are positions of the quadrature points. We give expressions for each FKL q x KL p k x q with fixed p, q and k, where k represents the x, y, or z direction. For simplicity of notation, we omit the subscripts p, q and superscript KL for now. Recall from the paper that we have where F = 3 g i ḡ i, with g α = a α + ξ 3 a 3,α, g 3 = a 3, i=1 a α = j Nj SD x j ξ 1, ξ 2, α = 1, 2 ξ α a 3 = a 1 a 2 a 3,α = I a 3 a 3 a 1,α a 2 + a 1 a 2,α = ã a 3 a 3 ã 1
2 in which we define ã to be Now we compute F. ã = a 1,α a 2 + a 1 a 2,α. F = 3 i=1 g i ḡ i, and where g α = a α + ξ 3 a 3,α g 3 = a 3 a α = N SD k ξ 1, ξ 2 ξ α e k summation convention does not apply here 1 and Finally, a 3 = a 1 a 2 + a 1 a 2 a 1 a 2 a 3 = a 3 a 1 a 2 + a 1 a 2, a 3,α = ã a 3 a 3 ã ã + a 3 a 3 a 3 ã, where ã = a 1,α a 2 + a 1,α a 2 + a 1 a 2,α + a 1 a 2,α a 1,α a 2 + a 1 a 2,α 2, in which a α,β = N SD k ξ 1, ξ 2 e k summation convention does not apply here. ξ β ξ α 2
3 2 Grid force computation The force on the MPM grid fi iii x computes as follows: f iii i x = p I iii = p I iii + p I iii χa pα ā pα + a E p3 ā p3 Vp 0 χa pα ā pα + a E p3 ā p3 : apα ā pα + a E p3 ā p3 F E a pβ χa pα ā pα + a E p3 ā p3 F E : apα ā pα + a E p3 ā p3 a E p3 Then, omitting the subscript p, we compute each term in the contraction: χa α ā α + a E 3 ā 3 F E = τ S a α ā α + a E 3 ā 3 T = τ S ã α ā α + ã 3 ā 3 where τ S is the Kirchhoff stress and ã α and ã 3 are the contravariant counterparts of a α and a E 3 respectively. And using index notation, we see that a α ā α + a E 3 ā 3 = a α i ā αj a β a βk = δ αβ δ ik ā αj = δ ik ā βj : a pβ Vp 0 x i : ae p3 x i V 0 p. Similarly, a α ā α + a E 3 ā 3 a E = δ ik ā 3j 3 Hence, contracting the first two terms in the summation, each term in the summation becomes τ S ã α ā α + ã 3 ā 3 ā β : a β + τ S ã α ā α + ã 3 ā 3 ā 3 : ae 3 = τ S ã β : a β + τ S ã 3 : ae 3 3
4 Note that a β = a β x p x p = a β x p w n ip, and the expression for a β x p is given equation 1. Ignoring further plastic flow, we have a E 3 x = j x j wn jp a E,n 3, and thus, a E 3 = w n ip ae,n 3 Therefore, we arrive at the final expression for the force of type iii: f iii i x = 3 Laminate Stress p I iii τ S p ãβ p : In this section we derive the expression for τ KL = τ αβ q KL,E α a pβ wip n x + τ S p ã3 p : wip n ae,n p3 p q KL,E β, τ KL αβ = 2µɛL αβ + λɛl γγδ αβ. 2 First notice that we may replace the right Hencky strain with left Hencky strain in the definition of energy because of the isotropic nature of the energy function. We now give the drivation of Equation 2 with index free notation assuming all variables are in 2D. ψf = ψuσv T because the energy is isotropic. Hence, PF = PUΣV T = UPΣV T PF = UPΣV T = U ψ Σ VT = U 2µ logσσ 1 + λtrlog ΣΣ 1 V T. 4
5 Therefore, τ KL = U 2µ logσσ 1 + λtrlog ΣΣ 1 V T F T = U 2µ logσ + λtrlog Σ U T = 2µɛ L + λtrɛ L 4 QR and Elastic Potential We can use QR orthogonalization of deformed material directions to define q i r ij = Fā j, F = r ij q i ā j, r ij = 0 for i > j Change of basis tensor Define the change of basis tensor Q = Q ij ā i ā j 4 with Q ij = q j ā i. With this convention we see that Qā i = q i and Q T Q = I. Furthermore, defining R = r ij ā i ā j we have F = QR. 4.2 Differentials The QR differential satisfies q k δq i r ij + δr kj = q k δfā j, δf = δr ij q i ā j + r ij δq i ā j 5 where q k δq i = q i δq k from orthogonality of the q i. And δf = δqr + QδR 6 where δq T Q = Q T δq from Q T Q = I. Furthermore, and the δr ij = 0 for i > j. δq = δq ij ā i ā j, δq ij = δq j ā i, δq i = δqā i 7 δr = δr ij ā i ā j 8 5
6 5 Elastic potential and stresses Define the hyperelastic potential as ψf = ˆψ[R] 9 where [R] = r 11 r 12 r 13 r 22 r The differential satisfies δψf = ψ F : δf = P : δf = [R]δr ij 11 F where P = ψ F F. Therefore Similarly, δr ij q i Pā j + r ij δq i Pā j = [R]δr ij. 12 P : δf = P : δqr + P : QδR = [R]δr ij 13 Choosing δf = δr ij q i ā j i.e. δq i = 0, we can conclude that q i Pā j δr ij = [R]δr ij 14 for arbitrary δr ij with i j. Therefore the q i Pā j = [R] for i j. Similarly, P : QδR = Q T P : δr = δr ij ā i Q T Pā j = δrij q i Pā j = [R]δr ij. Choosing δf = r ij δq i ā j i.e. δr ij = 0, we can conclude that 15 0 = r ij δq i Pā j. 16 6
7 Similarly, 0 = P : δqr = PR T : δq = PR T : δqq T Q = PR T Q T : δqq T = PF T : δqq T 17 In other words, the Kirchhoff stress τ = PF T is symmetric since δqq T is arbitrary skew. Furthermore, and we know P ij = τ 11 τ 12 τ 13 τ 21 τ 22 τ 23 = τ 31 τ 32 τ 33 P = P ij q i ā j, τ = P ij r kj q i q k = τ ik q i q k 18 = for i j from Equation 14. Thus P 11 P 12 P 13 P 21 P 22 P 23 P 31 P 32 P 33 r 11 r 12 r 22 r 13 r P 11 r 11 + P 12 r 12 + P 13 r 13 P 12 r 22 + P 13 r 32 P 13 P 21 r 11 + P 22 r 12 + P 23 r 13 P 22 r 22 + P 23 r 32 P 23 P 31 r 11 + P 32 r 12 + P 33 r 13 P 32 r 22 + P 33 r 32 P 33 20, and since τ = τ T and P ij = for i j, τ 11 τ 12 τ 13 τ 21 τ 22 τ 23 τ 31 τ 32 τ 33 = r 11 r 11 + r 12 r 12 + r 12 r 22 + r 13 r 13 r 32 r 13 r 13 r 12 r 22 + r 13 r 32 r 22 r 22 + r 23 r 23 r 32 r 13 r In particular, the matrix representation of τ S reads τ 11 τ 12 τ 13 τ 21 τ 22 τ 23 τ 31 τ 32 τ 33 = = 0 0 γs γs f s 3 γs 2 1 γs 1 s 2 γs 1 s 3 γs 1 s 2 γs 2 2 γs 2 s 3 γs 1 s 3 γs 2 s 3 f s s 1 s 2 s
8 6 Frictional Contact Yield Condition Coulomb friction places a constraint on the stress as t S c F σ n 24 where σ n = a KL 3 σa KL 3. Recall that a KL 3 = q 3 and thus σ n = q 3 σq 3. On the other hand, t S is the tangential component of the force density and has the form t S = cq 1 + sq 2 σq 3 for some c and s such that c 2 + s 2 = 1. Hence, we may rewrite the constraint on stress as cq 1 + sq 2 σq 3 + c F q 3 σq Using the fact that σ = detfτ, we rewrite the constraint as cq 1 + sq 2 τ q 3 + c F q 3 τ q Substituting in the expression for τ from equation 23, we find that the maximum on the left-hand-side is ±γs 3 s s2 2 + c F f s 3 We apply the particular form of f in the paper where fx = 1 3 kc 1 x 3 for x 1 and 0 otherwise. When s 3 > 1, the maximum is γs 3 s s 2 2. In this case the return mapping set s 1 and s 2 to 0. If 0 < s 3 1, the maximum is and thus we need γs 3 s s2 2 c F k c s s 3, s s2 2 c F k c γ 1 s 3 2. In this case we uniformly scale back s 1 and s 2 to satisfy the constraint. 7 Denting Yield Condition and Return Mapping We apply the von Mises yield condition to the Kirchhoff-Stress in Equation 2 This condition states that the deviatoric component of the stress is less than a threshold value c vm f vm τ = τ trτ 3 I F c vm. 27 8
9 This condition defines a cylindrical region of feasible states in the principal stress space since 2 f vm τ = 3 τ τ 22 + τ 32 τ 1 τ 2 + τ 2 τ 3 + τ 1 τ 3 28 where τ = i τ iu i u i with principal stresses τ i. The plane stress nature of τ KL = α τ α KL u α u α means that feasible stresses are those where the principal stresses are in the ellipsoidal intersection of the cylinder and the τα KL plane. The yield condition is satisfied via associative projection or return mapping of the stress to the feasible region. The elastic and plastic strains are then computed to be consistent with the projected stress. We use F KL,Etr, F KL,P tr to denote the trial state of elastoplastic strains with associated trial stress τ KLtr. We use F KL,E, F KL,P, τ KL to denote their projected counterparts. F KL,Etr, F KL,P tr, τ KLtr F KL,E, F KL,P, τ KL. 29 The deformation gradient constraint must be equal to the product of trial and projected elastic and plastic deformation gradients, creating the constraint on the projection F KL = F KL,Etr F KL,P tr = F KL,E F KL,P. 30 The projection is completed by first computing the trial state of stress τ KLtr from F KL,Etr using Equation 2. This is done by computing the QR decomposition of the trial elastic deformation gradient F KL,Etr = r KL,E tr αβ q KL,E α ā β + q KL,E 3 ā 3. Then we compute the SVD of matrix [r KL,Etr ] R 2 2 and the trial strain [ɛ Ltr ] [r KL,Etr ] = [U E ] [ɛ Ltr ] = [U E ] σ1 E tr logσ E 1 tr σ2 E tr [V E ] T 31 logσ2 E tr [U E ] T 32 From Equation 2 we see that the two non-zero principal stresses τ KLtr α of τ KLtr are equal to the eigenvalues of the matrix [τ KLtr ] [τ KLtr ] = 2µ[ɛ Ltr ] + λtr[ɛ Ltr ]I = [U E ] τ1 KL tr τ2 KL tr [U E ] T. 33 9
10 We therefore project the eigenvalues τ KLtr α τ KL α into the ellipsoidal intersection the von Mises yield surface and the τ 1, τ 2 plane in the direction that maximizes energy dissipation. We approximate this region by the diamond shaped region whose boundaries have slopes of ±1 to simplify the return mapping. Note that the direction of the return that maximizes energy dissipation is a function of the Cauchy-Green strain derivative of the Kirchhoff stress and thus is non-trivial to find in general. Fortunately, the quadratic Hencky strain model has the favorable property that the return direction is perpendicular to the yield surface [1] which greatly simplifies the return mapping. We illustrate this property in Figure 1. After projection, we rebuild the matrix without changing the eigenvectors and rebuild τ KL from the matrix [τ KL ] = [U E τ KL ] 1 τ KL 2 [U E ] T, τ KL = ταβ KL qkl,e α q KL,E β 34 where ταβ KL are the entries in the projected matrix [τ KL ] R 2 2. The projected strain [ɛ L ] is computed from the projected principal stresses from [ɛ L ] = [U E logσ E ] 1 logσ2 E [U E ] T 35 logσ E 1 1 2µ + λ λ τ KL logσ2 E = 1 36 λ 2µ + λ τ KL 2 and the projected elastic deformation gradient is F KL,E = F KL,E αβ q KL,E α ā β + q KL,E 3 ā 3 where [ˆF KL,E ] = [U E σ E ] 1 [V E ] T. 37 The projected plastic deformation gradient is computed from F KL,P = F KL,E 1 F KL in order to maintain the constraint in Equation 30. References [1] C. Mast. Modeling landslide-induced flow interactions with structures using the Material Point Method. PhD thesis, σ E 2 10
11 a b Figure 1: Return Mapping. In general in the return mapping direction is non trivial left. Quadratic Hencky strain energy density simplifies the return mapping right. 11
Anisotropic elastoplasticity for cloth, knit and hair frictional contact supplementary technical document
Anisotroic elastolasticity for cloth, knit and hair frictional contact sulementary technical document Chenfanfu Jiang, Theodore Gast, Joseh Teran 1 QR differentiation Here we discuss how to comute δq and
More informationConcept Question Comment on the general features of the stress-strain response under this loading condition for both types of materials
Module 5 Material failure Learning Objectives review the basic characteristics of the uni-axial stress-strain curves of ductile and brittle materials understand the need to develop failure criteria for
More informationConstitutive models: Incremental plasticity Drücker s postulate
Constitutive models: Incremental plasticity Drücker s postulate if consistency condition associated plastic law, associated plasticity - plastic flow law associated with the limit (loading) surface Prager
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 informationMODELING OF CONCRETE MATERIALS AND STRUCTURES. Kaspar Willam. Uniaxial Model: Strain-Driven Format of Elastoplasticity
MODELING OF CONCRETE MATERIALS AND STRUCTURES Kaspar Willam University of Colorado at Boulder Class Meeting #3: Elastoplastic Concrete Models Uniaxial Model: Strain-Driven Format of Elastoplasticity Triaxial
More informationA Brief Introduction to Tensors
A Brief Introduction to Tensors Jay R Walton Fall 2013 1 Preliminaries In general, a tensor is a multilinear transformation defined over an underlying finite dimensional vector space In this brief introduction,
More informationConstitutive models. Constitutive model: determines P in terms of deformation
Constitutive models Constitutive model: determines P in terms of deformation Elastic material: P depends only on current F Hyperelastic material: work is independent of path strain energy density function
More informationTheory of Plasticity. Lecture Notes
Theory of Plasticity Lecture Notes Spring 2012 Contents I Theory of Plasticity 1 1 Mechanical Theory of Plasticity 2 1.1 Field Equations for A Mechanical Theory.................... 2 1.1.1 Strain-displacement
More informationThe Finite Element Method for the Analysis of Non-Linear and Dynamic Systems. Prof. Dr. Eleni Chatzi Lecture ST1-19 November, 2015
The Finite Element Method for the Analysis of Non-Linear and Dynamic Systems Prof. Dr. Eleni Chatzi Lecture ST1-19 November, 2015 Institute of Structural Engineering Method of Finite Elements II 1 Constitutive
More informationContinuum mechanics V. Constitutive equations. 1. Constitutive equation: definition and basic axioms
Continuum mechanics office Math 0.107 ales.janka@unifr.ch http://perso.unifr.ch/ales.janka/mechanics Mars 16, 2011, Université de Fribourg 1. Constitutive equation: definition and basic axioms Constitutive
More information2.073 Notes on finite-deformation theory of isotropic elastic-viscoplastic solids
2.073 Notes on finite-deformation theory of isotropic elastic-viscoplastic solids Lallit Anand Department of Mechanical Engineering Massachusetts Institute of Technology Cambridge, MA 02139, USA April
More informationYou may not start to read the questions printed on the subsequent pages until instructed to do so by the Invigilator.
MATHEMATICAL TRIPOS Part III Thursday 1 June 2006 1.30 to 4.30 PAPER 76 NONLINEAR CONTINUUM MECHANICS Attempt FOUR questions. There are SIX questions in total. The questions carry equal weight. STATIONERY
More information(2) 2. (3) 2 Using equation (3), the material time derivative of the Green-Lagrange strain tensor can be obtained as: 1 = + + +
LAGRANGIAN FORMULAION OF CONINUA Review of Continuum Kinematics he reader is referred to Belytscho et al. () for a concise review of the continuum mechanics concepts used here. he notation followed here
More informationL8. Basic concepts of stress and equilibrium
L8. Basic concepts of stress and equilibrium Duggafrågor 1) Show that the stress (considered as a second order tensor) can be represented in terms of the eigenbases m i n i n i. Make the geometrical representation
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 informationOn the Numerical Modelling of Orthotropic Large Strain Elastoplasticity
63 Advances in 63 On the Numerical Modelling of Orthotropic Large Strain Elastoplasticity I. Karsaj, C. Sansour and J. Soric Summary A constitutive model for orthotropic yield function at large strain
More informationCONTINUUM MECHANICS. lecture notes 2003 jp dr.-ing. habil. ellen kuhl technical university of kaiserslautern
CONTINUUM MECHANICS lecture notes 2003 jp dr.-ing. habil. ellen kuhl technical university of kaiserslautern Contents Tensor calculus. Tensor algebra.................................... Vector algebra.................................
More informationElastic Fields of Dislocations in Anisotropic Media
Elastic Fields of Dislocations in Anisotropic Media a talk given at the group meeting Jie Yin, David M. Barnett and Wei Cai November 13, 2008 1 Why I want to give this talk Show interesting features on
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 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 informationPEAT SEISMOLOGY Lecture 2: Continuum mechanics
PEAT8002 - SEISMOLOGY Lecture 2: Continuum mechanics Nick Rawlinson Research School of Earth Sciences Australian National University Strain Strain is the formal description of the change in shape of a
More informationSEMM Mechanics PhD Preliminary Exam Spring Consider a two-dimensional rigid motion, whose displacement field is given by
SEMM Mechanics PhD Preliminary Exam Spring 2014 1. Consider a two-dimensional rigid motion, whose displacement field is given by u(x) = [cos(β)x 1 + sin(β)x 2 X 1 ]e 1 + [ sin(β)x 1 + cos(β)x 2 X 2 ]e
More informationContinuum Mechanics and the Finite Element Method
Continuum Mechanics and the Finite Element Method 1 Assignment 2 Due on March 2 nd @ midnight 2 Suppose you want to simulate this The familiar mass-spring system l 0 l y i X y i x Spring length before/after
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 information(MPa) compute (a) The traction vector acting on an internal material plane with normal n ( e1 e
EN10: Continuum Mechanics Homework : Kinetics Due 1:00 noon Friday February 4th School of Engineering Brown University 1. For the Cauchy stress tensor with components 100 5 50 0 00 (MPa) compute (a) The
More information06 - concept of stress concept of stress concept of stress concept of stress. me338 - syllabus. definition of stress
holzapfel nonlinear solid mechanics [2000], chapter 3, pages 109-129 holzapfel nonlinear solid mechanics [2000], chapter 3, pages 109-129 1 2 me338 - syllabus definition of stress stress [ stres] is a
More informationBAR ELEMENT WITH VARIATION OF CROSS-SECTION FOR GEOMETRIC NON-LINEAR ANALYSIS
Journal of Computational and Applied Mechanics, Vol.., No. 1., (2005), pp. 83 94 BAR ELEMENT WITH VARIATION OF CROSS-SECTION FOR GEOMETRIC NON-LINEAR ANALYSIS Vladimír Kutiš and Justín Murín Department
More information21 Symmetric and skew-symmetric matrices
21 Symmetric and skew-symmetric matrices 21.1 Decomposition of a square matrix into symmetric and skewsymmetric matrices Let C n n be a square matrix. We can write C = (1/2)(C + C t ) + (1/2)(C C t ) =
More informationProfessor George C. Johnson. ME185 - Introduction to Continuum Mechanics. Midterm Exam II. ) (1) x
Spring, 997 ME85 - Introduction to Continuum Mechanics Midterm Exam II roblem. (+ points) (a) Let ρ be the mass density, v be the velocity vector, be the Cauchy stress tensor, and b be the body force per
More informationTensor Transformations and the Maximum Shear Stress. (Draft 1, 1/28/07)
Tensor Transformations and the Maximum Shear Stress (Draft 1, 1/28/07) Introduction The order of a tensor is the number of subscripts it has. For each subscript it is multiplied by a direction cosine array
More informationIn this section, thermoelasticity is considered. By definition, the constitutive relations for Gradθ. This general case
Section.. Thermoelasticity In this section, thermoelasticity is considered. By definition, the constitutive relations for F, θ, Gradθ. This general case such a material depend only on the set of field
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 informationDurham Research Online
Durham Research Online Deposited in DRO: 15 May 2017 Version of attached le: Accepted Version Peer-review status of attached le: Peer-reviewed Citation for published item: Charlton, T.J. and Coombs, W.M.
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 informationComparative Analysis of Mesh Generators and MIC(0) Preconditioning of FEM Elasticity Systems
Comparative Analysis of Mesh Generators and MIC(0) Preconditioning of FEM Elasticity Systems Nikola Kosturski and Svetozar Margenov Institute for Parallel Processing, Bulgarian Academy of Sciences Abstract.
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 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 informationStress and fabric in granular material
THEORETICAL & APPLIED MECHANICS LETTERS 3, 22 (23) Stress and fabric in granular material Ching S. Chang,, a) and Yang Liu 2 ) Department of Civil Engineering, University of Massachusetts Amherst, Massachusetts
More informationLagrange Multipliers
Optimization with Constraints As long as algebra and geometry have been separated, their progress have been slow and their uses limited; but when these two sciences have been united, they have lent each
More informationVector and tensor calculus
1 Vector and tensor calculus 1.1 Examples Example 1.1 Consider three vectors a = 2 e 1 +5 e 2 b = 3 e1 +4 e 3 c = e 1 given with respect to an orthonormal Cartesian basis { e 1, e 2, e 3 }. a. Compute
More informationHÅLLFASTHETSLÄRA, LTH Examination in computational materials modeling
HÅLLFASTHETSLÄRA, LTH Examination in computational materials modeling TID: 2016-28-10, kl 14.00-19.00 Maximalt 60 poäng kan erhållas på denna tenta. För godkänt krävs 30 poäng. Tillåtet hjälpmedel: räknare
More informationClassification of Prostate Cancer Grades and T-Stages based on Tissue Elasticity Using Medical Image Analysis. Supplementary Document
Classification of Prostate Cancer Grades and T-Stages based on Tissue Elasticity Using Medical Image Analysis Supplementary Document Shan Yang, Vladimir Jojic, Jun Lian, Ronald Chen, Hongtu Zhu, Ming C.
More informationChemnitz Scientific Computing Preprints
Arnd Meyer The Koiter shell equation in a coordinate free description CSC/1-0 Chemnitz Scientific Computing Preprints ISSN 1864-0087 Chemnitz Scientific Computing Preprints Impressum: Chemnitz Scientific
More informationCode_Aster. Constitutive law élasto (visco) plastic in large deformations with metallurgical transformations
Titre : Loi de comportement élasto(visco)plastique en gran[...] Date : 10/08/2010 Page : 1/20 Constitutive law élasto (visco) plastic in large deformations with metallurgical transformations Abstract This
More informationLecture Notes 5
1.5 Lecture Notes 5 Quantities in Different Coordinate Systems How to express quantities in different coordinate systems? x 3 x 3 P Direction Cosines Axis φ 11 φ 3 φ 1 x x x x 3 11 1 13 x 1 3 x 3 31 3
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 informationMechanics of materials Lecture 4 Strain and deformation
Mechanics of materials Lecture 4 Strain and deformation Reijo Kouhia Tampere University of Technology Department of Mechanical Engineering and Industrial Design Wednesday 3 rd February, 206 of a continuum
More informationFunctions of Several Variables
Functions of Several Variables The Unconstrained Minimization Problem where In n dimensions the unconstrained problem is stated as f() x variables. minimize f()x x, is a scalar objective function of vector
More informationContinuum Mechanics Fundamentals
Continuum Mechanics Fundamentals James R. Rice, notes for ES 220, 12 November 2009; corrections 9 December 2009 Conserved Quantities Let a conseved quantity have amount F per unit volume. Examples are
More informationDIAGONALIZATION OF THE STRESS TENSOR
DIAGONALIZATION OF THE STRESS TENSOR INTRODUCTION By the use of Cauchy s theorem we are able to reduce the number of stress components in the stress tensor to only nine values. An additional simplification
More informationBone Tissue Mechanics
Bone Tissue Mechanics João Folgado Paulo R. Fernandes Instituto Superior Técnico, 2016 PART 1 and 2 Introduction The objective of this course is to study basic concepts on hard tissue mechanics. Hard tissue
More informationConstitutive Relations
Constitutive Relations Andri Andriyana, Ph.D. Centre de Mise en Forme des Matériaux, CEMEF UMR CNRS 7635 École des Mines de Paris, 06904 Sophia Antipolis, France Spring, 2008 Outline Outline 1 Review of
More informationLecture #6: 3D Rate-independent Plasticity (cont.) Pressure-dependent plasticity
Lecture #6: 3D Rate-independent Plasticity (cont.) Pressure-dependent plasticity by Borja Erice and Dirk Mohr ETH Zurich, Department of Mechanical and Process Engineering, Chair of Computational Modeling
More informationSimulation of elasticity, biomechanics and virtual surgery. Joseph M. Teran
Contents Simulation of elasticity, biomechanics and virtual surgery Joseph M. Teran 1 Simulation of elasticity, biomechanics and virtual surgery 3 Introduction 3 Real-time computing 4 Continuum mechanics
More informationarxiv: v1 [cond-mat.soft] 31 Dec 2016
A new class of plastic flow evolution equations for anisotropic multiplicative elastoplasticity based on the notion of a corrector elastic strain rate arxiv:1701.00095v1 [cond-mat.soft] 31 Dec 2016 Abstract
More informationPlasticity R. Chandramouli Associate Dean-Research SASTRA University, Thanjavur
Plasticity R. Chandramouli Associate Dean-Research SASTRA University, Thanjavur-613 401 Joint Initiative of IITs and IISc Funded by MHRD Page 1 of 9 Table of Contents 1. Plasticity:... 3 1.1 Plastic Deformation,
More informationIII. TRANSFORMATION RELATIONS
III. TRANSFORMATION RELATIONS The transformation relations from cartesian coordinates to a general curvilinear system are developed here using certain concepts from differential geometry and tensor analysis,
More informationMechanics PhD Preliminary Spring 2017
Mechanics PhD Preliminary Spring 2017 1. (10 points) Consider a body Ω that is assembled by gluing together two separate bodies along a flat interface. The normal vector to the interface is given by n
More informationLaw of behavior élasto (visco) plastic in great deformations with transformations metallurgical
Titre : Loi de comportement élasto(visco)plastique en gran[...] Date : 25/09/2013 Page : 1/20 Law of behavior élasto (visco) plastic in great deformations with transformations metallurgical Summary This
More informationMATH45061: SOLUTION SHEET 1 II
MATH456: SOLUTION SHEET II. The deformation gradient tensor has Cartesian components given by F IJ R I / r J ; and so F R e x, F R, F 3 R, r r r 3 F R r, F R r, F 3 R r 3, F 3 R 3 r, F 3 R 3 r, F 33 R
More informationThe Effect of Evolving Damage on the Finite Strain Response of Inelastic and Viscoelastic Composites
Materials 2009, 2, 858-894; doi:0.3390/ma204858 Article OPEN ACCESS materials ISSN 996-944 www.mdpi.com/journal/materials The Effect of Evolving Damage on the Finite Strain Response of Inelastic and Viscoelastic
More informationAA242B: MECHANICAL VIBRATIONS
AA242B: MECHANICAL VIBRATIONS 1 / 57 AA242B: MECHANICAL VIBRATIONS Dynamics of Continuous Systems These slides are based on the recommended textbook: M. Géradin and D. Rixen, Mechanical Vibrations: Theory
More informationFEM for elastic-plastic problems
FEM for elastic-plastic problems Jerzy Pamin e-mail: JPamin@L5.pk.edu.pl With thanks to: P. Mika, A. Winnicki, A. Wosatko TNO DIANA http://www.tnodiana.com FEAP http://www.ce.berkeley.edu/feap Lecture
More informationSome Aspects of a Discontinuous Galerkin Formulation for Gradient Plasticity at Finite Strains
Some Aspects of a Discontinuous Galerkin Formulation for Gradient Plasticity at Finite Strains Andrew McBride 1 and B. Daya Reddy 1,2 1 Centre for Research in Computational and Applied Mechanics, University
More informationGKSS Research Centre Geesthacht
Preprint 02-2009 GKSS Research Centre Geesthacht Materials Mechanics A variational formulation for finite deformation wrinkling analysis of inelastic membranes J. Mosler & F. Cirak This is a preprint of
More informationRANS Equations in Curvilinear Coordinates
Appendix C RANS Equations in Curvilinear Coordinates To begin with, the Reynolds-averaged Navier-Stokes RANS equations are presented in the familiar vector and Cartesian tensor forms. Each term in the
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 information2 Tensor Notation. 2.1 Cartesian Tensors
2 Tensor Notation It will be convenient in this monograph to use the compact notation often referred to as indicial or index notation. It allows a strong reduction in the number of terms in an equation
More informationCIVL4332 L1 Introduction to Finite Element Method
CIVL L Introduction to Finite Element Method CIVL L Introduction to Finite Element Method by Joe Gattas, Faris Albermani Introduction The FEM is a numerical technique for solving physical problems such
More informationThe Finite Element Method for Computational Structural Mechanics
The Finite Element Method for Computational Structural Mechanics Martin Kronbichler Applied Scientific Computing (Tillämpad beräkningsvetenskap) January 29, 2010 Martin Kronbichler (TDB) FEM for CSM January
More informationA Numerical Study of Finite Element Calculations for Incompressible Materials under Applied Boundary Displacements
A Numerical Study of Finite Element Calculations for Incompressible Materials under Applied Boundary Displacements A Thesis Submitted to the College of Graduate Studies and Research in Partial Fulfillment
More informationIntroduction to Geometry
Introduction to Geometry it is a draft of lecture notes of H.M. Khudaverdian. Manchester, 18 May 211 Contents 1 Euclidean space 3 1.1 Vector space............................ 3 1.2 Basic example of n-dimensional
More informationMODELING OF CONCRETE MATERIALS AND STRUCTURES. Kaspar Willam. Isotropic Elastic Models: Invariant vs Principal Formulations
MODELING OF CONCRETE MATERIALS AND STRUCTURES Kaspar Willam University of Colorado at Boulder Class Meeting #2: Nonlinear Elastic Models Isotropic Elastic Models: Invariant vs Principal Formulations Elastic
More informationLarge strain anisotropic plasticity including effects of plastic spin
377 Large strain anisotropic plasticity including effects of plastic spin Francisco Javier Monta ns a and Klaus-Ju rgen Bathe b, * a Universidad de Castilla-La Mancha, Escuela Te cnica Superior de Ingenieros
More informationarxiv: v2 [cond-mat.mtrl-sci] 23 Oct 2007 ANDREW N. NORRIS
EULERIAN CONJUGATE STRESS AND STRAIN arxiv:0708.2736v2 [cond-mat.mtrl-sci] 23 Oct 2007 ANDREW N. NORRIS Abstract. New results are presented for the stress conjugate to arbitrary Eulerian strain measures.
More information1 Useful Definitions or Concepts
1 Useful Definitions or Concepts 1.1 Elastic constitutive laws One general type of elastic material model is the one called Cauchy elastic material, which depend on only the current local deformation of
More informationMeasurement of deformation. Measurement of elastic force. Constitutive law. Finite element method
Deformable Bodies Deformation x p(x) Given a rest shape x and its deformed configuration p(x), how large is the internal restoring force f(p)? To answer this question, we need a way to measure deformation
More informationChemnitz Scientific Computing Preprints
Arnd Meyer The linear Naghdi shell equation in a coordinate free description CSC/13-03 Chemnitz Scientific Computing Preprints Impressum: Chemnitz Scientific Computing Preprints ISSN 1864-0087 (1995 2005:
More information16.21 Techniques of Structural Analysis and Design Spring 2003 Unit #5 - Constitutive Equations
6.2 Techniques of Structural Analysis and Design Spring 2003 Unit #5 - Constitutive quations Constitutive quations For elastic materials: If the relation is linear: Û σ ij = σ ij (ɛ) = ρ () ɛ ij σ ij =
More information1. Basic Operations Consider two vectors a (1, 4, 6) and b (2, 0, 4), where the components have been expressed in a given orthonormal basis.
Questions on Vectors and Tensors 1. Basic Operations Consider two vectors a (1, 4, 6) and b (2, 0, 4), where the components have been expressed in a given orthonormal basis. Compute 1. a. 2. The angle
More informationThe Finite Element Method II
[ 1 The Finite Element Method II Non-Linear finite element Use of Constitutive Relations Xinghong LIU Phd student 02.11.2007 [ 2 Finite element equilibrium equations: kinematic variables Displacement Strain-displacement
More informationBindel, Fall 2016 Matrix Computations (CS 6210) Notes for At a high level, there are two pieces to solving a least squares problem:
1 Trouble points Notes for 2016-09-28 At a high level, there are two pieces to solving a least squares problem: 1. Project b onto the span of A. 2. Solve a linear system so that Ax equals the projected
More informationAnd similarly in the other directions, so the overall result is expressed compactly as,
SQEP Tutorial Session 5: T7S0 Relates to Knowledge & Skills.5,.8 Last Update: //3 Force on an element of area; Definition of principal stresses and strains; Definition of Tresca and Mises equivalent stresses;
More informationCrystal Micro-Mechanics
Crystal Micro-Mechanics Lectre Classical definition of stress and strain Heng Nam Han Associate Professor School of Materials Science & Engineering College of Engineering Seol National University Seol
More informationSolutions for Fundamentals of Continuum Mechanics. John W. Rudnicki
Solutions for Fundamentals of Continuum Mechanics John W. Rudnicki December, 015 ii Contents I Mathematical Preliminaries 1 1 Vectors 3 Tensors 7 3 Cartesian Coordinates 9 4 Vector (Cross) Product 13 5
More informationYURI LEVIN, MIKHAIL NEDIAK, AND ADI BEN-ISRAEL
Journal of Comput. & Applied Mathematics 139(2001), 197 213 DIRECT APPROACH TO CALCULUS OF VARIATIONS VIA NEWTON-RAPHSON METHOD YURI LEVIN, MIKHAIL NEDIAK, AND ADI BEN-ISRAEL Abstract. Consider m functions
More informationSSNV231 Sphere digs under internal pressure in great deformations
Titre : SSNV231 Sphère creuse sous pression interne en g[...] Date : 25/02/2014 Page : 1/18 SSNV231 Sphere digs under internal pressure in great deformations Summary: The objective of this test is to validate
More informationStabilization and Acceleration of Algebraic Multigrid Method
Stabilization and Acceleration of Algebraic Multigrid Method Recursive Projection Algorithm A. Jemcov J.P. Maruszewski Fluent Inc. October 24, 2006 Outline 1 Need for Algorithm Stabilization and Acceleration
More informationCVEN 7511 Computational Mechanics of Solids and Structures
CVEN 7511 Computational Mechanics of Solids and Structures Instructor: Kaspar J. Willam Original Version of Class Notes Chishen T. Lin Fall 1990 Chapter 1 Fundamentals of Continuum Mechanics Abstract In
More informationElastoplastic Deformation in a Wedge-Shaped Plate Caused By a Subducting Seamount
Elastoplastic Deformation in a Wedge-Shaped Plate Caused By a Subducting Seamount Min Ding* 1, and Jian Lin 2 1 MIT/WHOI Joint Program, 2 Woods Hole Oceanographic Institution *Woods Hole Oceanographic
More informationLarge bending deformations of pressurized curved tubes
Arch. Mech., 63, 5 6, pp. 57 56, Warszawa Large bending deformations of pressurized curved tubes A. M. KOLESNIKOV Theory of Elasticity Department Southern Federal University Rostov-on-Don, 344, Russian
More informationA return-map algorithm for general associative isotropic elasto-plastic materials in large deformation regimes
International Journal of Plasticity 15 (1999) 1359±1378 return-map algorithm for general associative isotropic elasto-plastic materials in large deformation regimes F. uricchio a, *, R.L. Taylor b a Dipartimento
More informationAdditional corrections to Introduction to Tensor Analysis and the Calculus of Moving Surfaces, Part I.
Additional corrections to Introduction to Tensor Analysis and the Calculus of Moving Surfaces, Part I. 1) Page 42, equation (4.44): All A s in denominator should be a. 2) Page 72, equation (5.84): The
More informationIntroduction to Seismology Spring 2008
MIT OpenCourseWare http://ocw.mit.edu 12.510 Introduction to Seismology Spring 2008 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms. Stress and Strain
More informationChapter 0. Preliminaries. 0.1 Things you should already know
Chapter 0 Preliminaries These notes cover the course MATH45061 (Continuum Mechanics) and are intended to supplement the lectures. The course does not follow any particular text, so you do not need to buy
More information1.050 Engineering Mechanics. Lecture 22: Isotropic elasticity
1.050 Engineering Mechanics Lecture 22: Isotropic elasticity 1.050 Content overview I. Dimensional analysis 1. On monsters, mice and mushrooms 2. Similarity relations: Important engineering tools II. Stresses
More informationMathematical Modeling. of Large Elastic-Plastic Deformations
Applied Mathematical Sciences, Vol. 8, 04, no. 60, 99 996 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/0.988/ams.04.447 Mathematical Modeling of Large Elastic-Plastic Deformations L. U. Sultanov Kazan
More informationMechanics of solids and fluids -Introduction to continuum mechanics
Mechanics of solids and fluids -Introduction to continuum mechanics by Magnus Ekh August 12, 2016 Introduction to continuum mechanics 1 Tensors............................. 3 1.1 Index notation 1.2 Vectors
More informationOn the Rank 1 Convexity of Stored Energy Functions of Physically Linear Stress-Strain Relations
J Elasticity (2007) 86:235 243 DOI 10.1007/s10659-006-9091-z On the Rank 1 Convexity of Stored Energy Functions of Physically Linear Stress-Strain Relations Albrecht Bertram Thomas Böhlke Miroslav Šilhavý
More informationFinite Strain Elastoplasticity: Consistent Eulerian and Lagrangian Approaches
Finite Strain Elastoplasticity: Consistent Eulerian and Lagrangian Approaches by Mohammad Amin Eshraghi A thesis presented to the University of Waterloo in fulfillment of the thesis requirement for the
More information