1 Hooke s law, stiffness, and compliance
|
|
- Grant Barnaby Dalton
- 5 years ago
- Views:
Transcription
1 Non-quilibrium Continuum Physics TA session #5 TA: Yohai Bar Sinai Linear elasticity I This TA session is the first of three at least, maybe more) in which we ll dive deep deep into linear elasticity theory. Linear elasticity by itself can be the topic of a year-long course, and in these lectures we ll try to convey a significant fraction of the richness of the theory. While ran will continue to explore more complex constitutive behaviors non-linear elasticity, visco-elasticity, and so on), in the next few TA sessions we ll stay in the linear realm and try to acquire the knowledge we believe is fundamental for a well-educated physicist. Hooke s law, stiffness, and compliance Hooke s law is σ ij = ε kk δ ij + 2µε ij = + 23 ) µ trε)δ ij + 2µ ε ij ) 3 trε)δ ij ) Let s write it explicitly, to get a better feel of what s going on = 2µε xx + ε xx + ε yy + ε zz ) = 2µ + )ε xx + ε yy + ε zz ) σ yy = 2µε yy + ε xx + ε yy + ε zz ) = 2µ + )ε yy + ε xx + ε zz ) σ zz = 2µε zz + ε xx + ε yy + ε zz ) = 2µ + )ε zz + ε xx + ε yy ) σ ij = 2µε ij, i j 2) The off-diagonal terms shear terms) have a simple dependence, while the diagonal ones are a bit more complicated. Let s try to invert these relations to find εσ) - that is, let s find the compliance matrix for an isotropic linear elastic material. The same considerations that we used to derive q. ) i.e. that ε and σ are related by a 4-rank isotropic tensor) allow us to write ε ij = a σ kk δ ij + b σ ij 3) So finding the compliance reduces to finding a, b. If tr ε = 0, then q.) and 3) reduce to, respectively, σ ij = 2µ ε ij 4) ε ij = b σ ij 5) So we immediately find b = 2µ). Taking the trace of q. )and 3) gives, respectively trσ) = 3 + 2µ) trε) 6) trε) = 3a + b) trσ) 7) This tells us that 3a + b = 3 + 2µ) a = ) µ b = 8) 2µ3 + 2µ) ε ij = 2µ3 + 2µ) trσ)δ ij + 2µ σ ij 9)
2 Writing explicitly, we have ε xx = 2µ3 + 2µ) + ) 2µ ε yy = 2µ3 + 2µ) + ) σ yy 2µ ε zz = 2µ3 + 2µ) + ) σ zz 2µ 2µ3 + 2µ) σ yy + σ zz ) 0) 2µ3 + 2µ) + σ zz ) ) 2µ3 + 2µ) + σ yy ) 2) ε ij = 2µ σ ij i j 3) As discussed in class, the term in parentheses in q.0) is the inverse of the Young s modulus, and for uniaxial stress it reads = /ε xx = 2µ3 + 2µ) + ) 3 + 2µ = µ 2µ + µ. 4) It is the microscopic analogue of the spring constant. If the uniaxial stress is in the x direction, then we have ε yy = 2µ3 + 2µ) σ xx = 2µ3 + 2µ) ε xx = 2 + µ) ε xx, 5) and as discussed in class, ε yy /ε xx is known as the Poisson ratio ν = qs. 0) with these quantities yields a much nicer expression: ε xx = [ νσ yy + σ zz ),. Rewriting 2+µ) ε yy = [σ yy ν + σ zz ), ε zz = [σ zz ν + σ yy ), ε ij = + ν σ ij for i j. 6) 2 Green s function for an infinite medium It seems that the time is ripe to fully and completely solve a problem, with all the 2π s and everything, without resorting to hand waving and scaling arguments. While the emphasis will still be on the structure of the problem, we think it will be instructive, at least once, to write down a problem and solve it exactly. A nice problem to consider is the response of an infinite linear isotropic homogeneous elastic medium to a localized force f = F i δ r), i.e. finding the Green s function of an infinite medium. We define the Green function matrix) G ij r, r 2 ) as the displacement in the i direction at the point r as a response to a localized force in the j direction applied at r 2. For homogeneous materials we know that G ij r, r 2 ) = G ij r r 2 ). We therefore denote 2
3 r = r r 2. You all know well that, within the linear elastic theory, this will allow us to solve the problem of an arbitrary force distribution f r) by convolving f r) with the Green function. Conceptually, the structure is the following. We want to find a displacement field u i r), from which we calculate u i ε ij σ ij j σ ij = δ r) Of course, we will want to do the whole thing backwards. I stress no pun intended ) that we already know how to express divσ) in terms of u i - this is what we called the Navier-Lamé equation. But for completeness, let s do it again: σ ij = ε kk δ ij + 2µ ε ij = k u k δ ij + µ i u j + j u i ) j σ ij = j k u k δ ij + µ i u j + j u i )) = + µ) i j u j + µ j j u i, which is nothing but the u-dependent term of the Navier-Lamé equation. Since the equation is linear, it seems right to solve the problem by Fourier transform. We use the conventions u i q) = d 3 x e i q r u i r), 7) u i r) = d 3 q e i q r u 2π) 3 i q). 8) The equation we want to transform is + µ) j i u i + µ j j u i = F i δ x), 9) which readily gives This is a matrix equation: + µ)q j q i u j µ q j q j u i = F i. 20) [ + µ)q j q i + µ q k q k δ ij u j = F i. 2) Or even more explicitly: + µ) q 2 + µ q 2 + µ) q q 2 + µ) q q 3 u F + µ) q q 2 + µ) q2 2 + µ q 2 + µ) q 2 q 3 u 2 = F 2 22) + µ) q q 3 + µ) q 2 q 3 + µ) q3 2 + µ q 2 u 3 F 3 This matrix can be inverted by using any of your favorite methods, giving u i = [ δij q i q j F µ q k q k 2 ν) q k q k ) 2 j. 23) Where we used ν = the Green function: 2+µ) Just kidding, of course it s intended.. In other words, we have found the Fourier representation of G ij q) = [ δij q i q j. 24) µ q k q k 2 ν) q k q k ) 2 3
4 We now need to preform the inverse Fourier transform. We ll begin with the first term, and do it in spherical coordinates with the q z direction parallel to r: 2π) 3 d 3 e i q x q = q 2 2π) 3 = 2π) 2 e iqr e iqr iqr e iqr cos θ q 2 sin θ dθ dq dφ = q 2 2π) 2 dq = 2 2π) 2 0 sinqr) dq = qr 4πr e iqr cos θ sin θ dθ dq If this result surprises you, maybe you should remind yourself of the first linear field theory that you met in your life - Poisson s equation for a point charge 2 φ = δ r). I ll let you complete the analogy by yourselves. For the second term, we use a dirty trick: { } F qi q j = 2 { )} q k q k ) F q 2 i = i { )} F q j q k q k 2 x i q j q k q k = i { }) ix j F = xj ) 25) 2 x i q k q k 2 x i 4πr = δij 2 4πr x ) ix j 4πr 3 Plugging in 24) we get G ij r) = [ 3 4ν)δ ij + x ix j 6 ν)πµr r 2. 26) A more elegant way to go, is to write G ij as gradients of r I mean r, not r): G ij r) = [ k k r δ ij 8πµr 2 ν) i j r = [ I 2 r 8πµr r 2 ν) 27) 2. Notes about the solution. The displacements go as /r, which means that the strain/stress go as /r 2. Therefore the elastic energy density, which goes like ε 2, goes like /r 4 and its integral diverges. This is much like the case of electrostatics, where the total energy of the electrical field of a point charge diverges. 2. The scaling u /r could also have been obtained from simple dimensional analysis. It is quite common that dimensional considerations in elasticity take the dimension from the shear modulus µ, and then there s an unknown and usually uninteresting) function of ν. 3. You might have noticed that G ij is a symmetric matrix. This might look at first glance as a trivial property that stems from the translational symmetry or rotational symmetry =isotropy), but this is not the case. This symmetry property does not stem from any simple argument that I can think of). Instead, this symmetry is a special case of a more general property that is called reciprocity. For a general 4
5 Figure : Top: Displacement field lines in the x y plane for point F = F ˆx in the horizontal direction. From left to right, with ν = 0, 0.33, Bottom: deformation of a regular mesh under this motion. Note that crossing of two lines of the same color is physically forbidden why?). linear elastic solid, and by general I mean that C ijkl r) can have any symmetry and can even depend on space, the static Green function satisfies G ij r, r ) = G ji r, r) 28) The proof of this general property is given in Sec. 4 of this file. I m not sure I ll have time to go over it in class, but I wanted you to be able to read it even if we don t, in case you re interested. 3 2D elasticity Note: Most of this section appears word by word in ran s lecture notes As shown in class, the field equation of elasticity is the Navier-Lamé equation ρ tt u = + µ) u) + µ 2 u + b. 29) It is notoriously difficult to solve. However, things become much simpler in 2D. There are 3 generic ways in which one can obtain an effectively 2D elastic system, by ignoring the z dimension:. Plane strain: when the system is translationally invariant in z, and therefore z of anything vanishes. This typically holds in very thick in the z direction) systems. 5
6 2. Plane stress: when σ zi = 0. This typically holds in very thin systems in the z direction). 3. Flatland : If the world truly is two dimensional. We will not treat this case as it is a bit delicate. The third case is a bit delicate, and we will not discuss it in the course. We ll now develop the theory for the first two cases. To see how it is done, let us explicitly write Hooke s law 6) in terms of and ν. We note that although the stiffness matrix is a 4-rank tensor, it can be represented as a 6 by 6 matrix by rearranging the entries as σ yy σ zz σ zy = σ zx + ν) 2ν) ν ν ν ν ν ν ν ν ν ν ν ν ε xx ε yy ε zz ε zy ε zx ε xy. 30) We can invert this relation into the form ε = Sσ, where S = C is the compliance tensor if you noticed that C is called the stiffness tensor and S is called the compliance tensor and wondered about it, this is not a mistake and there is no intention to confuse you. It is a long-time convention that cannot be reverted anymore). We write the last relation explicitly as ε xx ν ν ε yy ν ν σ yy ν ν ε zz ε zy ε zx ε xy = ν ν ν We are now ready to perform the reduction to 2D. 3. Plane-stress σ zz σ zy σ zx. 3) We first consider objects that are thin in one dimension, say z, and are deformed in the xy-plane. What happens in the z-direction? Since the two planes z = 0 and z = h where h is the thickness which is much smaller than any other lengthscale in the problem) are traction-free, we approximate σ zz = 0 everywhere an approximation that becomes better and better as h 0). Similarly, we have σ zy = σ zx = 0. We can therefore set σ zz = σ zy = σ zx = 0 in q. 3) to obtain ε xx ε yy = ν 0 ν 0 ε xy ν σ yy 32) and ε zz = ν + σ yy ). 33) 6
7 To obtain the plane-stress analog of q. the Navier-Lamé equation, we need to invert q. 32), obtaining σ yy = ν 0 ε xx ν 0 ε ν σ 2 yy. 34) xy 0 0 ν ε xy Note that the last equation can not be obtained from q. 30) by simply removing columns and rows. We can now substitute the last relation in the 2D momentum balance equation σ = ρ tt u we stress again that σ and u are already 2D here). The resulting 2D equation reads [ [ ν ν + u) ν) 2 + ν) which is identical in form to q. 29) simply with a renormalized = The shear modulus µ remains unchanged 2 u = ρ tt u, 35) ν ν 2 = 2νµ ν = 2µ + 2µ. 36) µ = µ = 2 + ν). 37) Finally, we can substitute x, y) and σ yy x, y) inside q. 33) to obtain ε zz x, y). Note that u z x, y, z) = ε zz x, y)z is linear in z. 3.2 Plane-strain We now consider objects that are very thick in one dimension, say z, and are deformed in the xy-plane with no z dependence. These physical conditions are termed plane-strain and are characterized by ɛ zx = ɛ zy = ɛ zz = 0. liminating these components from q. 30) we obtain and σ yy = σ zz x, y) = + ν) 2ν) ν ν 0 ν ν ν ε xx ε yy ε xy 38) ν + ν) 2ν) [ε xxx, y) + ε yy x, y). 39) We can now substitute q. 38) in the 2D momentum balance equation σ = ρ tt u where again σ and u are 2D). The resulting 2D equation is identical to q. 29), both in form and in the elastic constants. With the solution at hand, we can use q. 39) to calculate σ zz x, y). Finally, we note that q. 38) can be inverted to ε xx ε yy = + ν ε xy ν ν 0 ν ν σ yy, 40) 7
8 which can not be simply obtained from q. 3) by eliminating columns and rows. Using the last relation we can rewrite q. 39) as σ zz x, y) = ν [ x, y) + σ yy x, y). 4) In summary, we see that in both plane-stress and plane-strain cases we can work with 2D objects instead of their 3D counterparts, which is a significant simplification. This allows to use the heavy mathematical tools available in 2D: complex analysis and conformal mapping. One has, though, to be careful with the elastic constants as explained above. 4 Clapeyron s & Betti s theorems and reciprocity time permitting) Clapeyron s theorem is a nice result in static elasticity that says that the total elastic energy stored in a body is equal to half the work done by the external forces computed assuming these forces had remained constant from the initial state to the final state. You already know a trivial version of this theorem, in the context of the energy stored in a simple spring. Assume you load a given spring with a force F. The energy stored in the spring is U = 2 kx2 = 2 kx)x = 2 F x. 42) If you assume that during the elongation of the spring the force had the constant value F although the force clearly started from zero and ramped up to F ) you get that the total work done was F x. The real energy, however, is exactly one half that value. The proof is pretty simple. Assume an elastic body is subject to volume forces br) and surface tractions tr). The total elastic energy in the body is U = [ 2) = 2 = 2 σ 2 ijr)ε ij r)d 3 r ) σ ij u j n i d 2 r + σ ij i u j d 3 r = 2 [ b j u j d 3 r = 2 [ i σ ij u j ) u j i σ ij d 3 r 43) b j u j d 3 r t j u j d 2 r + In the transition ) we used the symmetry of σ and in 2) we used Gauss theorem and the equilibrium condition j σ ij + b i = 0. Betti s theorem sometimes called Betti s reciprocal theorem ) is a general important result about the energetics in static elasticity. It is useful in theoretical analyses e.g. for proving the uniqueness of solutions to the static Navier-Lamè equation) and is also fo practical use if average quantities are calculated you might have an exercise about that in the HW). Suppose that when a set of body and traction forces b ) and t ) is applied to a body, the resulting deformation field is u ), and the stress and strain fields are σ ) and ε ). Suppose also that when a different set of body and traction forces, b 2) and t 2), is applied then the resulting deformation field is u 2), and the stress and strain fields are σ 2) and ε 2). Betti s theorem states that t ) i u 2) i d 2 r + b ) i u 2) i d 3 r = 8 t 2) i u ) i d 2 r + b 2) i u ) i d 3 r. 44)
9 That is, the work done by the set ) through the displacements produced by the set 2) is equal to the work done by the set 2) through the displacements produced by the set ). To see this, consider the total energy of the system if both sets of forces are applied: U = 2 = 2 ) ) σ ) ij + σ 2) ij ε ) ij + ε 2) ij d 3 r 45) [ σ ) ij ε) ij + σ 2) ij ε2) ij + σ ) ij ε2) ij + σ 2) ij ε) ij d 3 r 46) [ = U ) + U 2) + 2 σ ) ij ε2) ij }{{} A + σ 2) ij ε) ij }{{} B d 3 r. 47) The first two terms are the energies of each of the pure modes and the integral is the interaction energy. Using exactly the same arguments as we used in q. 43) it is possible and easy) to show that term A equals the left-hand-side of q. 44) and term B equals the right-hand-side. However, terms A and B are equal since A = σ ) ij ε2) ij = C ijkl ε ) kl ε2) ij = C klij ε ) kl ε2) ij = ε ) kl σ2) kl = B, 48) where the symmetry C ijkl = C klij was used, and the theorem is proven. Note that we did not invoke isotropy or translational invariance. The symmetry C ijkl = C klij is a thermodynamic one that holds for the most general linear elastic case. With Betti s theorem it is very easy to prove the reciprocity property, q. 28). If the set ) is a point force F ) applied at r ) and set 2) is a point force F 2) applied at r 2) then b ) i = F ) i δr r ) ) b 2) i = F 2) i δr r 2) ) 49) u ) i r) = G ij r, r ) )F ) j u 2) i r) = G ij r, r 2) )F 2) j 50) Homogeneous boundary conditions on mean either u = 0 or t = 0, or that some portions of the surface is subject to that and other portions to the other. In any case, this means that the boundary integrals in q. 44) vanish identically. Then, using the properties of the delta function, we immediately get F ) i F 2) j G ij r ), r 2) ) = F 2) i F ) j G ij r 2), r ) ) 5) Since F, F 2, r ) and r 2) were arbitrary the reciprocal theorem of q. 28) is obtained. 9
PEAT 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 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 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 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 information3D Elasticity Theory
3D lasticity Theory Many structural analysis problems are analysed using the theory of elasticity in which Hooke s law is used to enforce proportionality between stress and strain at any deformation level.
More information3D and Planar Constitutive Relations
3D and Planar Constitutive Relations A School on Mechanics of Fibre Reinforced Polymer Composites Knowledge Incubation for TEQIP Indian Institute of Technology Kanpur PM Mohite Department of Aerospace
More informationMechanics of Earthquakes and Faulting
Mechanics of Earthquakes and Faulting www.geosc.psu.edu/courses/geosc508 Overview Milestones in continuum mechanics Concepts of modulus and stiffness. Stress-strain relations Elasticity Surface and body
More informationChapter 1. Continuum mechanics review. 1.1 Definitions and nomenclature
Chapter 1 Continuum mechanics review We will assume some familiarity with continuum mechanics as discussed in the context of an introductory geodynamics course; a good reference for such problems is Turcotte
More informationMacroscopic theory Rock as 'elastic continuum'
Elasticity and Seismic Waves Macroscopic theory Rock as 'elastic continuum' Elastic body is deformed in response to stress Two types of deformation: Change in volume and shape Equations of motion Wave
More informationFinite Element Method in Geotechnical Engineering
Finite Element Method in Geotechnical Engineering Short Course on + Dynamics Boulder, Colorado January 5-8, 2004 Stein Sture Professor of Civil Engineering University of Colorado at Boulder Contents Steps
More informationMECH 5312 Solid Mechanics II. Dr. Calvin M. Stewart Department of Mechanical Engineering The University of Texas at El Paso
MECH 5312 Solid Mechanics II Dr. Calvin M. Stewart Department of Mechanical Engineering The University of Texas at El Paso Table of Contents Thermodynamics Derivation Hooke s Law: Anisotropic Elasticity
More information3.2 Hooke s law anisotropic elasticity Robert Hooke ( ) Most general relationship
3.2 Hooke s law anisotropic elasticity Robert Hooke (1635-1703) Most general relationship σ = C ε + C ε + C ε + C γ + C γ + C γ 11 12 yy 13 zz 14 xy 15 xz 16 yz σ = C ε + C ε + C ε + C γ + C γ + C γ yy
More informationFundamentals of Linear Elasticity
Fundamentals of Linear Elasticity Introductory Course on Multiphysics Modelling TOMASZ G. ZIELIŃSKI bluebox.ippt.pan.pl/ tzielins/ Institute of Fundamental Technological Research of the Polish Academy
More information1 Gauss integral theorem for tensors
Non-Equilibrium Continuum Physics TA session #1 TA: Yohai Bar Sinai 16.3.216 Index Gymnastics: Gauss Theorem, Isotropic Tensors, NS Equations The purpose of today s TA session is to mess a bit with tensors
More informationNDT&E Methods: UT. VJ Technologies CAVITY INSPECTION. Nondestructive Testing & Evaluation TPU Lecture Course 2015/16.
CAVITY INSPECTION NDT&E Methods: UT VJ Technologies NDT&E Methods: UT 6. NDT&E: Introduction to Methods 6.1. Ultrasonic Testing: Basics of Elasto-Dynamics 6.2. Principles of Measurement 6.3. The Pulse-Echo
More informationBasic Equations of Elasticity
A Basic Equations of Elasticity A.1 STRESS The state of stress at any point in a loaded bo is defined completely in terms of the nine components of stress: σ xx,σ yy,σ zz,σ xy,σ yx,σ yz,σ zy,σ zx,andσ
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 Introduction to mechanics
21 Motivation Thermodynamic bodies are being characterized by two competing opposite phenomena, energy and entropy which some researchers in thermodynamics would classify as cause and chance or determinancy
More informationConservation of mass. Continuum Mechanics. Conservation of Momentum. Cauchy s Fundamental Postulate. # f body
Continuum Mechanics We ll stick with the Lagrangian viewpoint for now Let s look at a deformable object World space: points x in the object as we see it Object space (or rest pose): points p in some reference
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 informationUnderstand basic stress-strain response of engineering materials.
Module 3 Constitutive quations Learning Objectives Understand basic stress-strain response of engineering materials. Quantify the linear elastic stress-strain response in terms of tensorial quantities
More informationLecture 8: Tissue Mechanics
Computational Biology Group (CoBi), D-BSSE, ETHZ Lecture 8: Tissue Mechanics Prof Dagmar Iber, PhD DPhil MSc Computational Biology 2015/16 7. Mai 2016 2 / 57 Contents 1 Introduction to Elastic Materials
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 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 informationStress/Strain. Outline. Lecture 1. Stress. Strain. Plane Stress and Plane Strain. Materials. ME EN 372 Andrew Ning
Stress/Strain Lecture 1 ME EN 372 Andrew Ning aning@byu.edu Outline Stress Strain Plane Stress and Plane Strain Materials otes and News [I had leftover time and so was also able to go through Section 3.1
More informationContinuum Mechanics. Continuum Mechanics and Constitutive Equations
Continuum Mechanics Continuum Mechanics and Constitutive Equations Continuum mechanics pertains to the description of mechanical behavior of materials under the assumption that the material is a uniform
More informationMechanics of Earthquakes and Faulting
Mechanics of Earthquakes and Faulting www.geosc.psu.edu/courses/geosc508 Standard Solids and Fracture Fluids: Mechanical, Chemical Effects Effective Stress Dilatancy Hardening and Stability Mead, 1925
More information3 2 6 Solve the initial value problem u ( t) 3. a- If A has eigenvalues λ =, λ = 1 and corresponding eigenvectors 1
Math Problem a- If A has eigenvalues λ =, λ = 1 and corresponding eigenvectors 1 3 6 Solve the initial value problem u ( t) = Au( t) with u (0) =. 3 1 u 1 =, u 1 3 = b- True or false and why 1. if A is
More informationRock Rheology GEOL 5700 Physics and Chemistry of the Solid Earth
Rock Rheology GEOL 5700 Physics and Chemistry of the Solid Earth References: Turcotte and Schubert, Geodynamics, Sections 2.1,-2.4, 2.7, 3.1-3.8, 6.1, 6.2, 6.8, 7.1-7.4. Jaeger and Cook, Fundamentals of
More informationVariational principles in mechanics
CHAPTER Variational principles in mechanics.1 Linear Elasticity n D Figure.1: A domain and its boundary = D [. Consider a domain Ω R 3 with its boundary = D [ of normal n (see Figure.1). The problem of
More informationLecture 7. Properties of Materials
MIT 3.00 Fall 2002 c W.C Carter 55 Lecture 7 Properties of Materials Last Time Types of Systems and Types of Processes Division of Total Energy into Kinetic, Potential, and Internal Types of Work: Polarization
More informationGroup Representations
Group Representations Alex Alemi November 5, 2012 Group Theory You ve been using it this whole time. Things I hope to cover And Introduction to Groups Representation theory Crystallagraphic Groups Continuous
More informationReview of Strain Energy Methods and Introduction to Stiffness Matrix Methods of Structural Analysis
uke University epartment of Civil and Environmental Engineering CEE 42L. Matrix Structural Analysis Henri P. Gavin Fall, 22 Review of Strain Energy Methods and Introduction to Stiffness Matrix Methods
More informationElements of Rock Mechanics
Elements of Rock Mechanics Stress and strain Creep Constitutive equation Hooke's law Empirical relations Effects of porosity and fluids Anelasticity and viscoelasticity Reading: Shearer, 3 Stress Consider
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 informationSurface force on a volume element.
STRESS and STRAIN Reading: Section. of Stein and Wysession. In this section, we will see how Newton s second law and Generalized Hooke s law can be used to characterize the response of continuous medium
More informationLecture 8. Stress Strain in Multi-dimension
Lecture 8. Stress Strain in Multi-dimension Module. General Field Equations General Field Equations [] Equilibrium Equations in Elastic bodies xx x y z yx zx f x 0, etc [2] Kinematics xx u x x,etc. [3]
More informationBasic concepts to start Mechanics of Materials
Basic concepts to start Mechanics of Materials Georges Cailletaud Centre des Matériaux Ecole des Mines de Paris/CNRS Notations Notations (maths) (1/2) A vector v (element of a vectorial space) can be seen
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 informationNumerical Modelling in Geosciences. Lecture 6 Deformation
Numerical Modelling in Geosciences Lecture 6 Deformation Tensor Second-rank tensor stress ), strain ), strain rate ) Invariants quantities independent of the coordinate system): - First invariant trace:!!
More information3.22 Mechanical Properties of Materials Spring 2008
MIT OpenCourseWare http://ocw.mit.edu 3.22 Mechanical Properties of Materials Spring 2008 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms. Quiz #1 Example
More informationConstitutive Relations
Constitutive Relations Dr. Andri Andriyana 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 field
More informationElasticity in two dimensions 1
Elasticity in two dimensions 1 Elasticity in two dimensions Chapters 3 and 4 of Mechanics of the Cell, as well as its Appendix D, contain selected results for the elastic behavior of materials in two and
More informationExercise: concepts from chapter 8
Reading: Fundamentals of Structural Geology, Ch 8 1) The following exercises explore elementary concepts associated with a linear elastic material that is isotropic and homogeneous with respect to elastic
More informationOther state variables include the temperature, θ, and the entropy, S, which are defined below.
Chapter 3 Thermodynamics In order to complete the formulation we need to express the stress tensor T and the heat-flux vector q in terms of other variables. These expressions are known as constitutive
More informationElectromagnetism HW 1 math review
Electromagnetism HW math review Problems -5 due Mon 7th Sep, 6- due Mon 4th Sep Exercise. The Levi-Civita symbol, ɛ ijk, also known as the completely antisymmetric rank-3 tensor, has the following properties:
More informationMechanics of Biomaterials
Mechanics of Biomaterials Lecture 7 Presented by Andrian Sue AMME498/998 Semester, 206 The University of Sydney Slide Mechanics Models The University of Sydney Slide 2 Last Week Using motion to find forces
More informationLinear Elasticity ( ) Objectives. Equipment. Introduction. ε is then
Linear Elasticity Objectives In this lab you will measure the Young s Modulus of a steel wire. In the process, you will gain an understanding of the concepts of stress and strain. Equipment Young s Modulus
More informationStress, Strain, Mohr s Circle
Stress, Strain, Mohr s Circle The fundamental quantities in solid mechanics are stresses and strains. In accordance with the continuum mechanics assumption, the molecular structure of materials is neglected
More informationMechanical Properties of Materials
Mechanical Properties of Materials Strains Material Model Stresses Learning objectives Understand the qualitative and quantitative description of mechanical properties of materials. Learn the logic of
More informationELASTICITY (MDM 10203)
LASTICITY (MDM 10203) Lecture Module 5: 3D Constitutive Relations Dr. Waluyo Adi Siswanto University Tun Hussein Onn Malaysia Generalised Hooke's Law In one dimensional system: = (basic Hooke's law) Considering
More informationGame Physics. Game and Media Technology Master Program - Utrecht University. Dr. Nicolas Pronost
Game and Media Technology Master Program - Utrecht University Dr. Nicolas Pronost Soft body physics Soft bodies In reality, objects are not purely rigid for some it is a good approximation but if you hit
More informationDynamics of the Mantle and Lithosphere ETH Zürich Continuum Mechanics in Geodynamics: Equation cheat sheet
Dynamics of the Mantle and Lithosphere ETH Zürich Continuum Mechanics in Geodynamics: Equation cheat sheet or all equations you will probably ever need Definitions 1. Coordinate system. x,y,z or x 1,x
More informationSoft Bodies. Good approximation for hard ones. approximation breaks when objects break, or deform. Generalization: soft (deformable) bodies
Soft-Body Physics Soft Bodies Realistic objects are not purely rigid. Good approximation for hard ones. approximation breaks when objects break, or deform. Generalization: soft (deformable) bodies Deformed
More information16.20 HANDOUT #2 Fall, 2002 Review of General Elasticity
6.20 HANDOUT #2 Fall, 2002 Review of General Elasticity NOTATION REVIEW (e.g., for strain) Engineering Contracted Engineering Tensor Tensor ε x = ε = ε xx = ε ε y = ε 2 = ε yy = ε 22 ε z = ε 3 = ε zz =
More informationGenerators for Continuous Coordinate Transformations
Page 636 Lecture 37: Coordinate Transformations: Continuous Passive Coordinate Transformations Active Coordinate Transformations Date Revised: 2009/01/28 Date Given: 2009/01/26 Generators for Continuous
More information12. Stresses and Strains
12. Stresses and Strains Finite Element Method Differential Equation Weak Formulation Approximating Functions Weighted Residuals FEM - Formulation Classification of Problems Scalar Vector 1-D T(x) u(x)
More informationChapter 1 Fracture Physics
Chapter 1 Fracture Physics In this chapter we give a brief account of the theory of fracture, restricting ourselves to brittle fracture. The structure of this chapter is as follows. The first part describes
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 informationLinearized theory of elasticity
Linearized theory of elasticity Arie Verhoeven averhoev@win.tue.nl CASA Seminar, May 24, 2006 Seminar: Continuum mechanics 1 Stress and stress principles Bart Nowak March 8 2 Strain and deformation Mark
More informationC:\Users\whit\Desktop\Active\304_2012_ver_2\_Notes\4_Torsion\1_torsion.docx 6
C:\Users\whit\Desktop\Active\304_2012_ver_2\_Notes\4_Torsion\1_torsion.doc 6 p. 1 of Torsion of circular bar The cross-sections rotate without deformation. The deformation that does occur results from
More information2.2 Relation Between Mathematical & Engineering Constants Isotropic Materials Orthotropic Materials
Chapter : lastic Constitutive quations of a Laminate.0 Introduction quations of Motion Symmetric of Stresses Tensorial and ngineering Strains Symmetry of Constitutive quations. Three-Dimensional Constitutive
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 informationChapter 9: Differential Analysis of Fluid Flow
of Fluid Flow Objectives 1. Understand how the differential equations of mass and momentum conservation are derived. 2. Calculate the stream function and pressure field, and plot streamlines for a known
More informationMAE 323: Lecture 1. Review
This review is divided into two parts. The first part is a mini-review of statics and solid mechanics. The second part is a review of matrix/vector fundamentals. The first part is given as an refresher
More informationPart 8: Rigid Body Dynamics
Document that contains homework problems. Comment out the solutions when printing off for students. Part 8: Rigid Body Dynamics Problem 1. Inertia review Find the moment of inertia for a thin uniform rod
More informationInverse Design (and a lightweight introduction to the Finite Element Method) Stelian Coros
Inverse Design (and a lightweight introduction to the Finite Element Method) Stelian Coros Computational Design Forward design: direct manipulation of design parameters Level of abstraction Exploration
More information20. Rheology & Linear Elasticity
I Main Topics A Rheology: Macroscopic deformation behavior B Linear elasticity for homogeneous isotropic materials 10/29/18 GG303 1 Viscous (fluid) Behavior http://manoa.hawaii.edu/graduate/content/slide-lava
More informationChapter 9: Differential Analysis
9-1 Introduction 9-2 Conservation of Mass 9-3 The Stream Function 9-4 Conservation of Linear Momentum 9-5 Navier Stokes Equation 9-6 Differential Analysis Problems Recall 9-1 Introduction (1) Chap 5: Control
More informationDistributed: Wednesday, March 17, 2004
MASSACHUSETTS INSTITUTE OF TECHNOLOGY DEPARTMENT OF MECHANICAL ENGINEERING CAMBRIDGE, MASSACHUSETTS 019.00 MECHANICS AND MATERIALS II QUIZ I SOLUTIONS Distributed: Wednesday, March 17, 004 This quiz consists
More informationChapter 2. Rubber Elasticity:
Chapter. Rubber Elasticity: The mechanical behavior of a rubber band, at first glance, might appear to be Hookean in that strain is close to 100% recoverable. However, the stress strain curve for a rubber
More informationIntroduction, Basic Mechanics 2
Computational Biomechanics 18 Lecture : Introduction, Basic Mechanics Ulli Simon, Lucas Engelhardt, Martin Pietsch Scientific Computing Centre Ulm, UZWR Ulm University Contents Mechanical Basics Moment
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 information1 Exercise: Linear, incompressible Stokes flow with FE
Figure 1: Pressure and velocity solution for a sinking, fluid slab impinging on viscosity contrast problem. 1 Exercise: Linear, incompressible Stokes flow with FE Reading Hughes (2000), sec. 4.2-4.4 Dabrowski
More informationPart 5 ACOUSTIC WAVE PROPAGATION IN ANISOTROPIC MEDIA
Part 5 ACOUSTIC WAVE PROPAGATION IN ANISOTROPIC MEDIA Review of Fundamentals displacement-strain relation stress-strain relation balance of momentum (deformation) (constitutive equation) (Newton's Law)
More informationTHE GENERAL ELASTICITY PROBLEM IN SOLIDS
Chapter 10 TH GNRAL LASTICITY PROBLM IN SOLIDS In Chapters 3-5 and 8-9, we have developed equilibrium, kinematic and constitutive equations for a general three-dimensional elastic deformable solid bod.
More informationMICROMECHANICS AND HOMOGENIZATION. Materials Containing Coated and Uncoated Spherical Inhomogeneities
MICROMECHANICS AND HOMOGENIZATION Materials Containing Coated and Uncoated Spherical Inhomogeneities Václav Nežerka, 2012 i Acknowledgement I would like to thank my supervisor, Jan Zeman, who supported
More informationNIELINIOWA OPTYKA MOLEKULARNA
NIELINIOWA OPTYKA MOLEKULARNA chapter 1 by Stanisław Kielich translated by:tadeusz Bancewicz http://zon8.physd.amu.edu.pl/~tbancewi Poznan,luty 2008 ELEMENTS OF THE VECTOR AND TENSOR ANALYSIS Reference
More informationMODELING OF CONCRETE MATERIALS AND STRUCTURES. Kaspar Willam
MODELING OF CONCRETE MATERIALS AND STRUCTURES Class Meeting #1: Fundamentals Kaspar Willam University of Colorado at Boulder Notation: Direct and indicial tensor formulations Fundamentals: Stress and Strain
More informationSTEP Support Programme. STEP 2 Matrices Topic Notes
STEP Support Programme STEP 2 Matrices Topic Notes Definitions............................................. 2 Manipulating Matrices...................................... 3 Transformations.........................................
More informationMATERIAL MECHANICS, SE2126 COMPUTER LAB 4 MICRO MECHANICS. E E v E E E E E v E E + + = m f f. f f
MATRIAL MCHANICS, S226 COMPUTR LAB 4 MICRO MCHANICS 2 2 2 f m f f m T m f m f f m v v + + = + PART A SPHRICAL PARTICL INCLUSION Consider a solid granular material, a so called particle composite, shown
More informationUniformity of the Universe
Outline Universe is homogenous and isotropic Spacetime metrics Friedmann-Walker-Robertson metric Number of numbers needed to specify a physical quantity. Energy-momentum tensor Energy-momentum tensor of
More informationLecture Note 1. Introduction to Elasticity Equations
ME34B Elasticity of Microscopic tructures tanford University Winter 24 Lecture Note. Introduction to Elasticity Equations Chris Weinberger and Wei Cai c All rights reserved January 2, 24 Contents Index
More informationLinear Algebra March 16, 2019
Linear Algebra March 16, 2019 2 Contents 0.1 Notation................................ 4 1 Systems of linear equations, and matrices 5 1.1 Systems of linear equations..................... 5 1.2 Augmented
More informationCH.6. LINEAR ELASTICITY. Multimedia Course on Continuum Mechanics
CH.6. LINEAR ELASTICITY Multimedia Course on Continuum Mechanics Overview Hypothesis of the Linear Elasticity Theory Linear Elastic Constitutive Equation Generalized Hooke s Law Elastic Potential Isotropic
More informationExample 3.7 Consider the undeformed configuration of a solid as shown in Figure 3.60.
162 3. The linear 3-D elasticity mathematical model The 3-D elasticity model is of great importance, since it is our highest order hierarchical model assuming linear elastic behavior. Therefore, it provides
More information202 Index. failure, 26 field equation, 122 force, 1
Index acceleration, 12, 161 admissible function, 155 admissible stress, 32 Airy's stress function, 122, 124 d'alembert's principle, 165, 167, 177 amplitude, 171 analogy, 76 anisotropic material, 20 aperiodic
More informationChapter 2 Governing Equations
Chapter Governing Equations Abstract In this chapter fundamental governing equations for propagation of a harmonic disturbance on the surface of an elastic half-space is presented. The elastic media is
More informationExact Solutions of the Einstein Equations
Notes from phz 6607, Special and General Relativity University of Florida, Fall 2004, Detweiler Exact Solutions of the Einstein Equations These notes are not a substitute in any manner for class lectures.
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 informationDROP-WEIGHT SYSTEM FOR DYNAMIC PRESSURE CALIBRATION
DROP-WIGHT SYSTM FOR DYNAMIC PRSSUR CALIBRATION Antti Lakka and Sari Saxholm MIKS, Centre for Metrology and Accreditation P.O. Box 9 (Tekniikantie 1) FI-2151 spoo Finland telephone: +358-29-554456, e-mail:
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 informationIntroduction to Seismology Spring 2008
MIT OpenCourseWare http://ocw.mit.edu 1.510 Introduction to Seismology Spring 008 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms. 1.510 Introduction to
More informationBasic Concepts of Strain and Tilt. Evelyn Roeloffs, USGS June 2008
Basic Concepts of Strain and Tilt Evelyn Roeloffs, USGS June 2008 1 Coordinates Right-handed coordinate system, with positions along the three axes specified by x,y,z. x,y will usually be horizontal, and
More informationAnalytical Mechanics: Elastic Deformation
Analytical Mechanics: Elastic Deformation Shinichi Hirai Dept. Robotics, Ritsumeikan Univ. Shinichi Hirai (Dept. Robotics, Ritsumeikan Univ.) Analytical Mechanics: Elastic Deformation 1 / 60 Agenda Agenda
More informationSEISMOLOGY I. Laurea Magistralis in Physics of the Earth and of the Environment. Elasticity. Fabio ROMANELLI
SEISMOLOGY I Laurea Magistralis in Physics of the Earth and of the Environment Elasticity Fabio ROMANELLI Dept. Earth Sciences Università degli studi di Trieste romanel@dst.units.it 1 Elasticity and Seismic
More informationFinite Element Analysis Prof. Dr. B. N. Rao Department of Civil Engineering Indian Institute of Technology, Madras. Lecture - 06
Finite Element Analysis Prof. Dr. B. N. Rao Department of Civil Engineering Indian Institute of Technology, Madras Lecture - 06 In the last lecture, we have seen a boundary value problem, using the formal
More informationMASSACHUSETTS INSTITUTE OF TECHNOLOGY DEPARTMENT OF MATERIALS SCIENCE AND ENGINEERING CAMBRIDGE, MASSACHUSETTS 02139
MASSACHUSETTS NSTTUTE OF TECHNOLOGY DEPARTMENT OF MATERALS SCENCE AND ENGNEERNG CAMBRDGE, MASSACHUSETTS 39 3. MECHANCAL PROPERTES OF MATERALS PROBLEM SET SOLUTONS Reading Ashby, M.F., 98, Tensors: Notes
More informationEE C245 ME C218 Introduction to MEMS Design Fall 2007
EE C245 ME C218 Introduction to MEMS Design Fall 2007 Prof. Clark T.-C. Nguyen Dept. of Electrical Engineering & Computer Sciences University of California at Berkeley Berkeley, CA 94720 Lecture 13: Material
More informationMECH 5312 Solid Mechanics II. Dr. Calvin M. Stewart Department of Mechanical Engineering The University of Texas at El Paso
MECH 5312 Solid Mechanics II Dr. Calvin M. Stewart Department of Mechanical Engineering The University of Texas at El Paso Table of Contents Preliminary Math Concept of Stress Stress Components Equilibrium
More information