3.1.
−=+−=
−++−=
=
+=−+=
=
==
−==
−==
=
cossin)(cossincossin
)sin(coscossincossin
cossincossin2cossin
: plane obliqueon Stresses
cos
sin
)cos,sin(,
000
00
00
22
2222
XYYX
YX
S
YnT
XnT
nT
Y
X
xyyxxyS
xyyxyS
yy
n
y
xx
n
x
jij
n
i
ij
n
3.2*.
:Plots MATLAB
2sin2cos12292sin62cos219
sin2cos4
sin6cossin4cos
sincos4sincos2
)0,sin,(cos,
120
261
014
(b))0,sin,(cos,
104
041
412
(a)
222222 −−=++=++=++=
=++=−=++=
−=++=+=++=
+=++=+=++=
=
−==
−
−
=
n
z
n
y
n
x
nn
z
n
y
n
x
n
zzyzyxzx
n
zzzyzyxzx
n
z
zyzyyxyx
n
yzyzyyxyx
n
y
zxzyxyxx
n
xzxzyxyxx
n
x
ijij
TTTTTT
nnnTnnnT
nnnTnnnT
nnnTnnnT
TT
nn
x
y
S
n
S
S
6
6.5
Angle,
part(a)
part(b)
3.3.
+
−
=
−++−=
=
−
−
−
+
=
−+=
=
+
−
+
+
=
++=
=
2cos2sin
2
)sin(coscossincossin
2sin2cos
22
cossin2cossin
2sin2cos
22
cossin2sincos
22
22
22
xy
xy
xyyxxyr
xy
yxyx
xyyxy
xy
yxyx
xyyxxr
3.4.
)sin(coscossincossin
cossin2cossin
cossin
sincos
cossin
sincos
cossin
sincos
:MatrixRotation D–2
22
22
22
−+−=
−+=
−
−
==
−
=
rrxy
rry
rrx
r
rr
T
QQ
Q
x
y
x’
y’
r
r
r
3.5.
)
22
)
22
)(4)()[(
2
1
0)(
00)det(
000
0
0
stress Plane
2
2
21
max
2
2
22
2,1
22
xy
yx
xy
yxyx
xyyxyxyx
xyyxyx
yxy
xyx
ijij
yxy
xyx
ij
+
−
=
−
=
+
−
+
=
−−++=
=−++−
=
−
−
=−
=
3.6.
2
13
2
32
2
21
2
2
122
2
1
2
321
2
3
2
2
2
1
2222
1321
)()()(
9
1
62
9
1
2)2(2
9
1
9
)(
9
1
)(
3
1
)(||
3
1
3
1
)(
3
1
3
1
,
3
1
,
3
1
with scoordinate principal Using
jiijkikjij
n
oct
kkjiijoct
IIIII
nnnnN
Inn
−+−+−=
−=−−=
++−++=
−=−=
==++==
=
T
n
3.7.
( )
2
2
2
2
max
2
2
2
2
2,1
2
2
1
21
max
2
2
112
2
11
22
2,1
22
2222
so system, coordinate theoft independen bemust thusand
invariants theof functions are stressesshear max and principal the
4
1
2
Equation) sticCharacteri D–2 start with could(or
4
1
2
1
4
4
1
2
1
)(4)()[(
2
1
+
−
=+
−
=
+
−
+
=+
−
+
=
−=
−
=
−=−=
−−++=
r
r
xy
yx
r
rr
xy
yxyx
xyyxyxyx
II
IIIIII
3.8.
)(
0
)(0
)(2
4
3
:Directions Principal
)(
4
3
22
0
2
4
3
)2,1(
2
22
22222
)2,1(
1
)2,1(
2
)2,1(
1
2222222
2222222
3
22222
3
2
2
2,1
22
22
3
−
+−
=
=
+−−−
−+−−
+−=+
−
+
=
−
−
=
=
n
yc
yxycxy
n
n
n
yxycxyyc
ycyxycxyxy
c
P
yxycxy
c
P
yc
ycxy
c
P
xy
yxyx
yxy
xyx
ij
3.9*.
Using MATLAB the first principal stress contours are:
3.10.
( )
( ) ( ) ( )
P
ngxfnPgxPfnxPf
P
n
P
xfnxPfn
xPf
jijkijjijkijjijkij
kk
jijkijjijkijjijij
kijij
on dependnot willdirections principal the thusand
0)(0)(0)(
from follow directions principal The
form theof be wouldproblem thisof valuesprincipal the thusand
0)(0)(0)(
)( with problem valueprincipal theup Setting
)3,2,1(
3,2,1
)3,2,1(
3,2,1
)3,2,1(
3,2,1
3,2,13,2,13,2,1
3,2,1
=−=−=−
=
−=−=−
=
3.11*.
: Plotand Code MATLAB
)()(
)2
)(
)(
)
2
)(
2
,
)(
2
,
)(
2
22
2
222
2
2
222
22
2
2
max
222
2
222
3
222
2
yx
Py
yx
Pxy
yx
yxPy
yx
Pxy
yx
Py
yx
yPx
xy
yx
xyyx
+
=
+
−
+
+
−−
=+
−
=
+
−=
+
−=
+
−=
3.12.
2
13
2
1
2
3
2
32
2
3
2
2
2
21
2
2
2
1
22
33
2
22
2
11
2
3
2
3
2
2
2
2
2
1
2
1
2
2
2
2
33
2
22
2
11
2
3
2
3
2
2
2
2
2
1
2
1
2n
)()()(
)(
|T|
−+−+−=
++−++=−=
++=
++=
nnnnnn
nnnnnnNS
nnnN
nnn
n
T
% Exercise 3-6
% Maximum Shear Stress Contours for Flamant Problem
clear;
% Generation of variable space
[x,y]=meshgrid(-3:0.05:3,0:0.05:3);
% Calculation of Nondimensional Maximum Shear Stress
sxy=y./(x.^2+y.^2);
hold on
axis equal
axis off
% Plot contours with reversed y-axis
contour(x,-y,sxy,20);
3.13.
2
1
))((
)(
2
1
)(
2
1
)(
4
1
))((
))((
0
))((
0
))((
))((
2
1
))((
)(
2
1
)(
2
1
)(
4
1
))((
))((
)9.4.3( Using
)(
2
1
0))(()(
4
1
0))(( :circle ofequation in the)(
2
1
Using
2313
231131
2
31
2313
21
2
2
3
12321232
13
2
2
2
3121
331231
2
31
3121
2
1
3113
2
31
13
2
31max
=
−−
−+
−++−
=
−−
−−+
=
=
−−
=
−−
−−+
=
=
−−
−+
−++−
=
−−
−−+
=
+==−−+−
=−−+−=
NNS
n
NNS
n
NNS
n
NNN
NNSS
3.14.
−
=
−
−
−
−==
−
−=
=
200
010
004
2/12/10
3/13/13/1
6/16/16/2
021
201
113
2/12/10
3/13/13/1
6/16/16/2
on theoryansformati tensor trstandard Using
matrixn nsformatioproper tra thebe would
2/12/10
3/13/13/1
6/16/16/2
,directions principal determined previously theusing andmatrix stressgiven theis
021
201
113
T
mnjnimij
ij
ij
QQ
Q
3.15.
.
ˆ
stress deviatoric theofdirection
principal thealso Therefore .0)
ˆ
(or 0)
3
1
( as written be
can thisand ,0)( then, tensor stress theofdirection principal theis If
.
3
1
or
3
1
that followsit
thusand ,0]det[ is tensor stress for theequation sticcharacteri theHowever,
.0])
3
1
(det[]
3
1
det[
as written becan equation this,
3
1
ˆ
Since
.0]
ˆ
det[by given is tensor stress deviatoric for theequation sticcharacteri The
ij
jijdijjijdkkij
jijijij
kkddkk
ijijij
ijdkkijijdijkkij
ijkkijij
ijdij
nn
n
=−=
+−
=−
−=+=
=−
=+−=−−
−=
=−
n
n
3.16.
−=−=
−−
−−
=−=
−
−
−=
−
−
=
3/420
23/171
013/13
3
1
ˆ
3/404
03/51
413/1
3
1
ˆ
3/100
003/1
3
1
~
3/700
003/7
3
1
~
120
261
014
(b)
104
041
412
(a)
ijkkijijijkkijij
ijkkijijkkij
ijij
3.17.
449.2654
3
1
])6()3()3[(
3
1
])()()[(
3
1
1)(
3
196.527])6()3()3[(
2
1
])()()[(
2
1
2,1,4
021
201
113
With
2/12222/12
13
2
32
2
21
321
2/12222/12
2
2
321
===++=−+−+−=
=++=
==++=−+−+−=
−===
=
oct
oct
ij
3.18.
3/6
3
2
0
33
)(
3
13
2
)(
3
1
)(
3
1
3/1 plane, octahedral For the
)0,1,1(2/10)(
:Case
)1,0,0(0)(
:Case0
)0,1,1(2/10)(
:Case
,0,
0)(
00
0
0
]det[
000
00
00
2
2
321
2
33
2
22
2
11
2222
3
2
2
2
1
2
3
2
3
2
2
2
2
2
1
2
1
2
321
)3()3(
3
3
)2()2(
2
2
)1()1(
1
1
321
22
==−=
=
−
=++=++=
=+=++=++=
===
−==−
−=
==−
=
==−
=
−===
=−=
−
−
−
=−
=
NS
nnnN
nnn
nnn
n
n
n
n
n
jijij
jijij
jijij
ijij
ij
T
T
n
n
n
3.19*.
0,)(
3
0
2
4
3
7,–3 Exercise From
22222
2
2
22
22
3
=+−=+
−
+
=
−
−
=
=
yxyx
yxy
xyx
ij
yxycxy
P
yc
ycxy
c
P
von Mises Stress Contours
3.20.
2
22 2
1
2
1
2
tan1
tan2
2tan
2
2
02sin2cos
22
−
+
−
−=
−
=
−
=
−
−
=
+
−
+
+
=
xy
yx
xy
yx
p
p
p
yx
xy
pxy
yx
xy
yxyx
x
dx
dy
dx
dy
dx
dy
d
d
d
d
3.21.
yxxy
xy
xy
xy
yx
yx
yx
y
yxy
yxy
xy
xy
y
y
y
y
x
yx
x
xyx
yx
yx
x
x
x
x
dydxdxdy
dx
dy
dx
dydx
x
dy
dx
dy
dxdy
y
centerM
F
yx
dydxdxdy
dxdyFdydydx
x
dxdxdy
y
F
F
yx
dydxdxdy
dxdyFdxdxdy
y
dydydx
x
F
=
→
=−
+−
+
+
=
=+
+
→
=+−
++
−
+
=
=+
+
→
=+−
++
−
+
=
0,let and by through Divide
0
22
)(
22
)(
)(@0
0
0,let and by through Divide
0)(
)(
0
0
0,let and by through Divide
0)(
)(
0
x
y
dx
x
x
x
+
xy
yx
dy
y
y
y
+
dx
x
xy
xy
+
dy
y
yx
yx
+
x
F
y
F
3.22.
satisfiedy identicall are equations mequilibriuboth
00
2
22
2
2
2
−
+
=
+
x
yyxxyx
yxy
y
xyx
yxy
yx
x
3.23.
rr
r
r
rr
r
r
rr
r
r
rr
r
r
r
rr
r
rr
r
rr
r
r
r
ddrrdrd
rd
rd
dr
dr
r
centerM
F
rrr
ddrrdrd
rdrdFddrddrd
ddrddrdrdddrrdr
r
F
F
rrr
ddrrdrd
rdrdFddrddrdddr
ddrdrdddrrdr
r
F
=
→
=
=+
+
+
→
=+−
+
+
++−+
+
=
=+
−
+
+
→
=+−
+−−
++−+
+
=
0,let and by through Divide
22
22
)(@0
02
1
0,let and by through Divide
0)2/cos()2/cos()(
)2/sin()2/sin()())((
0
0
1
0,let and by through Divide
0)2/sin()2/sin()()2/cos(
)2/cos()())((
0
r
dr
r
r
r
+
r
r
+
d
dr
r
r
r
+
+
d
r
r
r
F
F
dr
rd
d
3.24.
equations. mequilibriu esatisfy thnot do
stresses materials of mechanics eapproximat The
000
3
)(
0
3
2
0
0,
3
)(
,
22
22
==+
+
+
−
−−=+
+
+
====
−
=−=
z
z
yz
xz
zyyxy
x
zx
yx
x
yzxzzyxyx
F
zyx
dx
dV
I
yR
zyx
V
I
y
dx
dM
I
y
F
zyx
I
yRV
I
My
3.25.
000
000
)(00
0,,0,)(
==+
+
+
==+
+
+
−=−==+
=+
+
+
=========
z
z
yz
xz
y
zyyxy
x
l
x
x
x
zx
yx
x
zyxxzyzxyzyxx
F
zyx
F
zyx
xlggdxg
x
F
zyx
FFgFx
3.26.
F
p
F
F
y
p
F
zyx
x
p
zyx
FporFpFpF
p
p
p
p
zz
z
yz
xz
yy
zyyxy
xx
zx
yx
x
iiiijjijij
ijij
=
=+
+
+
=
=+
+
+
===+−=+
−
−
−
=−=
0
0
approachscalar a Using
00
00
00
00
,,,
3.27.
( )
gives forcebody edappropriat the
adding and direction coordinate each in termsCollecting
1
)(
1
)(
1
)(
1
)(
11
1
)(
1
16,–1 exerciseor (1.9.14) From
zrz
rz
r
r
+
+
+
+−
++
+
+++
+
+
=
++
+++
+
+++++
=
+
++
=
+
+
=
++=
++=
++=
++=
z
z
r
r
rzrrz
r
z
r
zzz
r
rrrr
σττ
z
τστ
r
ττσ
r
ττσ
r
zrrr
z
u
u
r
ru
rr
σττ
τστ
ττσ
zθr
zrθθr
zθrθr
zzθθzrrzzθzθθrrθ
zrzθrθrrzrzθrθrr
z
θ
r
zzθθzrrzz
zθzθθrrθθ
zrzθrθrrr
zzθθrr
eee
eeeee
eeeeee
eeeeee
eeeeee
T
T
T
T
σ
u
eeeT
eeeT
eeeT
TeTeTe