2.3 - Stress Strain Relations
2.3 - Stress Strain Relations
2.3 - Stress Strain Relations
CHAPTER
for Linearly Elastic
Solids 3
3.1 INTRODUCTION
In the preceding two chapters we dealt with the state of stress at a point and the
state of strain at a point. The strain components were related to the displacement
components through six of Cauchy’s strain-displacement relationships. In this
chapter, the relationships between the stress and strain components will be estab-
lished. Such equations are termed constitutive equations. They depend on the
manner in which the material resists deformation.
The constitutive equations are mathematical descriptions of the physical phe-
nomena based on experimental observations and established principles.
Consequently, they are approximations of the true behavioural pattern, since an
accurate mathematical representation of the physical phenomena would be too
complicated and unworkable.
The constitutive equations describe the behaviour of a material, not the
behaviour of a body. Therefore, the equations relate the state of stress at a point
to the state of strain at the point.
On the same lines, if e1, e2 and e3 are the principal strains, which are also along x,
y, z directions, the normal strain in x¢ direction, from Eq. (2.20), is
2 2 2
ex¢x¢ = nx1 ε1 + n y1 ε 2 + nz1 ε 3 (3.7a)
The shear strain gx¢y¢ is obtained from Eq. (2.36c) as
gx¢y¢ =
1
(
⎡ 2 nx1nx 2ε1 + n y1n y 2ε 2 + nz1nz 2ε 3 )
(1 + ε x′ ) (1 + ε y′ ) ⎣
+ nx1nx 2 + n y1n y 2 + nz1nz 2 ⎤⎦
Using Eq. (3.5), and observing that ex¢ and ey¢ are small compared to unity in the
denominator,
gx¢y¢ = 2 (nx1nx 2ε1 + n y1n y 2ε 2 + nz1nz 2ε 3 ) (3.7b)
Substituting the values of s 1 , s 2 and s 3 from Eqs (3.4a)–(3.4c) into
Eq. (3.6b)
tx¢y¢ = nx1nx 2 (λ ∆ + 2 µ ε1 ) + n y1n y 2 (λ ∆ + 2 µ ε 2 ) + nz1nz 2 (λ ∆ + 2 µ ε 3 )
2
( 2 2
)
2
( 2 2
= λ ∆ nx1 + n y1 + nz1 + 2µ ε1nx1 + ε 2 n y1 + ε 3 nz1 )
= l D + 2m ex¢x¢ (3.11a)
Similarly, one gets
sy¢ = l D + 2m ey¢y¢ (3.11b)
sz¢ = l D + 2m ez¢z¢ (3.11c)
Similar to Eq. (3.8),
ty¢z¢ = m gy¢z¢ (3.12a)
tx¢z¢ = m gz¢x¢ (3.12b)
Equations (3.11a)–(3.11c), (3.8) and (3.12a) and (3.12b) relate the six rectangular
stress components to six rectangular strain components and in these only two
elastic constants are involved. Therefore, the Hooke’s law for an isotropic mate-
rial will involve two independent elastic constants l and m (or G).
e 1 = 1 [σ 1 − v (σ 2 + σ 3 )]
E
where E is Young’s modulus, and n is Poisson’s ratio. Comparing this with
Eq. (3.15),
µ ( 3λ + 2 µ ) λ
E= ; ν= (3.16)
(λ + µ ) 2 (λ + µ )
K=
(3 λ
+ 2µ )
, G ∫ m, l=
νE , (3.18)
3 (1 + ν ) (1 − 2ν )
It should be observed from Eq. (3.17) that for the bulk modulus to be positive,
the value of Poisson’s ratio n cannot exceed 1/2. This is the upper limit for n.
For n = 1/2,
3G = E and K = •
Stress–Strain Relations for Linearly Elastic Solids 103
A material having Poisson’s ratio equal to 1/2 is known as an incompressible
material, since the volumetric strain for such an isotropic material is zero.
For easy reference one can collect the equations relating stresses and strains
that have been obtained so far.
(i) In terms of principal stresses and principal strains:
s 1 = l D + 2me1
s 2 = l D + 2me2 (3.19)
s 3 = l D + 2me3
where D = e1 + e2 + e3 = J1.
λ+µ ⎡ λ ⎤
e1 = σ − (σ + σ 3 ) ⎥
µ ( 3λ + 2µ ) ⎢⎣ 1 2 ( λ + µ ) 2 ⎦
λ+µ ⎡ λ ⎤
e2 = ⎢σ 2 − (σ 3 + σ 1 ) ⎥ (3.20)
µ ( 3λ + 2µ ) ⎣ 2 ( λ + µ ) ⎦
λ+µ ⎡ λ ⎤
e3 = ⎢σ 3 − (σ 1 + σ 2 ) ⎥
µ ( 3λ + 2µ )
⎣ 2 ( λ + µ ) ⎦
(ii) In terms of rectangular stress and strain components referred to an orthogo-
nal coordinate system Oxyz:
sx = l D + 2mexx
sy = l D + 2meyy
s z = l D + 2mezz (3.21a)
where D = exx + eyy + ezz = J1.
txy = mgxy, tyz = mgyz, tzx = mgzx (3.21b)
λ+µ ⎡ λ ⎤
exx =
µ ( 3λ + 2µ ) ⎢σ x −
2 (λ + µ )
(
σy +σz ⎥ )
⎣ ⎦
λ+µ ⎡ λ ⎤
eyy = ⎢σ y − (σ z + σ x ) ⎥ (3.22a)
µ ( 3λ + 2µ ) ⎣ 2 (λ + µ ) ⎦
λ+µ ⎡ λ ⎤
ezz = ⎢σ z −
µ ( 3λ + 2µ ) ⎣ 2 (λ + µ )
( σ x + σ y )⎥
⎦
sx = E ⎡ ν ⎤
∆ + ε xx ⎥
⎢
(1 + ν ) ⎣(1 − 2ν ) ⎦
= l J1 + 2Gexx
104 Advanced Mechanics of Solids
sy = E ⎡ ν ⎤
∆ + ε yy ⎥ (3.23a)
⎢
(1 + ν ) ⎣ (1 − 2ν ) ⎦
= l J1 + 2Geyy
sz = E ⎡ ν ⎤
∆ + ε zz ⎥
⎢
(1 + ν ) ⎣(1 − 2ν ) ⎦
= l J1 + 2Gezz
txy = Ggxy, tyz = Ggyz, txx = Ggzx (3.23b)
(
exx = 1 ⎡σ x − ν σ y + σ z ⎤
E⎣ ⎦ )
eyy = 1 ⎡⎣σ y − ν (σ z + σ x ) ⎤⎦ (3.24a)
E
(
ezz = 1 ⎡σ z − ν σ x + σ y ⎤
E ⎣ ⎦ )
gxy = 1 txy, gyz = 1 tyz, gzx = 1 tzx (3.24b)
G G G
∂σ x ∂τ xy ∂τ zx
+ + =0
∂x ∂y ∂z
For an isotropic material
sx = λ ∆ + 2 µ ε xx ; τ xy = µ γ xy ; τ xz = µ γ xz
Hence, the above equation becomes
⎛ ∂ε ∂γ xy ∂γ xz ⎞
λ ∂∆ + µ ⎜ 2 xx + + =0
∂x ⎝ ∂x ∂y ∂ z ⎟⎠
From Cauchy’s strain-displacement relations
∂ ux ∂ ux ∂ u y ∂ ux ∂ uz
exx = , gxy = + , gzx = +
∂x ∂y ∂x ∂z ∂x
Stress–Strain Relations for Linearly Elastic Solids 105
Substituting these
⎛ ∂ 2u
λ ∂∆ + µ ⎜⎜ 2x +
∂ 2ux ∂ 2ux ⎞ ∂ ⎛ ∂ ux + ∂ u y + ∂ uz ⎞
+ ⎟⎟ + µ
∂ x ⎜⎝ ∂ x ∂ z ⎟⎠
or =0
∂x ⎝ ∂x ∂y 2 2
∂z ⎠ ∂y
Observing that
∂ ux ∂ u y ∂ uz
D = ε xx + ε yy + ε zz = + +
∂x ∂y ∂z
⎛∂u ∂ u y ∂ uz ⎞ ⎛ ∂ 2 u x ∂ 2 u x ∂ 2u x ⎞
( λ + µ ) ∂∂x ⎜ ∂ xx + + ⎟ + µ ⎜⎜ + + ⎟=0
⎝ ∂y ∂z ⎠ ⎝ ∂x
2
∂ y2 ∂ z 2 ⎟⎠
This is one of the displacement equations of equilibrium. Using the notation
∂2 + ∂2 + ∂2
—2 =
∂ x2 ∂ y 2 ∂ z 2
the displacement equation of equilibrium becomes
( λ + µ ) ∂∂∆x + µ ∇ 2 u x =0 (3.25a)
Similarly, from the second and third equations of equilibrium, one gets
( λ + µ ) ∂∂∆y + µ ∇2 u y =0 (3.25b)
( λ + µ ) ∂∂∆z + µ ∇ 2 u z =0
Example 3.1 A rubber cube is inserted in a cavity of the same form and size
in a steel block and the top of the cube is pressed by a steel block with a
pressure of p pascals. Considering the steel to be absolutely hard and
assuming that there is no friction between steel and rubber, find (i) the
pressure of rubber against the box walls, and (ii) the extremum shear stresses in
rubber.
106 Advanced Mechanics of Solids
l
y
x
Fig. 3.1 Example 3.1
Solution
(i) Let l be the dimension of the cube. Since the cube is constrained in x and
y directions
exx = 0 and eyy = 0
and sz = –p
Therefore
( )
exx = 1 ⎡σ x − ν σ y + σ z ⎤ = 0
E⎣ ⎦
eyy = 1 ⎡⎣σ y − ν (σ x + σ z ) ⎤⎦ = 0
E
Solving
ν s =– ν
sx = sy = p
1− ν z 1 −ν
If Poisson’s ratio = 0.5, then
sx = sy = sz = –p
(ii) The extremum shear stresses are
σ1 − σ 3 σ − σ2 σ − σ3
t2 = , τ3 = 1 , τ1 = 2
2 2 2
If n £0.5, then sx and sy are numerically less than or equal to sz. Since sx,
sy and sz are all compressive
s1 = sx = – ν p
1 −ν
s2 = sy = – ν p
1 −ν
s 3 = sz = –p
\ (
t1 = p 1 − ν =
1 −ν 1 −ν)
1 − 2ν
p, τ 2 =
1 − 2ν
1 −ν
p, τ 3 = 0
Solution Since the shear stresses on x, y and z planes are zero, the given stresses are
principal stresses. Arranging such that s 1 > _ s 2 _> s 3
s 1 = 100 MPa, s2 = –20 MPa, s3 = – 40 MPa
The extremal shear stresses are
E 2 × 105
G= = = 8 ¥ 104 MPa
2 (1 + ν ) 2 × 1.25
The principal shear strains are therefore
τ1 10
g1 = = = 1.25 ¥ 10–4
G 8 × 104
τ2 70 = − 8.75 × 10−4
g2 = =−
G 8 × 104
τ3 60
g3 = = = 7.5 ¥ 10–4
G 8 × 104
From Eq. (1.44a), the octahedral shear stress is
⎡ Ans. (a) 120 ¥ 106 kPa (1.22 ¥ 106 kgf/cm2), 80 ¥ 106 kPa ⎡
⎢ (8.1680 ¥ 105 kgf/cm2) ⎢
⎢ (b) 7.8 ¥10 kPa (7.96 ¥ 10 kgf/cm2), 11.7 ¥ 106 kPa ⎢
6 4
⎢ ⎢
⎣ (1.2 ¥ 105 kgf/cm2) ⎣
3.2 For steel, the following data is applicable:
E = 207 ¥ 106 kPa (2.1 ¥ 106 kgf/cm2),
and G = 80 ¥ 106 kPa (0.82 ¥ 106 kgf/cm2)
For the given strain matrix at a point, determine the stress matrix.
⎡ 0.001 0 −0.002 ⎤
[eij] = ⎢ 0 −0.003 0.0003⎥⎥
⎢
⎢⎣ −0.002 0.003 0 ⎥⎦
⎡ ⎡ −68.4 0 −160⎤ ⎡
⎢ ⎢ 0 −708.4 24 ⎥ ⎢
⎢ Ans. [tij] = ⎢ ⎥ ¥ 10 kPa ⎢
3
⎢ ⎣⎢ −160 24 −228.4⎦⎥ ⎢
⎣ ⎣
3.3 A thin rubber sheet is enclosed between two fixed hard steel plates (see
Fig. 3.2). Friction between the rubber and steel faces is negligible. If
the rubber plate is subjected to stresses sx and sy as shown, determine the
strains exx and eyy, and also the stress ezz
z ⎤ Ans. sz = +n (sx + sy) ⎤
⎥ ⎥
⎥ exx = + 1 + ν [(1 – n)sx – nsy] ⎥
sy
y ⎥ E ⎥
⎥ ⎥
⎥ e = + 1 + ν [(1 – n)s ns ] ⎥
x
sx ⎥⎦ yy E y x ⎥⎦
Fig. 3.2 Example 3.2