Physics of Continuous media

Similar documents
PEAT SEISMOLOGY Lecture 2: Continuum mechanics

Continuum Mechanics. Continuum Mechanics and Constitutive Equations

Elements of Rock Mechanics

3D Elasticity Theory

Chapter 16: Elastic Solids

Lecture 8: Tissue Mechanics

Mathematical Background

Fundamentals of Linear Elasticity

MECH 5312 Solid Mechanics II. Dr. Calvin M. Stewart Department of Mechanical Engineering The University of Texas at El Paso

Chapter 12. Static Equilibrium and Elasticity

Stress Strain Elasticity Modulus Young s Modulus Shear Modulus Bulk Modulus. Case study

Introduction to Seismology Spring 2008

Chapter 1. Continuum mechanics review. 1.1 Definitions and nomenclature

Macroscopic theory Rock as 'elastic continuum'

Game Physics. Game and Media Technology Master Program - Utrecht University. Dr. Nicolas Pronost

20. Rheology & Linear Elasticity

A short review of continuum mechanics

Lecture 8. Stress Strain in Multi-dimension

Continuum mechanics V. Constitutive equations. 1. Constitutive equation: definition and basic axioms

Elasticity in two dimensions 1

Chapter 3 Stress, Strain, Virtual Power and Conservation Principles

3.22 Mechanical Properties of Materials Spring 2008

SEMM Mechanics PhD Preliminary Exam Spring Consider a two-dimensional rigid motion, whose displacement field is given by

Finite Element Method in Geotechnical Engineering

Statics. Phys101 Lectures 19,20. Key points: The Conditions for static equilibrium Solving statics problems Stress and strain. Ref: 9-1,2,3,4,5.

3 2 6 Solve the initial value problem u ( t) 3. a- If A has eigenvalues λ =, λ = 1 and corresponding eigenvectors 1

Mechanics PhD Preliminary Spring 2017

Continuum Mechanics Fundamentals

Basic Concepts of Strain and Tilt. Evelyn Roeloffs, USGS June 2008

Soft Bodies. Good approximation for hard ones. approximation breaks when objects break, or deform. Generalization: soft (deformable) bodies

Surface force on a volume element.

By drawing Mohr s circle, the stress transformation in 2-D can be done graphically. + σ x σ y. cos 2θ + τ xy sin 2θ, (1) sin 2θ + τ xy cos 2θ.

Basic Equations of Elasticity

16.21 Techniques of Structural Analysis and Design Spring 2003 Unit #5 - Constitutive Equations

Problem " Â F y = 0. ) R A + 2R B + R C = 200 kn ) 2R A + 2R B = 200 kn [using symmetry R A = R C ] ) R A + R B = 100 kn

Lecture 3: 1. Lecture 3.

Mechanics of Earthquakes and Faulting

Exercise: concepts from chapter 8

Understand basic stress-strain response of engineering materials.

CHAPTER 12 STATIC EQUILIBRIUM AND ELASTICITY. Conditions for static equilibrium Center of gravity (weight) Examples of static equilibrium

Chapter Two: Mechanical Properties of materials

Constitutive Equations

- Marine Hydrodynamics. Lecture 4. Knowns Equations # Unknowns # (conservation of mass) (conservation of momentum)

Engineering Sciences 241 Advanced Elasticity, Spring Distributed Thursday 8 February.

Symmetry and Properties of Crystals (MSE638) Stress and Strain Tensor

Basic concepts to start Mechanics of Materials

Fundamentals of Fluid Dynamics: Elementary Viscous Flow

Elastic Wave Theory. LeRoy Dorman Room 436 Ritter Hall Tel: Based on notes by C. R. Bentley. Version 1.

Mechanical Properties of Materials

PHYSICS. Course Structure. Unit Topics Marks. Physical World and Measurement. 1 Physical World. 2 Units and Measurements.

Tensor Transformations and the Maximum Shear Stress. (Draft 1, 1/28/07)

In this section, thermoelasticity is considered. By definition, the constitutive relations for Gradθ. This general case

SEISMOLOGY I. Laurea Magistralis in Physics of the Earth and of the Environment. Elasticity. Fabio ROMANELLI

NDT&E Methods: UT. VJ Technologies CAVITY INSPECTION. Nondestructive Testing & Evaluation TPU Lecture Course 2015/16.

Continuum mechanism: Stress and strain

3.2 Hooke s law anisotropic elasticity Robert Hooke ( ) Most general relationship

Mechanics of Biomaterials

Concept Question Comment on the general features of the stress-strain response under this loading condition for both types of materials

