How to Relate Stress and Strain

How to Relate Stress and Strain
I
We have considered deformation of continua and some of
the balance laws in them.
I
These kinematics and mechanics apply to all the
continuous media. Then, where do the characteristics of
individual material come from?
I
Properties unique to a certain material are determined by
the material’s internal constitution or physical make-up.
The quantitative expressions for such internal constitution
are called constitutive equations / laws / relations /
models.
I
The role of constitutive relations is to relate stress (a
quantity needed for stating the linear momentum balance
law) and strain (a quantity used to describe motion of
continua).
Rules of Constitutive Theory
I
In this lecture, we briefly review rules for formulating a
constitutive theory and then introduce linear elasticity.
I
The constitutive relations are those between the internal
deformation quantified by strain and interfal forces
quantified by stress.
I
The constitutive relations for a given material must give
meaningful response to all classes of deformations and
temperatures to which the body may be subjected.
I
In case we need to develop a totally new constitutive
theory for a newly synthesized material for instance, we
need some guidelines by which a consitutive model can be
formulated.
Rules of Constitutive Theory
I
Principle of Equipresence:
An independent variable assumed to be present in one
constitutive equation of a material should be assumed to
be present in all constitutive equations of the same
material unless its presence contradicts an assumed
symmetry of the material or contradicts of material
frame-indifference or some other fundamental principle.
I
As one can see, this “principle” comes with many
conditions. Indeed, the development of a constitutive theory
depends highly on experience: researchers tend to avoid
complicating their equations with additional variables
unless experience (e.g., through experiments) indicates
that the variable are needed.
I
Still, when applied, this principle can help revelaing the
coupling effects between different kinds of phenomena that
might otherwise be overlooked.
Rules of Constitutive Theory
Three fundamental postulates are assumed to be valid for any
constitutive theory of purely mechanical phenomena in a
continuous medium.
1. Principle of determinism for stress:
The stress in a body is determined by the history of the
motion of that body.
I
This principle includes as special cases the Newtonian
fluid, which is history-independent, and linear elasticity, for
which there is only one hisotry, the reference state.
Viscoelasticity and plasticity are examples of more
complicated history-dependent theories.
Rules of Constitutive Theory
2. Principle of local action:
In determining the stress at a given particle, the motion
outside an arbitrary neighborhood of the particle may be
disregarded.
I
Another seemingly obvious principle. However, there are
cases where we need to break this principle and consider
the responses from a finite neighborhood around a particle.
Such a theory is called “non-local”.
Rules of Constitutive Theory
3. Principle of material frame-indifference:
Constitutive equations must be invariant under changes of
frame of references. That is, two observers, even if in
relative motion (in the sense of the Galilean relativity) with
respect to each other, observe the same stress in a given
body.
I
The assumed situation is that one observer is moving
(translation + rotation) relative to another. Distinguish the
chage of frame of references from coordinate
transformation, the latter of which is time-independent and
does not need to assume multiple observers.
Rules of Constitutive Theory
More on the material frame indifference
I
Events (x,t):
x∗ = c(t) + Q(t)x
t∗ = t − a
I
I
I
(1)
Vectors (v):
v∗ = Q(t)v
(2)
T∗ = Q(t)TQ(t)T
(3)
Rank 2 tensors (T):
Deformation gradient tensor (F):
F∗ = Q(t)F
(4)
Rules of Constitutive Theory
I
Frame-indifference: Functions and fields whose values
are scalars, vectors, or tensors are called
frame-indifferent or objective if both the dependent and
the independent vector and tensor variables transform
according to Eq. (2) to (4), while the scalar variable are
unchanged (e.g., readings of a thermometer should be the
same for all the observers however they are moving).
I
What does it mean for constitutive relations?
I
If a constitutive equation is satisfied for a dynamical
process with a motion and a symmetric stress tensor, then
the constitutive equation must also be satisfied by the
motion and stress tensor that are transformed as in Eq. (1)
and (3),
Rules of Constitutive Theory
I
The tedious and other-worldly discussion on seemingly
obvious statements sometimes pose practical challenges.
I
Sometimes, the kinematics of a continuum is formulated in
terms of velocity rather than displacements. This point of
view makes a lot of sense when a body is deforming in a
time-dependent fashion. For instance, lithosphere that is
extending at a certain rate.
I
In such a case, the deformation is better represented by
strain rate, not strain. If we assume a linear constitutive
relation between stress and strain, as in linear elasticity
we’ll see next, the constitutive relation will return stress
rate, not stress:
σ = Cε ⇒ σ˙ = Cε˙
Rules of Constitutive Theory
I
Now the question is whether the stress rate is material
frame-indifferent. Predictably, the answer is no.
I
To show this, we start with
σ ∗ = QσQT
(5)
By taking material time derivative, we get
T
˙
˙T
σ˙ ∗ = QσQ
˙ T + QσQ
+ Q(t)σ Q
Because of the last two terms in Eq. (6), σ˙ is not
frame-indifferent.
(6)
Rules of Constitutive Theory
I
A slightly different definition of stress rate can be shown to
be frame-indifferent but there are more than one of such
stress rate forms.
I
The most popular (due to easy computation) choice is
called Jaumann stress rate and defined as follows
σ
˚ = σ˙ + σW − Wσ
(7)
where W is the spin tensor, the skew-symmetric part of the
velocity gradient (dv/dx).
I
Showing σ
˚ is frame-indifferent (also called objective) is
straightforward. How to get an objective stress rate is a
separate question and to answer it takes a whole new set
of mathematical machinery. So, we skip that part in this
course.
Linear Elasticity
I
Hooke’s law for a 1 dimensional mass-spring system:
F = −kx
I
If no damping force acts on it, the system is conservative,
meaning by definition that there is a potential function U(x)
such that F = −∇U.
I
In this 1D example, integration to get U is straightforward
and U = 12 kx 2 .
Linear Elasticity
I
A material is called ideally elastic when a body formed of
the material recovers its original form completely upon
removal of the forces causing the deformation, and there is
a one-to-one relationship between the state of stress and
the state of strain, for a given temperature.
I
The one-to-one relationship precludes behaviors like creep
at constant load or stress relaxation at constant strain.
I
The classical elastic constitutive equations, often called the
generalized Hooke’s law, are nine equations expressing
the stress components as linear homonenous (i.e., all the
terms are of the same power) functions of the nine strain
components:
σij = Cijkl εkl
(8)
Linear Elasticity
I
The rank 4 tensor, Cijkl , has 81(= 34 ) components.
I
However, recall that stress and strain tensor are
symmetric: i.e., σij = σji and εkl = εlk .
I
Thus,
Cijkl = Cjikl and Cijkl = Cijlk .
I
We further consider the case in which the material is
elastically isotropic, i.e., there are no preferred directions
in the material. Then, the elastic constants (Cijkl ) must be
the same at a given particle for all possible choices of
rectangular Cartesian coordinates in which stress and
strain components are evaluated.
Linear Elasticity
I
The most general rank 4 tensor that satisfy all of the above
symmetry and isotropy conditions is
Cijkl = λδij δkl + µ(δik δjl + δil δjk )
(9)
(see Malvern Sec. 6.1 and 6.2 for further details.)
I
I
The constitutive relation becomes
σij = λδij δkl + µ(δik δjl + δil δjk ) εkl
(10)
Finally, after some simplification, we reach the isotropic
generalized Hooke’s law:
σij = λεkk δij + 2µεij ,
where λ and µ are called Lamé’s constants.
(11)
Linear Elasticity
I
By setting i = j in Eq. (11), we find
σii = (3λ + 2µ)εii
I
By substituting εkk = σkk /(3λ + 2µ) (from Eq. (12)) into Eq.
(11), we obtain
εij = −
I
(12)
1
λ
σkk δij +
σij
2µ(3λ + 2µ)
2µ
(13)
Recall the definitions of Young’s modulus and Poisson’s
relation:
I
I
Hooke’s law: σxx = Eεxx , etc, where E is the Young’s
modulus.
Poisson’s relation: e.g., εyy = εzz = −νεxx = − Eν σxx , where
ν is the Poisson’s ratio.
Linear Elasticity
I
From these relations, we derive the following for 3D:
1
[σxx − ν(σyy + σzz )]
E
1
= [σyy − ν(σzz + σxx )]
E
1
= [σzz − ν(σxx + σyy )]
E
1
=
σxy , etc.
2µ
εxx =
εyy
εzz
εxy
I
(14)
The above set of equations can be generalized to the
following indicial notation:
ν
1+ν
εij = − σkk δij +
σij
E
E
(15)
Linear Elasticity
I
It is also convenient to decompose the stress and strain
tensor into deviatoric and volumetric parts:
1
sij = σij − σkk δij
3
1
ij = εij − εkk δij
3
I
(16)
Then, the whole relationship can be expressed by the two
equations:
sij = 2µij and p = −Ke,
(17)
where p = −σkk /3 is the pressure, e = εkk si the volume
strain, and K is the bulk modulus, related to Lamé’s
constants by the relation
2
K =λ+ µ
3
(18)
Linear Elasticity
I
In the case of isotropic elasticity, all the elastic constants
(λ, µ, E, ν, K ) can be expressed in terms of just two
independent constants.
I
µ, the shear modulus, is often denoted as G.
I
Poisson’s ratio: Which is better for wine bottle, cork or
rubber?
Negative Poisson’s ratio?
http://www.youtube.com/watch?v=HJ1Ck6FIqwU
Linear Elasticity: Plane Stress
I
Stress state in which σ3j = 0 with j = 1 . . . 3.
I
In terms of principal stress, σ1 6= 0, σ2 6= 0 and σ3 = 0.
I
From Eq. (15),
1
(σ1 − νσ2 )
E
1
ε2 = (σ2 − νσ1 )
E
ν
ε3 = − (σ1 + σ2 )
E
ε1 =
I
Let’s assume that the lithosphere is in a lithostatic
condition: i.e., diviatroic stresses are zero and
σ1 = σ2 = σ3 = ρgh.
(19)
Linear Elasticity: Plane Stress
I
Furthermore, deviatoric stresses are in the plane stress
condition and the two horizontal principal stresses are
equal in magnitude.
∆σ1 = ∆σ2 6= 0 but ∆σ3 = 0.
I
Then, we get
ε1 = ε2 =
ε3 = −
I
1−ν
∆σ1
E
2ν
∆σ1 .
E
(20)
Let’s think about the arguments given by Turcotte and
Schubert for the plane stress condition in the lithosphere.