Module 3: Analysis of Strain: Ohr'S Ircle For Train
Module 3: Analysis of Strain: Ohr'S Ircle For Train
g
below the e-line; and the y-axis point a distance
2
g
above the e-line; and vice versa when the shear strain is negative.
2
By analogy with stress, the principal strain directions are found from the equations
tan 2q =
g xy
ex -ey
(3.19)
ex +ey
ex -ey
e1,2 =
2
2
g xy
+
2
2
(3.20)
2e x
y 2
2g xy
xy
2
3u e y
3v
=
,
=
x 2 y
xy 2 x 2
3u
3v
+
xy 2
x 2 y
(3.21)
1
Applied Elasticity for Engineers
Module3/Lesson2
or
2e y
2g xy
2e x
+
=
y 2
xy
x 2
This is the condition of compatibility of the two dimensional problem, expressed in terms of
strain. The three-dimensional equations of compatibility are derived in a similar manner:
Thus, in order to ensure a single-valued, continuous solution for the displacement
components, certain restrictions have to be imposed on the strain components.
These resulting equations are termed the compatibility equations.
Suppose if we consider a triangle ABC before straining a body [Figure 3.4(a)] then the same
triangle may take up one of the two possible positions Figure 3.4(b) and Figure 3.4(c)) after
straining, if an arbitrary strain field is specified. A gap or an overlapping may occur, unless
the specified strain field obeys the necessary compatibility conditions.
Now,
u
x
v
ey =
y
w
ez =
z
v u
g xy =
+
x y
w v
g yz =
+
y z
ex =
(3.23)
(3.23a)
(3.23b)
(3.23c)
(3.23d)
2
Applied Elasticity for Engineers
Module3/Lesson2
g zx =
u w
+
z x
(3.23e)
Differentiating Equation (3.23) with respect to y and Equation (3.23a) with respect to x
twice, we get
2e x
3u
=
y 2
xy 2
2e y
3v
=
x 2
yx 2
(3.23f)
(3.23g)
(3.23h)
Taking the derivative of Equation (3.23c) with respect to x and y together, we get
2g xy
3v
3u
=
+
xy yx 2 xy 2
(3.23i)
(3.23j)
2e y
2
2e z g yz
=
y 2
yz
(3.23k)
2e z 2e x 2g zx
+
=
x 2
z 2
xz
(3.23l)
z 2
Now, take the mixed derivative of Equation (3.23) with respect to z and y,
2e x
3u
=
yz xyz
(3.23m)
And taking the partial derivative of Equation (3.23c) with respect to z and x, we get
2 g xy
xz
3u
3v
+
xyz zx 2
(3.23n)
Also taking the partial derivative of Equation (3.23d) with respect to x twice, we get
2 g yz
x 2
3w
3v
+
x 2 y x 2 z
(3.23p)
3
Applied Elasticity for Engineers
Module3/Lesson2
2 g zx
3u
3w
Thus,
=
+
xy xyz x 2 y
(3.23q)
Now, adding Equations (3.23n) and (3.23q) and subtracting Equation (3.23p), we get
2 g yz
-
x 2
2 g xz 2 g xy
2 3 u
+
+
=
xy
xz xyz
(3.23r)
2 2 e x
g yz g xz g xy
= +
+
yz x x
y
z
(3.23s)
2 2 e y
g zx g yx g yz
+
+
xz
y y
z
x
2 2 e z
g xy g yz g zx
= +
+
xy z z
x
y
=
(3.23t)
(3.23u)
Thus the following are the six compatibility equations for a three dimensional system.
2
2
2 e x e y g xy
+
=
xy
y 2
x 2
2e y
z 2
2
2 e z g yz
+
=
yz
y 2
2 e z 2e x 2g zx
+
=
zx
x 2
z 2
2
2 e x
g yz g xz g xy
= +
+
yz
x x
y
z
2 2 e y
(3.24)
g yz g zx g xy
+
y x
y
z
zx
2
2 e z
g yz g zx g xy
=
+
xy
z x
y
z
=
4
Applied Elasticity for Engineers
Module3/Lesson2
g
ex +ey ex -ey
+
cos 2q + xy sin 2q
e q =
2
2 2
g q = (e y - e x )sin 2q + g xy cos 2q
(3.25)
(3.26)
e +ey
e max or e min = x
2
1
2
(e
- e y ) + g xy2
2
(3.27)
g xy
tan 2q =
e -e
y
x
(3.28)
Also, the maximum shear strain at the point is given by following relation.
g max =
(e
- e y ) + g xy2
2
(3.29)
5
Applied Elasticity for Engineers
Module3/Lesson2
g
e + ey ex -ey
+
cos 2q1 + xy sin 2q1
e q1 = x
2
2 2
g
e + ey ex -ey
+
cos 2q 2 + xy sin 2q 2
e q 2 = x
2
2 2
g
e + ey ex -ey
+
cos 2q 3 + xy sin 2q 3
e q3 = x
2
2 2
(a)
(b)
(c)
ex +ey ex -ey
+
+ 0
e 0 =
2 2
1
= (e x + e y + e x - e y )
2
\e 0 = e x
g
ex +ey ex -ey
+
(0) + xy
e 45 =
2
2 2
1
= (e x + e y + g xy )
2
or 2e 45 = e x + e y + g xy
\ g xy = 2e 45 - (e x + e y )
ex + ey
2
g
ex -ey
+
(cos 1800 ) + xy (sin 1800 )
2
2
ex +ey ex -ey
-
=
2 2
1
= (e x + e y - e x + e y )
2
=ey
Also, e 90 =
\ e 90
e x = e 0 , e y = e 90 0
and g xy = 2e 45 - (e 0 + e 90 )
6
Applied Elasticity for Engineers
Module3/Lesson2
q 1 = 0, q 2 = 60 0 , q 3 = 120 0
Substituting the above values in (a), (b) and (c), we get
e x = e0, e y =
and g xy =
2
3
1
(2e 60 + 2e 120 - e 0 )
3
(e 60 - e 120 )
Hence, using the values of e x , e y and g xy , the principal strains e max and e min can be
computed.
ex = -20010-6
ey = 100010-6
gxy = 90010-6
(a) Find the strain components associated with a set of axes x y inclined at an angle of
30o clockwise to the x y set as shown in the Figure 3.5. Also find the principal
strains and the direction of the axes on which they act.
y
y
300
300
x
Figure 3.5
7
Applied Elasticity for Engineers
Module3/Lesson2
Solution: (a)
The transformation equations for strains similar to that for stresses can be written as below:
e x =
ex +ey ex -ey
g xy
+
cos2q +
sin2q
2
2
2
e y =
g xy
ex +ey ex -e y
cos2q sin2q
2
2
2
g x 'y '
ex -ey
= -
2
2
sin2q + xy cos2q
2
450
0
2q = tan-1
= 36.8
600
Therefore,
g x 'y '
= - (750 10 -6 )sin (60 0 - 36.8 0 )
2
= -295 10 -6
hence, g xy = - 590 10 -6
Solution: (b)
From the Mohrs circle of strain, the Principal strains are
e 1 = 1150 10 -6
e 2 = -350 10 -6
8
Applied Elasticity for Engineers
Module3/Lesson2
9
Applied Elasticity for Engineers
Module3/Lesson2
The directions of the principal axes of strain are shown in figure below.
e1
e2
71.6o
Figure 3.7
Example 3.2
By means of strain rosette, the following strains were recorded during the test on a
structural member.
e x = e0
e y = e 90
g xy = 2e 45 - (e 0 + e 90 )
e x = -13 10 -6
e y = 13 10 -6
g xy = 2 7.5 10 -6 - (- 12 10 -6 + 13 10 -6 )
\ g xy = 15 10 -6
ex +ey 1
(e x - e y )2 + g xy2
e max or e min =
2 2
- 13 + 13 -6 1
[
\ e max or e min =
(- 13 - 13)10 -6 ]2 + (15 10 -6 )2
10
2
2
-6
\ e max or e min = 15 10
10
Applied Elasticity for Engineers
Module3/Lesson2
tan 2q =
(e
g xy
-ey)
15 10 -6
(- 13 - 13)10 -6
tan 2q = -0.577
\ 2q = 150 0
\q = 75 0
=
e 45 = -110 micrometres/m
e 90 = 210 micrometres/m
Find the magnitudes and directions of principal strains.
Solution: Given e 0 = 750 10 -6
e 45 = -110 10 -6
e 90 = 210 10 -6
Now, for a rectangular rosette,
e x = e 0 = 750 10 -6
e y = e 90 = 210 10 -6
g xy = 2e 45 - (e 0 + e 90 )
) (
g xy = -1180 10 -6
\ The magnitudes of principal strains are
ex +ey
e max or e min =
2
1
2
(e
- e y ) + g xy2
2
11
Applied Elasticity for Engineers
Module3/Lesson2
i.e.,
750 + 210 -6 1
e max or e min =
10
2
2
= 480 10 -6
-6 2
1
(1297.7 )10 -6
2
= 480 10 -6 648.85 10 -6
\ e max = e 1 = 1128.85 10 -6
e min = e 2 = -168.85 10 -6
The directions of the principal strains are given by the relation
tan 2q =
(e
\ tan 2q =
g xy
x
-ey)
- 1180 10 -6
= -2.185
(750 - 210)10 -6
\ 2q = 114.6 0
ex =
u
= 2 x 10 -3 ,
x
ey =
v
= 6 yz 10 -3 ,
y
ez =
w
= 3 10 -3
z
g xy = (x 2 + 3) 10 -3 + (3 y 2 z ) 10 -3
x
y
g xy = 0
12
Applied Elasticity for Engineers
Module3/Lesson2
g yz = (3 y 2 z 10 -3 ) + ( x + 3 z ) 10 -3
y
z
g yz = 3 y 2 10 -3
and g zx = ( x + 3 z )10 -3 +
x 2 + 3 10 -3
z
x
g zx = 1 10 -3
Therefore at point (1, 2, 3), we get
e x = 2 10 -3 , e y = 6 2 3 10 -3 = 36 10 -3 , e z = 3 10 -3 ,
g xy = 0, g yz = 12 10 -3 , g zx = 1 10 -3
Example 3.5
The strain components at a point with respect to x y z co-ordinate system are
x
x
y
1
y
z
Here l1 =
2
1
1
2
0
1
2
2
1
1
l2 =
, m2 =
,
2
2
l 3 = 0, m3 = 0, n3 = 1
m1 = -
z(Z)
n1 = 0
45
n2 = 0
x
x
Now, we have,
Figure 3.8
13
Applied Elasticity for Engineers
Module3/Lesson2
[e ] = [a ][e ][a ]T
1
2
[a ][e ] = - 1
2
0.127
= - 0.014
0.08
0
0.1 0.08 0.08
0 0.08 0.2 0.08
0
.
08
0
.
08
0
.
3
0.198 0.113
0.085
0
0.08
0.3
1
2
1
2
0
0.08
0.3
0.08
1
2
1
2
0
1
2
1
2
0
g yz = 0, g zx = 0.113 2 = 0.226
Example 3.6
The components of strain at a point in a body are as follows:
14
Applied Elasticity for Engineers
Module3/Lesson2
ex
g
e ij = xy
2
g xz
g xy
2
g xz
2 0 .1
0.15 - 0.04
g yz
= 0.15 - 0.05 0.05
2
0
.
04
0
.
05
0
.
05
ez
ey
g yz
2
1 2
g xy + g 2 yz + g 2 zx
4
\ J 2 = -0.0291
J 3 = (0.1)(- 0.05)(0.05) +
J 3 = -0.002145
1
(0.3)2 + (0.1)2 + (- 0.08)2
4
1
(0.3)(0.1) (-0.08) - 0.1(0.1)2 + 0.05 (0.08)2 - 0.05 (0.3)2
4
(i)
3
1
cos q - cos 3q = 0
4
4
Let e = r cos q +
= r cos q +
(ii)
J1
3
0 .1
3
e = r cos q + 0.033
\ (i) can be written as
15
Applied Elasticity for Engineers
Module3/Lesson2
0.03251
0.00112
cos q =0
2
r
r3
(iii)
0.03251 3
=
r2
4
4 0.03251
i.e., r =
= 0.2082
3
cos 3q 0.00112
and
=
4
r3
4 0.00112
or cos 3q =
= 0.496 @ 0.5
(0.2082)3
\ 3q = 60 0
60
= 20 0
3
q 3 = 140 0
or q1 =
q 2 = 100 0
\e 1 = r1 cosq1 +
J1
3
= 0.2082 cos 20 0 +
e 1 = 0.228
0 .1
3
J1
0 .1
= 0.2082 cos 100 0 +
= -0.0031
3
3
J
0 .1
e 3 = r3 cos q 3 + 1 = 0.2082 cos 1400 +
= -0.126
3
3
e 2 = r2 cos q 2 +
16
Applied Elasticity for Engineers
Module3/Lesson2
0.15
- 0.04
(0.1 - e 1 )
0.15
(- 0.05 - e 1 )
0.05
- 0.04
(0.05 - e 1 )
0.05
0.15
- 0.04
(0.1 - 0.228)
=
0.15
(- 0.05 - 0.228)
0.05
- 0.04
0.05
(0.05 - 0.228)
0.15
- 0.04
- 0.128
= 0.15
- 0.278
0.05
- 0.04
0.05
- 0.178
Now, A1 =
- 0.278
0.05
0.05
- 0.178
\ A1 = 0.046984
B1 = -
0.15
0.05
- 0.04 - 0.178
\ B1 = 0.0247
C1 =
0.15
- 0.278
- 0.04
0.05
\C1 = -0.00362
A1 + B1 + C1 =
2
= 0.0532
\ l1 =
m1 =
A1
A1 + B1 + C1
2
B1
A1 + B1 + C1
2
0.046984
= 0.883
0.0532
0.0247
= 0.464
0.0532
17
Applied Elasticity for Engineers
Module3/Lesson2
C1
n1 =
A1 + B1 + C1
2
- 0.00362
= -0.068
0.0532
(0.1 + 0.0031)
0.15
- 0.04
0.15
(- 0.05 + 0.0031)
0.05
- 0.04
0.05
(0.05 + 0.0031)
0.15
- 0.04
0.1031
A2 + B2 + C 2 =
\ l2 =
m2 =
n2 =
A2
A22 + B22 + C 22
B2
A22 + B22 + C 22
C2
A + B +C
2
2
2
2
2
2
- 0.00499
= -0.399
0.0125
- 0.009965
= -0.797
0.0125
0.00562
= 0.450
0.0125
= 0.0125
(0.1 + 0.126)
0.15
- 0.04
0.15
- 0.04
0.05
0.05
(0.05 + 0.126)
(- 0.05 + 0.126)
0.15
- 0.04
= 0.15 0.076
- 0.04 0.05
0.05
0.176
0.226
18
Applied Elasticity for Engineers
Module3/Lesson2
Now, A3 =
0.076
0.05
C3 = Now,
\ l3 =
m3 =
n3 =
0.15
0.076
- 0.04
0.05
A3 + B3 + C 3 =
2
A3
A + B +C
2
3
2
3
2
3
B3
A32 + B32 + C32
C3
A32 + B32 + C32
0.0109
= 0.338
0.0322
- 0.0284
= -0.882
0.0322
0.01054
= 0.327
0.0322
= 0.0322
Example 3.7
The displacement components in a strained body are as follows:
19
Applied Elasticity for Engineers
Module3/Lesson2
1
1
g xy
g xz
ex
2
2
1
1
e ij = g yx
ey
g yz
2
2
1
1 g
g zy
ez
zx
2
2
1
3
(e n )oct
= e x l 2 + e y m 2 + e z n 2 + 2(e xy lm + e yz mn + e zx nl )
1
1
e yz = g yz and e zx = g zx
2
2
2
2
2
1
1
1
\ (e n )oct = 0.0001
+ 0.0003
+ 0.0005
+
3
3
3
1
1
1
2 0.0002 + 0.0004 + 0.0005
3
3
3
Here e xy =
1
g xy ,
2
\ (e n )oct = 0.001
Octahedral Shearing Strain is given by
20
Applied Elasticity for Engineers
Module3/Lesson2
\ (e R )oct =
1
(e x + e xy + e xz )2 + (e xy + e y + e yz )2 + (e xz + e yz + e y )2
3
1
(0.0001 + 0.0002 + 0.0005)2 + (0.0002 + 0.0003 + 0.0004)2 (0.0005 + 0.0004 + 0.0005)2
3
\ (e R )oct = 0.001066
\ g oct = 0.000739
(b) Deviator and Spherical strain tensors.
Here Mean Strain = e m =
ex +ey +ez
3
0.0001 + 0.0003 + 0.0005
=
3
\ e m = 0.0003
0.0002
(0.0001 - 0.0003)
(0.0003 - 0.0003)
\ Deviator Strain tensor =
0.0002
0.0005
0.0004
- 0.0002 0.0002 0.0005
i.e.,
E = 0.0002
0
0.0004
0.0005 0.0004 0.0002
0
e m 0
0
0 e m
i.e.,
0.0004
(0.0005 - 0.0003)
0.0005
0
0
0.0003
E = 0
0.0003
0
0
0
0.0003
Example 3.9
The components of strain at a point in a body are as follows:
e x = c1 z (x 2 + y 2 )
e y = x2 z
g xy = 2c 2 xyz
where c1 and c2 are constants. Check whether the strain field is compatible one?
21
Applied Elasticity for Engineers
Module3/Lesson2
Solution: For the compatibility condition of the strain field, the system of strains must
satisfy the compatibility equations
2
2
2 e x e y g xy
i.e.,
+
=
xy
y 2
x 2
e x
2e x
= 2c1 yz ,
= 2c1 z
y
y 2
e y
2e y
= 2 xz ,
= 2z
x
x 2
g xy
2g xy
= 2c2 yz ,
= 2c2 z
x
xy
2
2g xy
2e x e y
\ 2 +
= 2c1 z + 2 z = 2 z (1 + c1 ) and
= 2c 2 z
xy
y
x 2
2
2
2 e x e y g xy
Since
+
e x = 2axy 2 + by 2 + 2cxy
e y = ax 2 + bx
g xy = ax 2 y + bxy + ax 2 + hy
Solution: For compatibility, the strain components must satisfy the compatibility equation.
2
2
2e x e y g xy
i.e.,
+
=
y 2
x 2
xy
or
(i)
2
2
2e x e y g xy
+
=0
y 2
x 2
xy
(ii)
e x
= 4axy + 2by + 2cx
y
2e x
= 4ax + 2b
y 2
22
Applied Elasticity for Engineers
Module3/Lesson2
e y = ax 2 + bx
e y
x
2e y
x 2
= 2ax + b
= 2a
g xy = ax 2 y + bxy + ax 2 + hy
g xy
x
2g xy
xy
= 2axy + by + 2ax
= 2ax + b
\ (i) becomes
4ax + 2b + 2a = 2ax + b
(
)
v = c(4 x + 2 y + z )
u = c x 2 + 2x
w = 4cz 2
where c is a very small constant, determine the strain at (2,1,3), in the direction
1 1
0, ,
2 2
u
v u
= 2cx, g xy =
+
= 4c + 0 = 4c
x
x y
v
w v
ey =
= 4cy, g yz =
+
= 0+c = c
y
y z
w
u w
ez =
= 8cz , g zx =
+
= 2c + 0 = 2c
z
z x
\ At point (2,1,3),
Solution: e x =
23
Applied Elasticity for Engineers
Module3/Lesson2
e x = 2c 2 = 4c, g xy = 4c
e y = 4c 1 = 4c, g yz = c
e z = 8c 3 = 24c, g zx = 2c
\ The Resultant strain in the direction l = 0, m = -
1
1
, n=
is given by
2
2
e r = e x l 2 + e y m 2 + e z n 2 + g xy lm + g yz mn + g zx nl
2
1
1
1 1
= 0 + 4c + 24c
+ 4c(0) + c
+ 2c(0)
2
2 2
\ e r = 13.5 c
Example 3.12
The strain components at a point are given by
1 1 1
,
,
3 3 3
(e n )oct
\ (e n )oct
= e x l 2 + e y m 2 + e z n 2 + g xy lm + g yz mn + g zx nl
1
= [0.01 - 0.02 + 0.03 + 0.015 + 0.02 - 0.01]
3
= 0.015
ex
e xy
e
xz
e xy
ey
e yz
0.015
0.01
0.01
2
2 0.01
e xz
0.0075 - 0.005
0
.
015
0
.
02
= 0.0075 - 0.02
e yz =
- 0.02
0.01
2
2
e z 0.01 0.02
0.03
- 0.005 0.01
0.03
2
2
(e R )oct
1
(e x + e xy + e xz )2 + (e xy + e y + e yz )2 + (e xz + e yz + e z )2
3
1
(0.01 + 0.0075 - 0.005)2 + (0.0075 - 0.02 + 0.01)2 + (- 0.005 + 0.01 + 0.03)2
3
= 0.0004625
=
24
Applied Elasticity for Engineers
Module3/Lesson2
\ (e R )oct = 0.0215
and octahedral shearing strain is given by
(e S )oct = 2 (e R )2 - (e n )2
\ (e S )oct = 0.031
= 2 (0.0215) - (0.015)
2
Example 3.13
The displacement field is given by
u = K x 2 + 2 z , v = K 4 x + 2 y 2 + z , w = 4 Kz 2
where K is a very small constant. What are the strains at (2,2,3) in directions
1
1
(a) l = 0, m =
, n=
, (b) l = 1, m = n = 0, (c) l = 0.6, m = 0, n = 0.8
2
2
u
v
= 2 Kx, e y =
= 4 Ky,
x
y
v u
g xy = +
= 4K + 0 = 4K
x y
w v
g yz =
+
= 0+ K = K
y z
u w
g zx =
+
= 2K + 0 = 2K
z x
\ At point (2,2,3),
e x = 4 K , e y = 8 K , e z = 24 K
g xy = 4 K , g yz = K , g zx = 2 K
Solution: e x =
ez =
w
= 8 Kz
z
e r = e x l 2 + e y m 2 + e z n 2 + g xy lm + g yz mn + g zx nl
(i)
Case (a)
Substituting the values in expression (i), we get
2
1
1
1 1
e r = 4 K (0) + 8K
+ 24 K
+ 4 K (0) + K
+ 2 K (0)
2
2
2 2
1
\ e r = 4 K + 12 K + K
2
\ e r = 16.5 K
Case (b)
e r = 4 K (1) + 8 K (0 ) + 24(0 ) + 4 K (0 ) + K (0 ) + 2 K (0 )
\ e r = 4K
2
25
Applied Elasticity for Engineers
Module3/Lesson2
Case (c)
26
Applied Elasticity for Engineers