2.3 BASIC EQUATIONS OF CONTINUUM MECHANICS

Theories of Straight Beams

Variational principles in mechanics

Lecture contents. Stress and strain Deformation potential. NNSE 618 Lecture #23

DIVIDED SYLLABUS ( ) - CLASS XI PHYSICS (CODE 042) COURSE STRUCTURE APRIL

Lecture 20: Modes in a crystal and continuum

Review of Strain Energy Methods and Introduction to Stiffness Matrix Methods of Structural Analysis

Elasticity. A PowerPoint Presentation by Paul E. Tippens, Professor of Physics Southern Polytechnic State University Modified by M.

Losses mechanisms in quartz tuning forks vibrating at the fundamental and first overtone mode

Rock Rheology GEOL 5700 Physics and Chemistry of the Solid Earth

Getting started: CFD notation

Lecture #6: 3D Rate-independent Plasticity (cont.) Pressure-dependent plasticity

Elements of Continuum Elasticity. David M. Parks Mechanics and Materials II February 25, 2004

Conservation of mass. Continuum Mechanics. Conservation of Momentum. Cauchy s Fundamental Postulate. # f body

MASSACHUSETTS INSTITUTE OF TECHNOLOGY DEPARTMENT OF MATERIALS SCIENCE AND ENGINEERING CAMBRIDGE, MASSACHUSETTS 02139

PURE BENDING. If a simply supported beam carries two point loads of 10 kn as shown in the following figure, pure bending occurs at segment BC.

Elasticity: Term Paper. Danielle Harper. University of Central Florida

Chapter 2. Rubber Elasticity:

Mechanics of Earthquakes and Faulting

Dynamics of the Mantle and Lithosphere ETH Zürich Continuum Mechanics in Geodynamics: Equation cheat sheet

STANDARD SAMPLE. Reduced section " Diameter. Diameter. 2" Gauge length. Radius

Linearized theory of elasticity

INTRODUCTION TO STRAIN

Lecture 18. In other words, if you double the stress, you double the resulting strain.

Chapter 1 Fluid Characteristics

Chapter 5 Elastic Strain, Deflection, and Stability 1. Elastic Stress-Strain Relationship

D.A.V. PUBLIC SCHOOL, UPPAL S SOUTHEND, SECTOR 49, GURUGRAM CLASS XI (PHYSICS) Academic plan for

ELASTICITY (MDM 10203)

9 MECHANICAL PROPERTIES OF SOLIDS

Basic Energy Principles in Stiffness Analysis

ROTATING RING. Volume of small element = Rdθbt if weight density of ring = ρ weight of small element = ρrbtdθ. Figure 1 Rotating ring

NCHRP FY 2004 Rotational Limits for Elastomeric Bearings. Final Report. Appendix I. John F. Stanton Charles W. Roeder Peter Mackenzie-Helnwein

Elastoplastic Deformation in a Wedge-Shaped Plate Caused By a Subducting Seamount

Announcements. Suppose your scores are 12, 16, 15. (16/20)*20 + (15/20)*15 + (12/20)*15 = out of a possible

Lecture 7: The Beam Element Equations.

Introduction to Condensed Matter Physics

EMA 3702 Mechanics & Materials Science (Mechanics of Materials) Chapter 2 Stress & Strain - Axial Loading

1 (a) On the axes of Fig. 7.1, sketch a stress against strain graph for a typical ductile material. stress. strain. Fig. 7.1 [2]

INDEX 363. Cartesian coordinates 19,20,42, 67, 83 Cartesian tensors 84, 87, 226

Module 3 : Equilibrium of rods and plates Lecture 15 : Torsion of rods. The Lecture Contains: Torsion of Rods. Torsional Energy

Chapter 2: Basic Governing Equations

NORMAL STRESS. The simplest form of stress is normal stress/direct stress, which is the stress perpendicular to the surface on which it acts.

Transcription:

Physics of Continuous media Sourendu Gupta TIFR, Mumbai, India Classical Mechanics 2012 October 26, 2012

Deformations of continuous media If a body is deformed, we say that the point which originally had coordinates x has coordinates x are the deformation. The deformation itself is given by u = x x. The strain is a measure of the deformation of the body. In order to make a quantitative definition, we begin by measuring lengths of line elements in the body. The length of a piece of the material may change: dl 2 = dx 2 = dx i dx i +2 u i x k dx i dx k + u i x j u i x k dx j dx k The factor 2 in the term in the middle comes after interchanging dummy indices. The strain tensor measures changes of lengths in various direction: dl 2 = dl 2 dl 2 = 2u ij dx i dx j, u ij = 1 2 ( ui x j + u j x i + u i x j u i x k ).

