16.1.
functions ionapproximator ioninterpolat theare)(
2
1
),( where
),(
)(
2
1
)(
2
1
)(
2
1
])()()[(
2
1
),(
),(for form original theinto results thesengsubstitutiBack
,, and element, theofarea theis wherewhere
)(
2
1
,)(
2
1
,)(
2
1
constants for the equations three theSolving
),(
),(
),(
node eachat Evaluating
),( :ionApproximatLinear
3
1
333322221111
332211332211332211
332211333221123322111
33321333
23221222
13121111
321
yx
A
yx
yxu
uyx
A
uyx
A
uyx
A
yuuuxuuuuuu
A
yxu
yxu
xxyyyxyxA
uuu
A
cuuu
A
cuuu
A
c
ycxccuyxu
ycxccuyxu
ycxccuyxu
ycxccyxu
iii
e
i
iii
eee
e
jkikjijkkjie
eee
++==
=
++++++++=
++++++++=
−=−=−=
++=++=++=
++==
++==
++==
++=
=
16.2.
=
=
=
=
=
663663662662661661
223123222122221121
123113122112121111
33
33
22
22
11
11
332211
321
321
66
2212
1211
33
33
22
22
11
11
66
2212
1211
332211
321
321
0
0
0
0
0
0
4
000
000
00
0
0
0
0
0
0
0
0
4
][
00
0
0
][,000
000
2
1
][
with,]][[][][
element,r triangulastrainconstant For the
CCCCCC
CCCCCC
CCCCCC
A
h
C
CC
CC
A
h
K
C
CC
CC
C
A
B
BCBAhK
e
e
e
e
e
T
ee
16.3.
T
yxyxyx
ee
exxx
y
x
y
x
y
x
y
x
T
y
x
T
y
x
T
e
FFFFFF
Ah
F
AFdxdyFdxdyF
F
F
F
F
F
F
F
F
F
F
dxdy
F
F
hF
ee
e
3
}{
componentsother allfor likewise and,
3
1
forces,body constant wise–element of case For the
000
000
][
][}{by given is vector forcebody thegeneral, In
11
3
3
2
2
1
1
321
321
=
==
=
=
=
16.4.
12
12
1233
21
12
11111
12
12
3
3
2
2
1
1
0
0
2
][
sevaluation integral theof all Collecting
. on 0 since vanish will terms thecontaining integrals while
patterns, identicalfollow and involving componentsother theof Integrals
2
)(
2
1
Now
0
0
2
][
ractionsboundary tconstant withelement,angular linear tri For the
12
121212
1212
=
=
=++==
=
=
n
y
n
x
n
y
n
x
e
n
y
n
x
T
e
e
e
n
xe
n
xe
n
xe
n
y
n
x
n
y
n
x
e
n
y
n
x
n
y
n
x
n
y
n
x
e
n
y
n
x
T
e
T
T
T
T
Lh
ds
T
T
h
Lh
dsyx
A
ThdsThdsTh
T
T
T
T
Lh
ds
T
T
T
T
T
T
hds
T
T
h
16.5.
−
−
−
−
−
−
−−
−
−−
−−
−
−
−
−
−
=
=+=+
+
=+
====−=−===
−
−
−−
−
−
+
−
−
−
−
−
−
−−
−−
−
=
−
−=+=+
−
−=+
===−====−=
+
==
−
=
−
==
2
32
1
2
32
1
2
1
2
1101
101 2
1
2
1
2
1
00
2
1
)1(2
][
matrix stiffness thein componentsother allfor similarly and
...,0,0,
)1(2
1,2/1,1,0,1,1,1,0 , 2element For
1
0
2
1
1
2
1
2
32
1
2
1
2
3
0
2
1
2
1
2
1
2
10101
)1(2
][
matrix stiffness thein componentsother allfor similarly and
...,
1
,0,
1
1,2/1,1,1,0,0,1,1 1,–15 Examplein 1element For
)1(2
,
1
,
1
:case stress plane isotropic For the
2
)2(
662111216611121166
2
111
2
1
1321321
2
)1(
2
6621112166111211
2
66
2
111
2
1
1321321
66
2
12
2
2211
E
CCCC
E
CC
hA
E
E
CCCC
E
CC
hA
E
C
E
C
E
CC
K
K
16.6.
mT
V
U
V
U
T
T
V
U
V
U
T
T
V
U
V
U
E
T
T
V
U
V
U
KK
KKKK
KKK
KKKK
11
3
3
2
2
3
3
2
2
9
3
3
2
2
3
3
2
2
)2(
44
)1(
66
)2(
34
)1(
56
)2(
33
)1(
55
)1(
46
)1(
45
)1(
44
)1(
36
)1(
35
)1(
34
)1(
33
10
030.0
441.0
081.0
492.0
system theSolving
0
2/
0
2/
8.151
08.151
4.1104.418.151
6.274.41698.151
10
0
2/
0
2/
2
3
0
2
3
1
2
1
2
32
1
2
1
2
3
)1(2
0
2/
0
2/
bygiven matrix was system global reduced the1,–15 Example From
−
−
=
=
−
−−
=
−
−
−
−−
−
−
+
−
−
−
=
+
++
16.7.*
MATLAB PDE Toolbox: FEA Solution to Biaxial Loading Problem in Exercise 8-14
1000
1000
x
y
Comparison of x Results:
FEA Solution: (x)A 200 , (x)B 1900
Analytical Solution: (x)A 0 , (x)B 2000
•
•
A
B
Comparison of y Results:
FEA Solution: (y)A 1900 , (y)B 300
Analytical Solution: (y)A 2000 , (y)B 0
•
•
A
B
x – Contours
y – Contours
16.8.*
MATLAB PDE Toolbox: FEA Solution to Biaxial-Shear Loading Problem in
Exercise 8-15 (Figure 8-15)
1000
1000
x
y
x – Contours
Comparison of x Results:
FEA Solution: (x)A -270 , (x)B 3300
Analytical Solution: (x)A 0 , (x)B 4000
•
•
A
B
Comparison of y Results:
FEA Solution: (y)A -3500 , (y)B 280
Analytical Solution: (y)A -4000 , (y)B 0
•
•
A
B
y – Contours
16.9.*
MATLAB PDE Toolbox: FEA Solution to Curved Beam Problem shown in Figure 8-32.
Choose case b/a = 4, to provide simple match with analytical results given in Figure 8-33
1000
Analytical Solution: (x)A -10,000 , (x)B 2571
16.10.*
MATLAB PDE Toolbox: FEA Solution to Torsion Problem in Exercise 9-17/18
A
FEA Model:
(b = 0.267, a = 0.8, 2129 Nodes, 4096 Elements)
•
B
FEA | | – Contours
75.0,2.1 : ResultsFEA; 2 : ResultsAnalytical
)( :stressshear Resultant 2
2
2
22
==
+
=+=
BA
BA
yzxz yx
16.11.
r
rrrrrr
r
rrrrrr
r
rrrrrr
r
GGGT
rrrrrr
r
rr
r
r
G
r
x
r
rr
r
G
ikj
kjiikjjikkji
kjiikjjkiijk
ljllljjllljl
ikjkijiklljikj
kjiikjjkiijk
kjiij
k
kij
i
i
jiijij
−−
=
−++−−
−
+
−++−−
−
+
−++−−
−
=
++=
−++−−
−
=
+−−
−
=
=
+
−
−
=
,
,,,,,,
,,,,,,
,,,,,,
,,,
,,,,,,
,,,
,
,
,
,,
)21(2
2
2)43(
)1(8
1
2)43(
)1(8
1
2)43(
)1(8
)(
by given is form stress The
2)43(
)1(8
1
)()43(
)1(8
1
fact that theUsing
1
ln)43(
)1(8
1
:case strain plane l,dimensiona– twoFor the