Finite Element Methods Reg 2024 Week 4
Finite Element Methods Reg 2024 Week 4
(MN234504 / 3 SKS)
Example:
E = 210 Gpa
A = 5.0 x104 m2
Element 1
2
DEVELOPMENT OF TRUSS EQUATIONS
Element 2
Element 3
Boundary Conditions
3
DEVELOPMENT OF TRUSS EQUATIONS
Displacement
4
DEVELOPMENT OF TRUSS EQUATIONS
Example:
The truss is composed of four nodes and three elements. A =
10x104 m2 and E=210 GPa for all the elements. A load of 20 kN
is applied at node 1 in the global x-direction. Nodes 2–4 are pin
supported and constrained from movement in the x, y, and z.
5
DEVELOPMENT OF TRUSS EQUATIONS
6
DEVELOPMENT OF TRUSS EQUATIONS
Boundary Conditions
7
DEVELOPMENT OF TRUSS EQUATIONS
Displacements
Element Stress
The Transformation :
9
Where :
The Transformation :
→ → →
Therefore :
Example :
Element 1 :
10
Element 2 :
Element 3 :
11
The Global Stiffness Matrix :
12
The Next Step :
Then
and
13
The Global Nodal Forces :
where :
Therefore :
14
Exercises:
For the plane trusses with inclined supports shown in Figures
a, b and c, solve for the nodal displacements and element
stresses in the bars. Let A = 2 in2, E = 30 106 psi, and L= 30
in for each truss.
a b c
15
SYMMETRY : correspondence in
a. Size
b. Shape
c. Position of Loads
d. Material Properties
e. Boundary Conditions
Axial Stiffness :
1, 2, 7, 8
3, 4, 5, 6
Element s 4 and 5
Half Areas → 22
The Boundary Conditions
The Force-Displacement Equation
or
Therefore and
24
On Applying the Boundary Conditions
and
25
The Governing Equation
Solving
or
26
Using Partitioning Matrix → Rearranging the Equation
Solving
27
SYMMETRY AND BANDWIDTH
Exercises :
For the truss shown in figures below, use symmetry to
determine the displacements of the nodes and the stresses
in each element.
a. All elements have E = 30 106 psi. Elements 1, 2, 4, and
5 have A =10 in2 and element 3 has A = 20 in2.
b. All elements of the structure in have the same AE
except element 1, which has an axial stiffness of 2AE.
a b
28
29
30
Banded Symmetric Matrices and Bandwidth
The coefficient matrix (stiffness matrix) for the linear equations
in structural analysis is always symmetric and banded
The skyline
An Envelope begins with the first
nonzero coefficient in each column
of the stiffness matrix.
Only the coefficients between the main
diagonal and the skyline are stored.
33
Banded Symmetric Matrices and Bandwidth
nb = nd(m + 1)
m = 4 - 1 = 3 and nd = 2, → nb = 2(3 + 1) = 8.
34
The Execution Time (Primarily, Equation-Solving Time)
→ A function of the number of equations to be solved
Without using Banded Storage of Global Stiffness Matrix K,
the Execution Time is Proportional to (1/3)n3.
n = the Number of Equations to be Solved.
Using Banded Storage of K,
The Execution Time is Proportional to n(nb)2.
The Ratio of Time of Execution without Banded Storage to
that using Banded Storage is (1/3)(n/nb)2.
For the Plane Truss Example, the Ratio (1/3)(24/8)2 = 3.
Therefore, it takes about Three Times as long to execute the
Solution of the Example Truss if Banded Storage is not used.
35
SYMMETRY AND BANDWIDTH
Exercises :
How to sequence node number of the system
36