Strain The instantaneous deformation of a body defines a tensor field u ij.

Strain The instantaneous deformation of a body defines a tensor field u ij.

Strain The instantaneous deformation of a body defines a tensor field u ij.

Strain The instantaneous deformation of a body defines a tensor field u ij.

Strain The instantaneous deformation of a body defines a tensor field u ij.

Properties of the strain tensor Since strain is a symmetric second rank tensor, it is possible to choose axes so that it is diagonal. With these principal axes one has dl 2 = (δ ij +2u ij )dx i dx j = (1+2u 1 )dx 2 1+(1+2u 2 )dx 2 2+(1+2u 3 )dx 2 3. In other words, the strain corresponds to stretching or compression along certain principal directions. Clearly, the change in the volume due to a strain is given by dv = dv dv = detudv = u 1 u 2 u 3 dv. If the dimensions of the body are large compared with the strain at any point, then one can linearize the expression for the strain: u ij 1 ( ui + u ) j. 2 x j x i The linear expression may break down for thin rods or plates.

Analogy to a metric A metric tensor is a symmetric tensor of rank two which allows us to define lengths, through the formula dl 2 = g ij dx i dx j. In Euclidean space with Cartesian coordinates the metric tensor is g ij = δ ij. A body under strain can be thought of as a curved space with a changed metric tensor, dl 2 = g ijdx i dx j, where g ij = g ij +2u ij. Problem 65: Exploring the analogy In a curved space defining a derivative operator involves introducing parallel transporters called Christoffel symbols. Explore the definition of derivatives in the curved space of the strained body. What is the physical interpretation of the Christoffel symbols?

The stress tensor When a body is strained there are clearly internal forces which try to return it to its original shape: if the external forces causing the strain are removed, then the body returns to its original shape. These forces are short ranged (except when electric fields are produced by the strain). Since the net force on every small volume vanishes, we can write this as an integral formula by introducing the force on a small volume element, f, through 0 = F i = dv f i = σ ij ds i. Here we used Stoke s theorem to convert the volume integral into an integral over the surface bounding that volume. In order to do this, we need the divergence formula f i = j σ ij. This new tensor of rank two is called the stress tensor.

Internal moments The equilibrium of every element of the body also requires that the net moment at each point vanishes. In other words ( ) σ jk σ ik 0 = M ij = dv(x i f j x j f i ) = dv x i x j x k x k = dv ( xi (x i σ jk x j σ ik ) dv σ jk x ) j σ ik dx k x k x k = ds k (x i σ jk x j σ ik ) dv(σ ij σ ik ) This is a surface integral if the tensor is symmetric. More generally, the second integral vanishes if the difference is a divergence of a rank-3 tensor with suitable symmetries. However, since the only physical quantity that we construct from σ ij is the net force, these symmetries can always be imposed to make the stress tensor symmetric.

Equilibrium of bodies If external forces F act on every unit volume of a body, then the condition for the equilibrium of the body is σ ij x j = F i. The diagonal elements of σ ij correspond to hydrostatic pressure. Assume that the force acting along the normal to a surface element ds i is pds i, then the stress tensor must be given by σ ij = pδ ij. If one works in a frame where there are off-diagonal elements of the stress tensor, then these correspond to shear forces acting along the surface of the body. In the most general case one has pressure as well shear acting on deformable bodies.

The energetics of a deformed body If the displacements u i in a deformed body change by an amount δu i, then the work done by the internal forces is δw = dvδu i j σ ij = ds j δu i σ ij dvσ ij j δu i. If there are no surface deformations, then δw = 1 dvσ ij ( j δu i + i δu j ) = 2 dvσ ij u ij. As a result, one may write the internal energy du = TdS +σ ij du ij, the Helmholtz free energy, df = SdT +σ ij du ij, and the Gibbs free energy dg = SdT u ij dσ ij. As a result, in a state of thermodynamic equilibrium one may write σ ij = F, and u ij = G. T T u ij σ ij

Thermodynamics of a deformed body Thermodynamics is applicable if the stresses and strains are variables of state, i.e., on reversing the strain the body returns to an earlier state. This regime of deformations is called elastic. If the strains are large, then the body may permanently deform. This regime of deformations is called plastic. In the elastic regime, σ ij = 0 when u ij = 0, so the expression for the free energy cannot contain a linear term. The quadratic term is expected to be a rotational scalar, so one may write F = F 0 + λ 2 u2 ii +µu 2 ij, where the constants λ and µ are called Lamé coefficients. The state of thermodynamic equilibrium is obtained by minimizing the free energy. When the strains vanish, this corresponds to F = F 0. For finite strains, the free energy must be higher. This implies that λ > 0 and µ > 0 for thermodynamic stability.

