Finite Element Method
Finite Element Method
Fig. 1
The nodes of the spring element are indicated by points i and
j in Fig. 1.
Let
k1 = AE
L
where
[1]
f1 = k1 u1 u2 = k1u1 k1u2
f 2 = k1 u2
u1 = k1u1 + k1u2
[2]
[3]
f
2
k1
k
1
k1 u1
u
k1 2
local stiffness
matrix
[4]
the
e
e
K1 = K1
k1
k
1
k1
k1
f ji ,
k1 u1
u
k1 2
f 2 k
1 2
2=
f k2
2
k2 u2
k2 3
Fig. 2
Assembly Method 1
We expand the two matrix equations so that they are in
equivalent form
f 1 k
1 1
1
f = k1
2
0 0
K1e =
k1
k
1
0
k1 0 u1
k1 0 u2
0 0 3
k1 0
k1 0
0 0
0 0
2
f1 = 0
f 2 0
2
K2e =
0
0 u1
k2 k2 u2
k2 k2 u3
0 0
k2 k2
k2 k2
f11
1
2
f + f
1
2
2
k1
k
1
0
k1
0
k1 + k2 k2
k2 k2
u
1
u
2
u
3
[5]
global stiffness
matrix
K = K1e +K2e
[6]
Note:
1) It is imperative that the element stiffness matrices are first
expanded before proceeding with the addition operation.
2) The global stiffness matrix will have dimensions n n ,
where n is the number of degrees of freedom (DOF) of
the system.
Assembly Method 2
An alternative way of assembling the stiffness matrix for the
system is to consider the equilibrium of forces for each node of
the system.
f11 = F1
F1 = k1u1 k1u2
f 21 + f12 = F2
F2 = k1u1 + ( k1 + k2 ) u2 k2u3
f 22 = F3
F3 = k2u2 + k2u3
Assembly Method 3
A 3rd method for performing the assembly of elementary
stiffness matrices is to use a table of connectivity.
Node 2
3
global node
numbers
u1 = 0 (node 1 is fixed)
F2 = F3 = P
We have
f11
1
2
f + f
1
2
2
F
1
F2
F
3
k1
k
1
0
k1
0
k1 + k2 k2
k2 k2
u
1
u
2
u
3
Therefore
F
1
P
P
k1
k
1
0
k1
0
k1 + k2 k2
k2 k2
0
u
2
u
3
which reduces to
P
P
k1 + k2
k
2
k2 u2
k2 3
and
F1 = k1.u2
u
= u2
3
2P
u
k1
2
=
P
u3 2P
+
k2
k1
F1 = 2P
Problem 1
Problem 2
Truss
russ element having 2 DOF / NODE
Note that the analysis is based on the fact that the elements
used sustain only axial forces (bars or rods). The terms bar,
rod or truss member
member can be used indifferently in what follows.
For this particular configuration,
configuration we consider the
e following:
a) each node of a bar has 2 degrees of freedom
b) the bar is no longer aligned with the global x coordinate
axis
Fig. 3
5.
Discretization
Formulation: model applied to individual elements (in
local coordinates)
Choice of a common reference system (global
coordinates)
Use of coordinate transformation to express the model
for individual elements in terms of global coordinates.
(globalization)
Assembly of elements.
y)
1 , u1 , u1
u11x , u1y
2x
2y
(nodal displacements)
1 , f1 , f1
f11x , f1y
2x
2y
(nodal forces)
k1
u
f
K1e
(stiffness constant)
(nodal disp. vector)
(nodal force vector)
(element stiffness matrix)
Global coordinates: ( X , Y )
(nodal displacements)
(nodal forces)
U
F
K1G
Coordinate Transformations
u1
1x
1
u
1y =
1
u2 x
u1
2 y
cos = c , sin = s
the
s
c
0
0
0
0
c
s
local
0
0
and
U 1
1X
1
U
1Y
U 1
2X
1
U 2Y
global
nodal
[7]
f1
1x
0
1 =
f 2x
s
c
0
0
0
0
c
s
0
0
F1
1X
1
F
1Y
F1
2X
1
F2Y
[8]
0 =
1
f 2x
k
1
0
k
1
0
0
0
0
k1 0
0 0
k1 0
0 0
u1
1x
1
u
1y
1
u2 x
u1
2 y
[9]
s
c
0
0
0
0
c
s
0
0
F1
1X
1
F
1Y
F1
2X
1
F2Y
k
1
0
k
1
0
0
0
0
k1
0
k1
0
0
0
0
0
s
c
0
0
0
0
c
s
0
0
U 1
1X
1
U
1Y
U 1
2X
1
U 2Y
T .F = K1e .T . U
[10]
T T .T .F = T T .K1e .T . U
T T .T = I
F = T T .K1e .T . U
Stiffness matrix
written in global
coordinates : K1G
[11]
K1G =
s
c
0
0
0
0
c
s
0 k1
0 0
s k1
c 0
0
0
0
0
k1
0
k1
0
0
0
0
0
s
c
0
0
0
0
c
s
0
0
1
0
1
0
0
0
s
c
0
0
0
0
c
s
0
0
k1 s
0
s
c
0
0
0
0
c
s
0
0
0
0
0
0
K1G =
k1 s
0
s
c
0
0
0
0
c
s
0
0
0
c
s
0
s
0
K1G =
c2
k1 cs2
c
cs
cs
s2
cs
s2
c2
cs
c2
cs
cs
s 2
cs
s2
c
0
c
0
s
0
s
0
Example
Let us consider the case of a three member truss.
P.
Problem
Solution :
Problem
A1 = A2 = A3 = 80mm2
E = 200GPa
Solution :