Hooke s law Alternatively, one may decompose the strain into a pure shear and a pure compression through the formula u ij = (u ij 13 ) δ iju kk + 1 3 δ iju kk. Then the free energy may be written as F = F 0 + K ( 2 u2 ii +µ u ij 1 ) 2 3 δ iju kk. K is called the bulk modulus or the modulus of compressibility and µ is called the shear modulus. We can now write df = [Ku kk δ ij +2µ(u ij u kk δ ij /3)]du ij. Using this we obtain Hooke s Law σ ij = Ku kk δ ij +2µ (u ij 13 ) u kkδ ij. This is valid only within the quadratic approximation to F.

Hooke s law Alternatively, one may decompose the strain into a pure shear and a pure compression through the formula u ij = (u ij 13 ) δ iju kk + 1 3 δ iju kk. Then the free energy may be written as F = F 0 + K ( 2 u2 ii +µ u ij 1 ) 2 3 δ iju kk. K is called the bulk modulus or the modulus of compressibility and µ is called the shear modulus. We can now write df = [Ku kk δ ij +2µ(u ij u kk δ ij /3)]du ij. Using this we obtain Hooke s Law σ ij = Ku kk δ ij +2µ (u ij 13 ) u kkδ ij. This is valid only within the quadratic approximation to F.

A simple case Solving a problem of elastic deformations means one should be able to write down all the components of the stress and strain tensors using information on external forces and the elastic moduli. Hydrostatic equilibrium Assume that a body is subjected to an uniform external pressure P, so that the tensor fields u ij and σ ij are constant through the body. Since the shear on the body vanishes, u ij = Uδ ij. The strain tensor is given by σ ij = 3KUδ ij. The condition of equilibrium of forces is σ ij = Pδ ij, since the force on every surface within the solid is given by the inward directed force equal to the pressure times the area of the surface. The solution of the problem is U = P/(3K). The strain is known once P and K are given.

A rod under tension (1) Take a rod with cross sectional area A, and choose the z-axis to be aligned in the long direction. Assume that the rod is under constant tension due to applied outward forces PA on each end of the rod. Assume also that the stresses and strains are uniform through the body. Since there are no sideward forces, one has σ ij n j = 0 for every sideward normal. Also, equilibrium at the ends implies σ zz = P. So the strain tensor has the form 0 0 0 σ = 0 0 0. 0 0 P From Hooke s law we find u ii = σ ii /(3K). Substituting this back gives σ ij = 1 ( 3 σ kkδ ij +2µ u ij 1 ) 9K σ kkδ ij. This allows us to solve for u ij.

A rod under tension (2) Solving for u ij gives u ij = ( σkk 9K ) δ ij + 1 (σ ij 13 ) 2µ σ kkδ ij. This formula is more general than the specific case we will apply it to. Sometimes this form is called Hooke s law. In the case at hand we find u zz = P 3 ( 1 3K + 1 µ ), u xx = u yy = P 3 The Young s modulus is defined to be Y = P u zz, hence Y = 9Kµ 3K +µ. Poisson s ratio is defined as σ = u xx u zz = 1 2 3K 2µ 3K +µ. ( 1 3K 1 ). 2µ

Some problems Problem 66: Elevator cables Take an elevator cable of transverse area A and a mass per unit length ρ. Assume that it is of length L, being hung from a support at the top and hanging vertically due to gravity. Neglect transverse forces on the cable. The vertical downward force at any point along the cable is due to the mass of all the material below it. At the point of suspension the upward force is just balanced by the mass of the full cable. Solve for the stresses and strains given that these vary along the length of the cable. Problem 67: A space elevator A suggestion by Arthur Clarke for a space elevator is essentially that one end of an elevator cable is in circular geostationary orbit and the other end reaches down to the surface of the earth. Compute the stresses and strains on this elevator cable.

Keywords and References Keywords strain tensor, principal axes, thin rods, plates, parallel transporters, Stoke s theorem, stress tensor, hydrostatic forces, shear forces, elastic deformation, plastic deformation, Lamé coefficients, thermodynamic stability, pure shear, pure compression, bulk modulus, modulus of compressibility, shear modulus, Hooke s Law, Young s modulus, Poisson s ratio, elevator cable, space elevator References Landau and Lifshitz (volume 7) Theory of Elasticity : sections 1 5.